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In  this  report  a  computational  and  experimental  investigation  of  the  fragmentation  of 
Armco  iron  and  HF-1  steel  explosively  loaded  cylinders  and  an  HF-1  steel  projectile  is  presented. 
The  computational  simulations  were  performed  with  SRI  International  Lagrangian 
finite-difference  computer  programs.1’2 Cylinder  and  projectile  fragment  mass  distributions 
were  experimentally  determined  for  comparison  with  the  computations.  The  PUFF  81  one- 
dimensional  and  TROTT2  two-dimensional  stress  wave  propagation  codes  were  used  for  the 
cylinder  and  projectile  computations,  respectively.  The  SHEAR21,3,4  shear  band  model  and  the 
BFRACT21  brittle  fracture  model  simulated  the  fracture  processes.  These  nuclcation  and  growth 
models  are  subroutines  in  the  stress  wave  propagation  codes. 

The  original  version  of  the  computational  shear  band  model  SHEAR  1 5  was  developed  by 
SRI  as  part  of  a  joint  NSWC-SRI  experimental  and  analytical  effort  to  simulate  the  fragmentation 
behavior  of  naturally  fragmenting  munitions.6  In  that  work  the  PUFF  8  computer  program 
containing  the  SHEAR1  and  BFRACT2  fracture  models  was  used  by  SRI  to  compute  fragment 
mass  distributions  for  cylinders  of  Armco  iron  and  a  selected  heat  treatment  of  HF-1  steel.5  The 
HF-1  steel  brittle  fracture  parmietcrs  that  were  used  as  input  to  the  BFRACT2  fracture  model 
were  determined  at  NSWC  via  gas  gun  impact  and  soft  recovery  experiments.  Fragment  mass 
distributions  for  the  cylinders  were  experimentally  determined  at  NSWC  for  comparison  with  the 
SRI  computations.6  SRI  refined  the  original  shear  band  model  under  joint  sponsorship  of  the 
Army  Ballistic  Research  Laboratory,  the  Advanced  Research  Projects  Agency,  and  NSWC.3  The 
new  shear  band  model  SHEAR2  was  then  used  with  the  PUFF  8  and  TROTT  computer  programs 
to  calculate  fragment  mass  distributions  for  Armco  iron  and  HF-1  steel  cylinders  and  an  HF-1  steel 
projectile,  respectively.3 

In  the  present  work,  computations  were  made  for  a  different  heat  treatment  of  HF-1  steel7 
and  also  for  the  previous  heat-treatment  that  was  used  in  the  earlier  work.  The  two  heat  treatments 
produced  different  microstructures  and  resulted  in  different  fracture  parameters.  This  allowed  the 
predictive  capability  of  the  BFRACT2  model  to  be  tested  for  small  changes  in  the  fracture 
parameters  for  HF-1  steel.  The  fracture  parameters  for  the  previous  heat  treatment  of  HF-1 
steel6  were  redetermined  for  these  calculations.  Computations  were  also  made  for  an  Armco  iron 
cylinder.  The  SHEAR2  and  BFRACT2  fracture  models  were  used  together  in  some  of  the  com- 
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putations,  since  the  recovered  fragments  had  both  shear  and  brittle  fracture  surfaces.  For  those 
cases  where  tire  models  could  not  be  used  together,  the  BFRACT2  model  was  used  to  simulate 
the  fracture  processes  because  the  recovered  fragments  had  more  brittle  than  shear  fracture  surface 
area.  Some  of  the  results  presented  in  the  earlier  report  (Reference  6)  are  presented  in  this  report 
in  more  detail  for  completeness. 

The  materials  are  described  in  Section  II.  Static  tensile-pull  measurements  and  ultrasonic 
velocity  measurements  for  HF-I  steel  are  presented  in  Section  III.  HF-1  steel  Hugoniot  measure¬ 
ments  are  presented  in  Section  IV.  Section  V  describes  the  HF-I  steel  gas  gun  impact  and  soft 
recovery  experiments.  The  calculated  dynamic  fracture  parameters  for  HF-1  steel  that  were  used  as 
input  for  the  cylinder  and  projectile  simulations  are  given  in  Section  VI.  The  results  of  a  framing 
camera  experiment  for  an  explosively  loaded  Armco  iron  cylinder  are  presented  in  Section  VII. 
Section  VIII  contains  the  fragment  mass  distribution  results  of  the  Armco  iron  and  HF-1  steel 
explosively  loaded  cylinder  experiments  and  the  HF-1  steel  projectile  experiment.  A  new 
characterization  scheme  for  grouping  the  recovered  fragments  with  respect  to  their  fracture 
surfaces  is  presented  in  this  section.  The  fragment  mass  distribution  computations  are  presented  in 
Section  IX.  Section  X  contains  the  summary.  Appendixes  A  through  P  contain  tables,  photo¬ 
graphs,  computer  listings,  computer  plots,  and  other  information  concerning  this  work. 


II.  MATERIALS 


The  Armco  iron  material  used  for  the  exploding  cylinder  experiments  was  purchased  from 
Corey  Steel  Company,  Chicago,  Illinois,  as  a  127-mm-diameter  hot-rolled  bar.  It  was  used  in  the 
as-received  condition.  The  chemical  analysis  is  given  in  Table  1.  The  microstructure  of  the  material 
shown  in  Figure  1  consists  primarily  of  a  single-phase  polycrystalline  aggregate  with  ASTM  grain 
size  number  of  approximately  2.5  to  3.  A  few  inclusions  are  also  present;  these  are  probably 
oxides. 


The  HF-1  steel  material  that  was  used  in  this  investigation  was  purchased  from  Norris 
Industries,  Pomona,  California,  in  cylinder  form  and  as  a  127-mm-diaineter  hot-rolled  bar.  The 
chemical  analysis  of  this  silico-manganese  steel  is  given  in  Table  1.  Two  heat  treatments  were 


Table  1.  Chemical  analysis  of  Armco  iron  and  HF-1  steel. 


Element 

Armco  Iron3 
(wt  %) 

HF-1  Steelb 
(wt  %) 

C 

0.004 

1.08 

Mn 

Trace 

1.83 

Si 

Trace 

0.86 

P 

Trace 

0.011 

S 

0.025 

0.0001 

Al 

0.0002 

0.006 

Mo 

Trace 

0.07 

Ni 

0.0001 

Trace 

Cr 

Trace 

0.15 

0 

0.010 

— 

N 

0.001 

- 

Cd 

0.003 

- 

Mg 

0.001 

- 

“Determined  by  Corey  Steel  Company,  Chicago,  Illinois. 
“Determined  by  NSWC,  Dahlgren,  Virginia. 


t  200  nm  | 

Figure  1.  Microstructure  of  Armco  iron. 
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considered  to  investigate  microstructural  changes  on  fragmentation  behavior.  Heat  treatment  A 
material  was  austenized  at  843±4°C  for  4  hr  and  then  quenched  in  warm  oil  at  60°C  to  form 
martensite.  The  martensite  was  tempered  at  566±4°C  for  2  brand  then  air-cooled.  Figure  2  shows 
the  mierostructure  for  this  heat  treatment  at  two  different  magnifications.  It  consists  mainly  of 
tempered  martensite  with  some  dispersed  carbide  precipitates  and  retained  austenite.  Heat  treat¬ 
ment  A  material  was  used  in  the  work  of  Reference  6.  Heat  treatment  B  material  was  austenized  at 
•>27i4°C  for  1  hr.  transferred  to  782±4°C  for  1  hr  to  form  grain  boundary  cementite,  and  then 
an -cooled  to  form  pearlite.7  The  microstructure  for  this  heat  treatment  is  shown  in  Figure  3. 
It  consists  mainly  of  a  cementite  network  and  a  fine  pearlite  structure. 


III.  MECHANICAL  AND  ACOUSTICAL  EXPERIMENTS 

In  this  section  tensile-pull  and  ultrasonic  velocity  measurements  are  presented  for  the 
materials.  This  information  will  be  useful  for  comparison  with  the  shock-wave  measurements  to 
be  presented  in  later  sections.  Tensile-pull  results  are  presented  for  Armco  iron  and  the  two  heat 
treatments  of  HF-1  steel.  Results  are  presented  for  longitudinal  (parallel  to  the  material  rolling 
direction)  and  transverse  (perpendicular  to  the  material  rolling  direction)  tensile  specimens. 
Figures  4  and  5  give  stress-strain  curves  for  longitudinal  tensile  specimens  of  Armco  iron  and  heat 
ticatment  A  of  HF-1  steel,  respectively.  Both  true  and  engineering  stress-strain  curves  are  shown  in 
these  figures.  The  crosses  indicate  specimen  fracture.  Stress-strain  curves  for  transverse  tensile 
specimens  of  heat  treatments  A  and  B  of  HF-1  steel  are  shown  in  Figure  6.  Two  experiments  were 
peilormed  lot  each  HF-1  steel  heat  treatment.  This  figure  was  constructed  from  the  measurements 
listed  in  Table  A-l.  The  specimen  dimensions  and  experimental  details  for  the  transverse  tensile- 
pull  measurements  are  given  in  this  table.  Table  2  summarizes  the  transverse  tensile-pull  measure¬ 
ments  for  the  two  HF-1  steel  heat  treatments.  The  heat  treatment  B  transverse  tensile  specimens 
fractuied  at  a  lower  stress  (about  207(  less)  and  a  lower  strain  (about  70'/  less)  than  the  heat 
ticatment  A  transverse  tensile  specimens.  Photographs  of  the  fracture  surfaces  for  two  specimens 
arc  shown  in  Figure  7.  Table  3  summarizes  the  mechanical  property  measurements  for  the  ma- 
tenals.  File  extent  of  the  directional  character  of  the  mechanical  properties  for  heat  treatment  A 
of  IIF-I  steel  is  shown  in  this  table. 
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Figure  4.  True  stress-true  strain  and  engineering  stress-strair  curves 
for  longitudinal  tensile  specimens  of  Armco  iron. 


Figure  5.  True  stress-true  strain  and  engineering  stress-strain  curves  for  longitudinal 
tensile  specimens  of  HF-1  steel  with  heat  treatment  A. 
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Table  2.  Summary  of  tensile-pull  measurements  fot  transverse-direction  spedr*  is  of  HF-1  steel. 


Experiment 

No. 

Heat 

Treatment8 

Initial 

Specimen 

Diameter 

(mm) 

Young’s 

Modulus1* 

(MPa)c 

0.2%  Offset 
Yield 
Strength 
(MPa) 

Ultimate 

Tensile 

Strengthd 

(MPa) 

Percent 

Elongation6 

Percent 
Reduction 
in  Areaf 

1 

A 

4.70 

19  f 

947 

1090 

2.7 

3.3 

A 

4.71 

191 

1013 

1 120 

4.0 

2.7 

3 

B 

4.70 

192 

855 

885 

- 

- 

4 

B 

4.70 

176 

858 

879 

1.0 

0.8 

d  Specimens  I  and  2  had  a  hardness  of  Rf  39.  Specimens  3  and  4  had  a  hardness  of  Rp  43.  Two  hardness  values  were  measured 
lor  each  disk  prior  (o  fabricating  a  tensile-pull  specimen. 

Obtained  from  the  slope  of  the  el  istic-region  part  of  the  stress-strain  curves  in  Figure  6.  The  value  for  experiment  4  is  about 
8r/t  less  than  the  values  for  the  other  experiments  and  is  therefore  probably  more  uncertain, 
c  I  megapascal  (MPa)  =  145.0  psi.  1  gigapascal  (GPa)  =  1(P  MPa. 

dThc  ultimate  tensile  strength  is  tile  stress  at  fracture.  The  average  strain  at  fracture  for  heat  treatments  A  and  B  is  0.0247  and 
0.0074.  respectively 

u  The  final  central-region  length  for  specimens  I,  2,  and  4  was  19.6,  13.2,  and  12.8  mm,  respectively.  The  final  length  for  specimen 
3  was  no'  measured. 

1  The  linal  central-region  diameter  of  the  necked-down  fracture  region  for  specimens  1,  2,  and  4  was  4.62,  4.64,  and  4.69  mm, 
respectively.  The  final  diameter  for  specimen  3  was  not  measured.  The  fracture  regions  for  specimens  3  and  4  occurred  at  the 
edge  ol  the  central  region  and  m  the  middle  of  the  central  region,  respectively. 
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Figure  6.  Engineering  stress-strain  curves  for  transverse  tensile 
specimens  of  HF-1  steel  with  heat  treatments  A  and  B. 
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(a)  (b) 

Figure  7.  Macrophotographs  of  the  fracture  surfaces  for  the  transverse  tensile  specimens  of 
HF-1  steel  with  (a)  heat  treatment  A  in  experiment  1  and  (b)  heat  treatment  B  in 
experiment  3. 


TV  •V’i 


Table  3.  Summary  of  mechanical  properties  for  Armco  iron  and  HF-1  Steel. 


Material 

Heat 

Treatment 

Hardness 

Tensile-Pull 

Direction3 

0.2%  Offset 
Yield 
Strength 
(MPa) 

Ultimate 

Tensile 

Strength 

(MPa) 

Percent 

Elongation 

Percent 
Reduction 
in  Area 

Armco  iron 

None 

Rb  40 

Longitudinal 

180 

300 

34 

57 

HF-1- Steel 

A 

Rc  40 

Longitudinal 

1040 

1320 

8.4 

14 

A 

Rc  39 

Transverse 

980 

1105 

3.4 

3.0 

IIF-l  Steel 

B 

Rc  43 

Transverse 

857 

882 

1.0 

0.8 

3  The  longitudinal  tensile  properties  for  Armco  iron  and  for  heat  treatment  A  of  HIM  steel  arc  averages  for  three  and  five  measure¬ 
ments,  respectively.  Tile  longitudinal  tensile  measurements  were  performed  on  standard  6.35-mm  ASTM  E-68  tensile  specimens; 
the  percent  elongation  results  are  for  25.4-mm  gage.  The  transverse  tensile  properties  for  heat  treatments  A  and  B  steel 

arc  both  averages  for  two  measurements.  The  details  of  the  transverse  tensile  measurements  arc  given  in  Tables  2  and  A-l. 


Zero-stress  longitudinal  and  shear  wave  velocity  measurements  have  been  performed  for  the 
two  heat  treatments  of  HF-1  steel  using  the  ultrasonic  pulse-echo  technique.  The  longitudinal 
wave  velocity  measurements  were  performed  with  a  6.5-MHz  Dapco  transducer.  A  1.6-MHz 
Panametrics  transducer  was  used  for  the  shear  wave  velocity  measurements.  The  transducers  were 
used  with  an  ultrasonic  pulser/rcceiver  (Panametrics  5050  PR).  A  high-vjscosity  fluid  (Dapco 
couplant  paste)  was  used  to  couple  the  transducers  to  the  specimens.  The  time  between  the  echos 
was  measured  with  an  oscilloscope  (Tektronix  R7704)  and  a  digital  delay  plug-in  unit  (Tektronix 
7D1I).  The  unit  was  checked  for  accuracy  with  a  digital  delay  pulse  generator  (Berkeley 
Nucleonics  7040).  Table  4  summarizes  the  ultrasonic  velocity  measurements  and  the  calculated 
elastic  constants  for  the  two  heat  treatments  of  HF-1  steel. 

Figure  8  shows  ultrasonic  pulses  for  HF-1  steel  specimens  with  heat  treatment  A.To  measure 
the  time  difference  between  echos,  the  first  positive  peak  in  the  pulse  is  aligned  with  a  selected 
graticule  line.  The  digital  delay  unit  is  then  used  to  move  the  pulse  train  across  the  display  in  1-ns 
increments  until  the  same  part  of  the  next  pulse  is  aligned  with  the  fiducial  graticule.  The  wave 
velocity  data  for  each  measurement  was  obtained  by  averaging  the  time  differences  for  five  or 
more  successive  echos. 
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Table  4.  Ultrasonic  velocity  measurements  and  elastic  constants  for  HF-1  steel. 


Longitudinal 

Heat 

Density 

Wave  Velocity' 

Treat¬ 

P°  3 

Cl. 

ment 

(Mg/nr ) 

(km/s) 

A 

7 . 7  7 11 

5.93 's' 

It 

7,;7h 

5.92JJ 

Shear  Wave 

Hulk  Wave 

Velocity  *’ 

Velocity0 

Poisson's 

CS 

(km/s) 

CB 

(km/s) 

R.itioJ 

r 

3  25l,k 

4  59 

0.285 

3.24,lk 

4.58 

0.286 

Young's 

Hulk 

Shear 

Modulus'' 

f 

Modulus 

Modulus® 

K 

K 

(5 

(Ofa) 

(GPa) 

211 

164 

82.1 

210 

164 

81.6 

3  Measured  lor  ;i  nominal  center  frequency  of  6.5  Mil/  of  a  broadband  pulse. 
*’  Measured  tor  a  nominal  center  frequency  of  1 .6  Mil/  ot  a  broadband  pulse. 

CCf}  *  C'2-<4/3  K'| 

d>-  =  KcL^s)2-’ri2<cL/t'sr-2i 

c  i;  =2d  +  i-)c|  p0 

f  K  =  l;<|3(l-2u)| 
s0-  l;'|  2d  +i>)| 

11  Average  for  two  specimens.  I  udi  specimen  v\as  6.36-mm  thick. 

1  Average  of  four  measurements  on  two  specimens. 

1  lire  scatter  in  the  measurements  was  0.5'/  or  less. 
k  Average  ot  sis  measurements  on  two  specimens. 


(a)  (b) 

Figure  8.  Longitudinal  and  shear  echo  pulses  for  HF-1  steel  specimens  with  heat  treatment  A. 

steel,  (a)  Longitudinal  echo  pulses.  The  vertical  scale  is  50  mV/div  and  the  horizontal 
scale  is  0.5  ps/div,  (b)  Shear  echo  pulses.  The  vertical  scale  is  50  mV/div  and  the 
horizontal  scale  is  1  ps/div. 


10 


IV.  HUG0N10T  EQUATION  OF  STATE  EXPERIMENTS 


A  series  of  Hugoniot  experiments  has  been  performed  for  heat  treatments  A  and  B  of  HF-1 
steel.  The  resulting  Hugoniot  equation  of  state  was  used  as  input  to  the  computer  programs  for 
calculating  the  shock  stress  in  HF-1  steel.  The  measurements  were  performed  with  a  40-mm-bore 
diameter  gas  gun.8  Both  direct-impact  and  transit-time  experiments  were  performed.  In  the 
direct-impact  experiment,  an  HF-1  steel  disk  was  impacted  directly  onto  a  quartz  gauge  to 
measure  the  HF-1  steel  shock  stress  and  particle  velocity.  In  the  transit-time  experiments,  an  HF-1 
steel  disk  was  impacted  onto  an  HF-1  steel  specimen  that  was  backed  with  a  quartz  gauge. 
Figure  9  is  a  schematic  of  the  muzzle  region  of  the  gas  gun  for  this  type  of  experiment.  The 
impact  time  and  impact  planarity  arc  measured  with  the  four  charged  tilt  pins  that  are  centered 
around  the  specimen.  The  tilt  pin  ends  are  positioned  in  the  plane  of  the  impact  face  of  the 
specimen  to  within  1  jam.  The  average  projectile  velocity  at  impact  is  measured  with  the  three 
charged  pins  in  the  side  of  the  barrel  and  the  specimen  tilt  pins.  The  quartz  gauge  measures  the 
stress-time  profile  at  the  specimen-quartz  interface.  Elastic  and  plastic  wave  velocities  in  the 
specimen  are  determined  by  measuring  the  time  between  the  impact  of  the  front  surface  of  the 
specimen  and  the  arrival  of  the  stress  waves  at  the  specimen  back  surface.  A  digital  delay  generator 
(Berkeley  Nucleonics  7040)  is  used  for  producing  two  reference  pulses  that  are  time  delayed 
with  respect  to  each  other  by  a  predetermined  amount.  The  initial  pulse  is  recorded  on  the  tilt 
data  trace,  and  the  delayed  pulse  is  recorded  on  the  quartz  gauge  data  trace. 

The  oscilloscope  records  for  the  five  Hugoniot  experiments  are  given  in  Appendix  B.  Shots 
1 17,  1 18,  and  120  were  performed  for  heat  treatment  A  material;  Shots  152  and  153  were  per¬ 
formed  for  heat  treatment  B  material.  All  the  experiments  were  shock  transit-time  experiments 
except  Shot  120,  which  was  a  direct  impact  experiment.  The  tilt  data  record  for  the  transit-time 
experiments  is  a  series  of  voltage  pulses  corresponding  to  the  contact  of  the  tilt  pins  by  the 
impactor  disk  The  tilt  angle  or  impact  planarity  is  measured  from  this  record.  The  quartz  gauge 
current  record  for  each  shot  shows  an  HF-1  steel  elastic  wave  followed  by  a  slower-moving  plastic- 
wave.  A  rising  transition  region  occurs  between  the  elastic  wave  front  and  the  peak  of  the  plastic- 
wave  for  both  HF-1  steel  heat  treatments.  A  similar  phenomenon  has  been  observed  by  Jones 
et  a!.9  for  SAE  4340  steel  and  attributed  to  work  hardening  and  wave  dispersion.  This  is  in 
contrast  to  Armco  iron,  which  exhibits  an  elastic  wave  followed  by  stress  relaxation  and  then  a 
gradual  rise  to  the  plastic  wave.1 0 
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Table  5  summarizes  the  shock  wave  measurements  and  calculations  for  HF-1  steel.  A  brief 
description  of  each  shot  is  given  in  the  third  column  of  the  table.  The  average  measured  density 
for  the  HF-1  steel  disks  is  pQ  =  7.774  Mg/m3.  Impact  tilt  was  measured  from  the  rise  time  of  the 
quartz  gauge  pulse  for  Shot  120.  The  average  tilt  value  for  the  shots  is  0.84  mrad.  Values  for  wave 
velocity,  stress  particle  velocity,  and  strain  are  given  for  the  elastic  and  plastic  waves.  An  average 
value  of  2.2  GPa  was  obtained  for  the  Hugoniot  elastic  limit  (HEL)  for  heat  treatment  A  material. 
The  hardness  for  this  material  is  Rc  40.  An  average  HEL  value  of  2.4  GPa  was  obtained  for  heat 
treatment  B  material  with  a  hardness  of  Rc  43.  Jones  et  al.9  obtained  a  value  of  2.0  GPa  for  the 
HEL  of  SAE  4340  steel  with  a  thickness  of  19.9  mm  and  a  hardness  of  Rc  40.  Minshall1 1  measured 
an  elastic  wave  velocity  of  5.95  km/s  and  an  HEL  of  2.5  GPa  for  25.4-mm-thick  4340  steel  with  a 
hardness  of  Rc  35. 
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Table  5.  Summary  of  Hugoniot  exf 
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\  value  for  the  one-dimensional  static  yield  strength  ay  can  be  calculated  lor  heat 
treatment  A  material  using  the  equation  aY  =1(1  -»))( 1  -  2t>)|  Y  where  v  is  Poisson's  ratio  and  Y 
is  the  longitudinal  yield  strength.  For  this  material  Y=  1.04  GPa  from  Table  3  and  u  =  0.285 
from  Table  4.  The  calculated  value  of  =  1.7  GPa  is  not  as  large  as  the  measured  HFL  value  of 
2.2  GPa  for  heat  treatment  A  material,  indicating  some  strain-rate  dependence  for  HF-1  steel. 
The  ultrasonic  longitudinal  wave  velocity  of  5.93  km/s  from  Table  4  is  in  good  agreement  with  the 
average  elastic  wave  velocity  value  of  5.89  km/s. 

Figure  10  gives  the  shock  velocity-particle  velocity  data  points  for  HF-1  steel.  The  Hugoniot 
relationship  Us  =  4.60  +  1.61  ii  was  obtained  from  a  straight-line  least-squares  fit  to  the  three 
Us-up  plastic  wave  data  points  from  Table  5  and  the  average  value  for  the  ultrasonic  bulk  wave 
velocity  from  Table  4.  The  Hugoniot  curves  for  Armco  iron1 0  and  304  stainless  steel1 2  are  shown 
for  comparison. 


HF-1  STEEL 
Us  =4.60+1.61  IV 
DENSITY  7.77  Mf/m3 
•  HEAT  TREATMENT  A 
▲  HEAT  TREATMENT  B 
■  ULTRASONIC  BULK  \ 
WAVE  VELOCITY _ 


ARMCO  IRON 

U$  =4.63+1.33  UP 
DENSITY  7.85  M|/m3  A 
(REF.  10) 


304  STAINLESS  STEEL 
Us  =4.57+1.49  IV 
DENSITY  7.90  Me/m3 


PARTICLE  VELOCITY  (km/s) 


Figure  10.  Shock  velocity-particle  velocity  relationship  for  HF-1  steel. 
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The  stress-particle  velocity  Hugoniot  relation  and  data  points  for  HF-1  steel  are  shown  in 
Figure  1  I .  A  maximum  stress  of  7. 1 1  GPa  was  obtained  for  these  shots.  The  stress-particle  velocity 
relationship  was  taken  to  be  a  straight  line  up  to  a  particle  velocity  of  u£  =  0.051  km/s,  the 
average  particle  velocity  at  the  elastic  limit.  The  stress  at  this  particle  velocity  is  oc  «  2.3  GPa, 
the  average  HEL  value  for  the  two  heat  treatments  of  HF-1  steel.  Above  a  particle  velocity  of 
0.05 1  km/s  the  stress-particle  velocity  relationship  was  obtained  from  the  equation  oH  =  oe  + 
p.(U  -  u  )  (u- u)  where  p=  p. /(I  -  u/C.  ).  Here  p  and  C,  =  5.89  km/s  are  the  density  and 
wave  velocity  at  the  elastic  limit,  respectively.  Using  the  HF-1  steel  shock  velocity-particle  velocity 
relationship  U$  =  4.60  +  1 .61  up  and  the  measured  values  for  the  other  parameters  gives  oH  =  0.5 1 
+  35.0  u  +  1 2.6  u~  for  up  >  0.05 1  km/s.  Here  oH  is  in  gigapascals. 
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Figure  12  gives  the  Hugoniot  stress  and  pressure  relationships  for  HF-I  steel. 
Iu| nation  (C-22)  of  Appendix  C  was  used  to  transform  the  linear  shock  velocity-particle  velocity 
relationship  into  the  stress-compression  relationship.  The  compression  p  =  Vfl/V  -  1  where  V0 
and  V  are  initial  and  final  specific  volumes,  respectively.  The  stress  relationship  is  valid  above 
the  elastic  limit  compression  /a.  =  0.0087.  In  the  elastic  region  the  stress  is  given  by 
o„  -  270  p/  ( 1  +  p).  At  the  elastic  limit  compression  pe.  the  stress  is  ae  =  2.3  GPa.  The  HF-1 
steel  hydrostat  was  obtained  by  subtracting  the  deviator  stress  o()  =  0.91  GPa  from  the  stress 
relationship.  The  hydrostat  is  part  of  the  input  data  for  the  SRI  computer  programs.  The  Armco 
iron  hydrostat13  that  was  used  in  the  computer  programs  is  shown  in  Figure  12  lor  comparison 
with  the  IIF-I  steel  hydrostat. 


n 


0  0.01  0.02  0.03  0.04 

COMPRESSION  n 


Figure  12.  Stress-compression  and  pressure-compression  relationships  for  HF-1  steel. 
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V.  DYNAMIC  FRACTURE  EXPERIMENTS 


A  series  of  spall  fracture  experiments  has  been  performed  for  heat  treatments  A  and  B  of 
UF-1  steel.  The  purpose  of  the  experiments  was  to  obtain  a  quantitative  measure  of  fracture 
damage  in  this  material  for  use  in  calculating  dynamic  fiacture  parameters.  The  procedure  is  as 
follows.  Differing  amounts  of  spall  fracture  are  produced  in  the  specimens  under  known  impact 
conditions.  The  specimens  are  soft  recovered,  sectioned  on  a  diameter,  and  polished  to  reveal  the 
microscopic  spall  fracture.  The  numbers  of  cracks,  their  lengths,  and  orientations  with  respect  to 
the  impact  plane  are  digitized.  The  surface  distributions  are  converted  to  volume  distributions 
using  the  SRI  statistical  transformation  computer  program  BABS3.14,15 

Figure  13  shows  the  muzzle  region  of  the  gas  gun  for  a  spall  fracture  experiment.  The 
specimen  is  held  in  a  stainless  steel  mounting  ring  with  a  low-melting-point  alloy  (Cerrobend 
alloy).  For  spall  fracture  to  occur,  the  experimental  conditions  must  cause  tensile  stresses  in  the 
specimen.  Figure  14  is  a  simplified  schematic  of  the  spall  process.  The  impactor  thickness  is  less 
than  the  specimen  thickness.  In  the  first  step  after  impact,  a  plane  shock  wave  propagates  into 
|  both  the  impactor  and  specimen.  In  the  next  step,  one  wave  reaches  the  back  surface  of  the 

impactor,  while  the  other  wave  .ontinues  into  the  specimen.  The  wave  in  the  impactor  is  reflected 
as  a  raiefaetion  wave  and  propagates  back  through  the  impactor  toward  the  specimen.  Subseq¬ 
uently,  a  rarefaction  wave  moves  into  the  specimen  from  its  free  surface.  When  the  two  plane 
rarefaction  waves  interact,  the  result  is  a  tensile  stress  in  the  specimen.  The  tensile  stress  occurs 
because  in  rarefaction  waves,  ni:.erial  is  accelerated  m  a  direction  opposite  to  the  rarefaction 
propagation  direction.  If  the  tensile  stress  amplitude  and  duration  are  sufficient,  spall  fracture 
occurs.  The  compressive  stress  amplitude  is  determined  from  the  impact  velocity  and  the  Hugomot 
equation  of  state  of  the  materials  used.  The  stress  duration  is  controlled  by  the  thickness  of  the 
impactor. 
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Figure  13.  Schematic  of  muzzle  region  of  gas  gun  for  a  spall  fracture  experiment. 


BEFORE  IMPACT 


Figure  14.  Simplified  schematic  of  the  spall  fracture  process.  After  impact  the  time  sequence 
of  wave  interactions  is  shown  from  left  to  right.  The  solid  and  dashed  lines  indicate 
the  propagation  of  compression  and  rarefaction  waves,  respectively. 
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When  performing  the  spall  fracture  experiments,  the  impacted  specimen  is  soft  recovered  to 
minimi/.e  any  unintentional  damage.1 6  Figure  15  shows  the  operation  of  the  soft  recovery  system. 
The  system  consists  of  a  steel  box  mounted  on  a  track  and  a  rag  cage.  The  box  is  filled  with 
laminated  sheets  of  plywood.  A  hole  through  the  plywood  core  is  larger  than  the  specimen 
diameter  but  smaller  than  the  diameter  of  the  impactor  and  projectile.  Figure  15(a)  shows  con¬ 
ditions  prior  to  impact.  Figure  15(b)  shows  conditions  shortly  after  impact.  The  projectile  and 
impactor  are  caught  in  the  separator  box.  and  the  box  moves  along  the  track  as  momentum  is 
transferred  to  it.  The  specimen  passes  unhindered  through  the  separator  box  and  moves  toward  the 
rag  cage.  In  Figure  15(c)  the  box  is  stopped  by  shock  absorbers  at  the  end  of  the  track,  and  the 
specimen  is  soft  recovered  in  the  rags.  The  recovered  specimen  is  sectioned  along  a  diameter, 
mounted,  and  metallographically  polished  to  reveal  the  crack  distribution  as  schematically  shown 
in  Figure  16. 


BARREL  SPECIMEN  SEPARATOR  RAGS 


Figure  15.  Schematic  of  soft  recovery  technique,  (a)  Before  impact,  (b)  In-flight  capture  of 
projectile  and  impactor,  (c)  Soft  recovery  of  specimen. 
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Figure  16.  Schematic  of  sectioned  specimen  to  reveal  fracture  damage. 


A  summary  of  the  spall  fracture  experiments  performed  for  HF-1  steel  is  presented  in 
Table  6.  Twelve  experiments  were  performed  for  heat  treatment  A  material,  and  eleven  experi¬ 
ments  were  performed  for  heat  treatment  B  material.  The  impaetor  disks  were  nominally  1.1 5-, 
1.59-.  and  2.37-mm  thick  The  specimen  disks  were  nominally  3.18-  and  6.35-mm  thick.  The 
impaetor  velocities  ranged  from  0.120  to  0.276  km/s.  The  last  column  of  the  table  gives  a 
qualitative  degree  of  f'actnre  damage  for  the  specimens  as  determined  by  microscopic  observation. 
Typical  comoressive  stress  amplitudes  ranged  from  about  2  to  5  GPa.  and  typical  stress  durations 
were  fractions  of  a  mic.oxeeond. 

To  obtain  toe  crack  distributions  the  sectioned  specimens  are  lirst  macrophotographed  at 
about  5X  magnification  using  a  Bausch  and  Lomb  L  camera.  During  this  procedure  it  is  very 
important  to  illuminate  the  specimen  with  diffuse  light  so  that  all  the  cracks  are  visible  in  the 
photograph.  Appendix  D  contains  these  photographs  for  all  the  sectioned  specimens  of  HF-I  steel. 
Figures  17  and  18  show  the  effect  of  increasing  impaetor  velocity  on  the  specimen  fracture 
damage  for  heat  treatments  A  and  B  of  HF-1  steel,  respectively.  For  these  shots  the  impaetor  and 
specimen  thicknesses  were  2.37  and  6.35  mm.  respectively.  The  specimens  in  Figures  17(a)  and 
18(b)  were  impacted  at  the  same  velocity.  The  increased  fracture  damage  shown  in  Figure  18(b). 
compared  with  Figure  17(a),  suggests  that  the  fracture  threshold  stress  is  less  for  heat  treatment  B 
material  than  heat  treatment  A  material. 
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Table  6.  Summary  of  spall  fracture  experiments  for  HF-1  steel. 
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For  crack  counting,  each  5X  photograph  was  enlarged  again  and  arranged  together  to  form 
a  24X  composite  photograph  of  the  specimen  surface.  Cracks  were  marked  as  straight-line  segments 
on  a  clear  plastic  sheet  that  was  placed  over  the  photograph.  Using  this  procedure  the  original  crack¬ 
counting  photograph  was  not  marked.  Only  cracks  in  the  approximate  central  one-thi:d  region  of 
the  specimen  were  counted.  This  was  done  to  ensure  that  the  measured  cracks  were  not  influenced 
by  stress  relief  wave  effects  from  the  edge  of  the  specimen.  The  general  crack  counting  criterion 
used  was  to  count  the  longest  segment  as  a  crack  and  each  intersecting  segment  as  a  crack.  If  the 
intersecting  crack  extended  through  the  primary  crack,  its  length  was  taken  as  extending  from  tip 
to  tip.  If  the  intersecting  crack  stopped  at  the  primary  crack  or  appeared  to  extend  through  the 
primary  crack  with  a  relative  displacement,  each  segment  was  treated  as  an  individual  crack.  Table 
7  summarizes  the  crack  measurements  for  the  spall  fracture  specimens  of  HF-1  steel.  The  specimen 
magnification  factors,  region  of  damage,  and  number  of  measured  cracks  are  lisied  for  each  shot. 

The  cracks  for  each  specimen  were  digitize!1  by  mounting  the  plastic  sheet  with  the  indicated 
cracks  on  a  digitizer  table  (Computer  Equipment  Corporation  Model  TF-10C  Digitizer)  to 
deteimine  the  coordinates  of  the  cracks.  The  digitized  crack  data  was  then  used  as  input  to  a 
computer  program  that  divided  the  damage  region  into  a  series  of  zones  parallel  to  the  specimen 
faces.  The  average  zone  width  was  0.40  mm.  The  length,  orientation,  and  location  of  the  midpoint 
of  each  crack  relative  to  the  impact  face  was  calculated.  Each  crack  was  assigned  to  a  zone  based 
on  its  midpoint.  This  computer  program  was  written  to  interface  with  the  SRI  BABS3  program 
that  calculated  surface  and  volume  crack-size  distributions  for  each  damage  zone. 

The  parameters  Ny  and  R,  are  obtained  for  each  damage  zone  containing  cracks  by  fitting 
the  equation  Ng  =  N0  exp(-R/Rj )  to  the  volume  crack-size  distribution.  Here  Ng  is  the  number 
of  cracks  per  unit  volume  with  radius  greater  than  R,  NQ  is  the  total  number  of  cracks  per  unit 
volume,  and  R(  is  the  average  crack  radius.  The  zone  of  maximum  damage  is  taken  as  that  zone 
in  which  NQR^  is  a  maximum.  The  NQ  and  Rj  values  for  the  zone  of  maximum  damage  are  used 
for  each  specimen  to  calculate  the  dynamic  fracture  parameters.  Appendix  E  contains  computer¬ 
generated  plots  for  each  sectioned  specimen  of  HF-1  steel  (excluding  the  no-damage  Shots  89  and 
147).  Each  figure  shows  plots  of  the  digitized  cracks,  and  the  computed  surface  and  volume  crack- 
size  distributions  for  each  damage  zone.  Appendix  F  lists  the  digitized  crack  coordinates  for  all 
the  counted  cracks.  Table  8  summarizes  the  N()  and  R(  values  for  each  spall  fracture  specimen  of 
HF-1  steel.  The  N0  and  Rj  values  determined  for  Shots  95  and  88  were  not  used  for  calculating 
dynamic  fracture  parameters. 
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Table  8.  Summary  of  experimental  volume  crack-size  distribution  parameters  for  spall  fracture 
specimens  of  HF-1  steel.* 


Distance  from 


Shot 

No. 

Heat 

Treatment 

Number  of 

Cracks  per 

Unit  Volume 

4N°  3 
( 1 04  /cm3 ) 

Average 

Crack 

Radius 

R. 

(pm) 

Specimen  Impact 
Surface  to  Midpoint 
of  Zone  of 
Maximum  Damage 
(mm) 

95b 

A 

0.03 

-1822.2 

2.73 

110 

A 

0.74 

150.3 

4.04 

98 

A 

2.55 

121.1 

3.96 

92 

A 

2.89 

103.9 

4.40 

91 

A 

8.49 

54.2 

4.00 

90 

A 

4.96 

99.7 

4.19 

94 

A 

7.37 

94.2 

4.17 

89c 

A 

- 

- 

- 

113 

A 

7.92 

66.6 

1.84 

88d 

A 

3.92 

107.6 

1.51 

111 

A 

5.36 

75.5 

2.19 

112 

A 

7.48 

84.7 

2.09 

154 

B 

1.45 

200.3 

3.15 

145 

B 

1.24 

190.5 

3.73 

144 

B 

4.26 

127.0 

3.36 

143 

B 

7.10 

88.7 

4.58 

142 

B 

3.68 

132.5 

4.32 

150 

B 

4.25 

114.9 

4.22 

148 

B 

4.58 

99.1 

2.01 

149 

B 

6.64 

97.7 

1.58 

147c 

B 

- 

- 

— 

146 

B 

9.64 

99.3 

2.03 

151 

B 

6.32 

115.8 

4.51 

3  The  parameters  Ng  and  R  j  were  obtained  for  each  damage  zone  of  each  specimen  by  fitting  the  equation  =  Ng  exp  (~R/Rj) 
to  the  volume  crack-stzc  distribution.  Here  N  is  the  number  of  cracks  per  unit  volume  with  radius  greater  than  R.  Plots  of  the 
volume  crack-size  distributions  for  each  damage  zone  of  each  specimen  are  given  in  Appendix  E.  The  parameters  Ng  and  Rj  listed 
in  this  table  are  the  values  for  the  zone  of  nia.\i:>tuiii  damage. 

b This  shot  was  not  used  for  determining  dynamic  fracture  parameters  due  to  the  small  number  of  measured  cracks,  hence,  the 
unrealistic  negative  value  for  the  average  crack  radius. 

c  Cracks  were  not  measured  for  the  no-damage  Shots  89  and  147. 

dThis  full-spall  shot  was  not  used  for  determining  dynamic  fracture  parameters.  Only  cracks  in  one  half  ol  the  specimen  were 
digitized. 
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Scanning  electron  microscope  (SHM)  photographs  of  the  fracture  surface  for  the  lull-spalled 
heat  treatment  A  specimen  in  Shot  88  are  shown  in  Figure  19.  These  micrographs  can  be  com¬ 
pared  with  the  fracture  surface  for  the  low-strain-rate  heat  treatment  A  transverse  tensile  specimen 
m  Figure  7(a).  Micrographs  of  cracks  were  made  for  some  of  the  sectioned  specimen  surfaces  in 
order  to  obtain  a  detailed  comparison  of  the  spall  fracture  for  the  two  heat  treatments  of  HI'-l 
steel.  The  specimens  were  etched  with  either  2'/<  Nital  or  4'/  Picral  solutions.  Figures  20  and  2 1  are 
crack  micrographs  for  the  heat  treatment  A  specimens  in  Shots  1 10  and  90.  respectively.  Figures 
22  and  23  are  crack  micrographs  for  the  heat  treatment  B  specimens  in  Shots  144  and  142.  re¬ 
spectively.  The  impactor  velocity  for  Shots  1 10  and  144  was  0.149  km/s.  The  impactor  velocity  for 
Shots  90  and  142  was  0.191  km/s.  Cracks  occurred  primarily  through  the  grain  boundary 
cementite  for  heat  treatment  B  material,  although  some  cracks  occurred  in  the  pearlite  micro- 
structure. 


Figure  19.  SEM  micrographs  of  the  fracture  surface  for  the  full-spalled  HF-1  steel,  heat  treatment  A, 
specimen  in  Shot  88. 

C 
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Figure  21.  Micrographs  of  cracks  for  the  HF-1  steel,  heat  treatment  A,  specimen  in  Shot  90. 

(a)  View  of  multiple  cracks,  (b)  High-magnification  view  of  a  crack,  (c)  High-mag- 
nificatiori  view  of  cracks  around  material  impurity. 
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Figure  23.  Micrographs  of  cracks  for  the  HF-1  steel,  heat  treatment  B,  specimen  in  Shot  142. 

(a)  View  of  multiple  cracks,  (b)  High-magnification  view  of  grain  boundary  cementite 
cracks  extending  into  an  adjacent  grain,  (c)  High-magnification  view  of  grain  boundary 
crack. 
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VI.  COMPUTATION  OF  DYNAMIC  FRACTURE  PARAMETERS 
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In  this  section  the  dynamic  fracture  parameters  are  calculated  for  heat  treatments  A  and  B 
of  HF-1  steel.  Initially,  PUFF  8  computations  assuming  no  fracture  damage  were  made  for  the 
twelve  heat  treatment  A  and  the  eleven  heat  treatment  B  spall  fracture  experiments  of  the  last 
section.  The  material  properties  determined  in  the  previous  sections  were  used  as  input  in  the 
computations.  For  a  computation,  the  impactor  thickness  was  divided  into  8  computational 
cells,  numbered  1  through  8,  beginning  at  the  impactor  free  surface.  The  specimen  thickness  was 
divided  into  20  computational  cells,  numbered  10  through  29,  beginning  at  the  specimen  impact 
surface.  The  initial  -maximum  tensile  stress  o  for  each  specimen  was  determined  from  these 
computations.  (Compressive  stresses  have  positive  values;  tensile  stresses  have  negative  values.) 
A  maximum  tensile  stress  of  -2.5  GPa  was  obtained  for  the  heat  treatment  A  specimen  in  Shot  95. 
This  value  was  takemas  the  nucleation  threshold  stress  ano  for  heat  treatment  A  material,  since  the 
experimental  crack  density  for  this  shot  was  only  300  cracks/cm3  (see  Table  8),  a  value  con¬ 
sidered  negligible.1 7  For  the  no-damage  Shot  89,  a  maximum  tensile  stress  of -2.4  GPa  was 
calculated,  a  value  less  than  the  nucleation  threshold  stress.  For  heat  treatment  B  material,  a  shot 
near  the  threshold  stress  was  not  obtained  for  the  2.37-mm-thick  impactor  disks.  All  the  experi¬ 
mental  crack  density  values  in  Table  8  for  this  heat  treatment  are  larger  than  1000  cracks/cm3 , 
except  Shot  147,  which  was  a  no-damage  shot  fora  1 . 1 5-mm-thick  impactor.  In  Shot  154  (the 
lowest  velocity  shot)  a  6.35-mm-thick  specimen  was  impacted  at  0.094  km/s  with  a  2.37-mm-thick 
impactor.  The  maximum  tensile  stress  for  this  shot  was -2.1  GPa,  The  nucleation  threshold  stress 
for  heat  treatment  B  material  was  estimated  by  plotting  the  experimental  crack  density  values 
from  Table  8  versusithe  maximum  tensile  stress  values  for  the  no-damage  PUFF  8  computations. 
The  stress-axis  intercept  of  this  plot  was  -1.8  GPa.  This  value  was  therefore  taken  as4he  nucle¬ 
ation  threshold  stress  for  heat  treatment  -'immaterial.  The  growth  threshold  stress  ogo  for  both  heat 
treatments  of  HF-1  steel  was  taken  as  -0:1  GPa.6  The  nucleation  and  growth  threshold  stress 
values,  ano  and  ago,  respectively,  are  two  of  the  six  dynamic  fracture  parameters.  To  obtain  initial 
estimates  for  the  remaining  parameters,  it  is  necessary  to  obtain  estimates  for  the  spall  fracture 
nucleation  and  growth  times,  Atn  and  Atg,  respectively,  for  each-of  the  shots  for  the  zone  of 
maximum  stress.  The  nucleation  time  is  the  time  interval  for  which  the  initial  tensilestress  pulse 
exceeds  the  nucleation  threshold  stress.  The  growth  time  is  the  time  interval  for  which  the  initial 
tensile  stress  pulse  exceeds  the  nucleation  threshold  stress  until  it  decreases  to  the  growth 
threshold  stress. 
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The  nucleation  rate  for  cracks  N  and  the  growth  rate  for  cracks  R  are  given,  respectively, 


(1) 


and 


R  =  Ti(o-  aKo>R  ,  (2) 

where  a  is  the  applied  tensile  stress,  N0  is  the  nucleation  threshold  rate,  Oj  is  the  nucleation 
sensitivity  parameter,  R  is  the  crack  radius,  and  T,  is  the  growth  coefficient.  Integrating 
Equation  (2)  gives 

/ 

InR  =  lnR0 +T,ot  ,  (3) 

where  RQ  is  the  nucleation  size  parameter.  Here  the  growth  threshold  stress  ago  is  assumed  to  be 
negligible.  N0 ,  a, ,  R0 ,  and  T,  arc  the  remaining  foe.*  dynamic  fracture  parameters. 

A  firsticstimatc  for  the  nucleation  parameters  is  obtained  by  plotting  the  values  ln(N0/Atn) 
versus  the  stress  above  threshold  0  "  ano  f°r  each  of  the  shots  for  eachiheat  treatment.  The 

NQ  values  were  obtained  from  Table  8,  and  the  o|n  and  Atn  values  were  obtained  from  the  initial 
no-damage  PUFF  computations.  A  straight-line  least-squares  fit  to  the  nucleation  data  gives  initial 
estimates  for  the  intercept  N0  and  slope  o" 1 .  A  first  estimate  for  the  growth  parameters  is  obtained 
by  plotting  the  values  lnRj  versus  am  Atg  for  each  of  the  shots  for  each  heat  treatment.  The  R, 
values  were  obtained  from  Table  8,  and  the  Atg  values  were  obtained  from  the  initial  no-damage 
PUFF  computations.  A  straight-line  least-squares  fit  to  the  growth  data  gives  initial  estimates 
for  the  intercept  R0  and  the  slope  T, .  Data  for  Shots  95,  89,  88,  and  147  was  not  used  in  these 
estimates  or  any  later  computations. 
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Tlie  parameters  c/no  and  og0  and  the  initial  estimates  for  the  parameters  N0,  Oj ,  R0,  and  Tj 
are  then  used  in  a  first  iteration  of  PUFF  8  computations  using  the  BFRACT2  subroutine,  which 
allows  fracture  to  occur.  In  these  computations,  crack  densities  N0  and  radii  Rj  are  computed 
for  each  specimen  cell  in  which  fracture  occurs.  The  cell  of  maximum  damage  for  a  specimen  is 
taken  as  the  cell  in  which  the  computed  value  NQ  R3  is  a  maximum  at  a  time  when  the  tensile  pulse 
has  decreased  to  the  --0.1  GPa  threshold  growth  stress.  The  cell  of  maximum  damage  is  not 
necessarily  the  same  as  the  cell  of  maximum  stress  for  the  initial  no-damage  PUFF  computations. 
Using  the  results  of  these  computations,  nucleation  and  growth  plots  are  again  constructed. 
Straight-line  least-squares  fits  to  the  experimental  and  computational  data  are  obtained  for  each 
plot.  The  NQ  and  Rj  values  from  Table  8  were  again  used  to  construct  the  ln(NQ/Atn)  and  lnRj 
values,  respectively  for  the  experimental  curves.  For  the  computational  curves,  the  computed  NQ 
and  Rj  values  were  used.  A  table  is  then  constructed  consisting  of  the  initial  set  of  dynamic 
fracture  parameters  Nfl,  ,  Rfl,  and  Tj  and  the  experimental  and  computed  intercept  and  slope 
values  from  the  nucleation  and  growth  plots.  Ratios  of  the  corresponding  experimental  and 
computed  intercept  and  slope  values  are  determined.  A  new  set  of  fracture  parameters  is  obtained 
by  multiplying  the  appropriate  parameter  by  the  appropriate  ratio.  The  new  set  of  dynamic 
fracture  parameters  is  then  used  in  a  second  series  of  PUFF  8  computations  using  the  BFRACT2 
subroutine.  Nucleation  and  growth  plots  are  again  obtained  and  then  a  new  set  of  dynamic 
fracture  parameters.  This  iterative  process  is  continued  until  the  experimental  and  computed 
nucleation  and  growth  curves  agree.  The  dynamic  fracture  parameters  used  for  this  final  iteration 
are  themtaken  as  the  dynamic  fracture  parameters  for  the  material. 

Six  iterations  with  the  PUFF  8  computer  program  using  the  BFRACT2  subroutine  were 
performed  to  obtain  the  dynamic  fracture  parameters  for  heat  treatments  A  and  B  of  HF-1  steel. 
Tables  9  and  10  give  the  experimental  and  computed  nucleation- and  growth  values  that  were 
obtained  using  the  final  set  of  dynamic  fracture  parameters  for  heat  treatments  A  and  B, 
respectively.  These  values  are  plotted  in  Figures  24  and  25.  Straight-line  least-squares  fits  to 
the  experimental  and  computed  valuescare  shown  in  the  figures.  Theexperimental  and  computed 
intercept  and  slope  values  for  the  heat  treatment  A  nucleation  rate  curves  in  Figure  24(a)  are  15.1 
and  16.9  cracks/cm3 -ns  and  - 12.9  and  - 12.0  GPa"1 ,  respectively.  The  experimental  and  computed 
intercept  and  slope  values  for  the  heat  treatment  A  growth  curves  in;  Figure  24(b)  are  68.3  and 
68.4  um  and  -0.137  and  -0.141  GPa-1  us-1,  respectively.  The  experimental  and  computed 
intercept  and  slope  values  for  the  heat  treatment  B  nucleation  rate  curves  in  Figure  25(a)  are  12.3 
and  12.5  cracks/cm3-ns  and  -5.56  and  -5.54  GPa"1,  respectively.  The  experimental  and 
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Summary  of  experimental  and  computed  nucleation  and  growth  values  for  spall  fracture  specimens  of  HF-1  steel  with  heat  treatment  B 
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Figure  24.  Experimental  and  computed  nucleation  rate  and  growth  curves  for  HF*1  steel  with 
heat  treatment  A.  (a)  Nucleation  rate  curves,  (b)  Growth  curves.  The  points  in  these 
plots  are  from  Table  9  and  were  obtained  using  the  final  set  of  dynamic  fracture 
parameters  for-  this  heat  treatment  of  HF-1  steel  in  a  PUFF  8  computation. 


computed  intercept  and  slope  values  for  the  heat  treatment  B  growth  curves  in  Figure  25(b) 
are  82.7  and  82.8  jam  and  -0.148  and  -0.152  GPa" 1  ps" 1 ,  respectively.  For  the  nucleation  rate 
curves,  all  the  experimental  points  were  used  in  the  curve  fits,  since  considerable  scatter  was 
evident  in  the  data.  For  the  growth  curves,  the  experimental  points  for  Shots  1 10,  91,  and  1 13  for 
heat  treatment  A  material  and  the  points  for  Shots  154,  145,  and  143  for  heat  treatment  B 
material  were  not  used  in  the  curve  fits,  since  they  were  out  of  range  of  the  majority  of  the 
experimental  points. 

Table  1 1  summarizes  the  dynamic  fracture  parameters  for  heat  treatments  A  and  B  of  HF-1 
steel.  Also  shown  for  comparison  are  the  previously  determined  dynamic  fracture  parameters 
for  Armco  iron.1  For  the  two  heat  treatments  of  HF-1  steel,  larger  differences  occurred  in  the 
nucleation  parameters  than  in  the  growth  parameters.  The  nucleation  threshold  stress  is  about 
40%  higher  and  the  nucleation  threshold  rate  is  about  30%  lower  for  heat  treatment  A  material 
than  for  heat  treatment  B  material.  This  suggests  that  heat  treatment  A  material  is  a  somewhat 
tougher  material  with  respect  to  crack  nucleation  than  heat  treatment  B  material. 

Table  11.  Dynamic  fracture  parameters  for  heat  treatments  A  and  B  of  HF-1  steel  and  Armco  iron. 


Parameter3 

Description 

Units 

Heat 

Treatment  A 

of  HF-1  Steel 

Heat 

Treatment  B 

of  HF-1  Steel 

Armco  Ironb 

Tl 

Growth  coefficient 

GPa" 1  fis~ 1 

-0.377 

-0.387 

-5.5 

atw 

Growth  threshold  stress 

GPa 

-0.1 

-0.1 

-0.02 

Ro 

Nucleation  size  parameter 

jum 

65.6 

81.8 

0.5 

No 

Nucleation  threshold  rate 

no.  cm-3  ns-1 

90.0 

125 

4600 

ano 

Nucleation  threshold  stress 

GPa 

-2.5 

-1.8 

-0.3 

°1 

Nucleation  sensitivity  parameter  GPa 

-0.0588 

-0.174 

-0.53 

a  In  llic  IIFRACT2  brittle  fracture  subroutine,  these  parameters  are  listed  as  T1  ,T2,T3,T4,  T5,  andT6,  respectively,  with  different 
units.  For  example,  for  heat  treatment  A  material,  the  values  arc  T1  =  -  3.77  .\  10-3  cin^/dync-s,  T2  =  -1.0  x  10^  dync/cni”, 
T3  =  6.56x  10-3  cm,T4  =  9.00  x  10 10  cm-3  s-1,T5  =  -2.50 x  101®  dync/cm2,  and  T6  =  -5.88  x  It)8  dyne/cm2. 
b  These  dynamic  fracture  parameters  for  Armco  iron  were  taken  from  Reference  1  and  are  shown  here  for  comparison  with  the 
IIF-I  steel  values. 


Input  listings  for  the  PUFF  8  computations  using  the  final  set  of  dynamic  fracture  para¬ 
meters  are  given  in  Appendix  G  for  the  HF-1  steel  spall  fracture  shots.  A  detailed  explanation  of 
the  meaning  of  each  of  the  input  parameters  can  be  found  in  Reference  1.  A  brief  description  is 
as  follows:  The  impactor  velocity  is  given  by  UZERO.  RHOS  =  7.774  g/cm3  is  the  average 
impactor  and  specimen  density.  The  EQST  parameters  (EQSTC,  EQSTD  EQSTE,  EQSTG, 
EQSTH,  and  EQSTS  from  left  to  right  in  the  EQST  line)  give  the  equation  of  state  parameters. 
The  parameters  EQSTC,  EQSTD,  and  EQSTS  represent  C,  D,  and  S,  respectively,  in  the  Hugoniot 
pressure-compression  relation  P{1  =  C/a  +  Dp2  +  S/a3  where  /a  =  VQ/V- 1.  VQ  and  V  are  the  initial 
and  final  specific  volumes,  respectively.  For  HF-1  steel  (see  Figure  12),  EQSTC  =  1.605  x  1012 
dyne/cm2,  EQSTD  =  3.561  x  101'2 dyne/cm2,  and  EQSTS  =  3.683  x  101 2  dyne/cm2.  The 
parameters  EQSTE,  EQSTG,  and  EQSTH  are  the  sublimation  energy,  Gruneisen  parameter,  and 
the  vapor  expansion  coefficient,  respectively.  The  sublimation  energy  was  taken  as  7.36  x  1010 
erg/g.1  ’* 8  The  Gruneisen  parameter  70  is  calculated  from  the  equation  yQ  ~  PCg/Cp  where  0  is  the 
volume  thermal  expansion  coefficient,  Cp  is  the  specific  heat  at  constant  pressure,  and  CQ  is  the 
bulk  sound  speed.  For  iron  at  25°C,  the  linear  thermal  expansion  coefficient  is  1 1.7  x  10-6  °K-1 
and  the  specific  heat  at  constant  pressure  is  0.448  x  107  erg/g-°K.1?  From  Figure  10,  CQ  =  4.60 
km/s  for  HF-1  steel.  Using  these  values  the  Gruneisen  parametter  yg  is  1.66.  The  parameter 
EQSTH  =  0.25  for  metals.1  The  TSR1  and  TSR2  lines  give  the  dynamic  fracture  parameters  (see 
Table  11)  and  other  input  parameters  for  the  BFRACT2  subroutine.1  The  two  YO  values 
represent  the  HF-1  steel  yield  strength  and  shear  modulus.  The  yield  strength  Y  was  calculated 
from  the  deviator  stress  oD  =  0.907  GPa  (see  Figure  12)  using  the  equation  Y  =  (3/2 )aD.  This 
gives  Y  =  1.361  x  101 0  dyn/cm2.  The  shear  modules  G  is  obtained  from  the  equation  aD  =  (4/3)G 
Me/(l  +  PC)  where  pe  =  0.00867  is  the  HF-1  steel  elastic  limit  compression.  Using  this  and  the 
value  for  pD  gives  G  =  7.912  x  101 1  dyne/cm2.  The  number  of  computational  cells  and  thickness 
for  an  impactor  and  specimen  are  given  in  the  last  two  lines  of  a  PUFF  listing,  respectively. 

Appendix  H  gives  computed  stress,  particle  velocity,  and  crack  concentration  versus  time 
plots  for  the  cell  of  maximum  damage  for  the  HF-1  steel  spall  fracture  shots.  Appendix  I  gives 
computed  stress,  particle  velocity,  andscrack  concentration  versus  time  plots  for  each  impactor 
and  specimen  computational  cell  for  thesHF-1  steel,  heat  treatment  B,  Shot  144.  This  appendix  is 
included  so  that  selected  shock  histories  can  be  observed  for  all  computational  cells  for  a  typical 
shot.  Figures  26  and  27  give  propagation-direction  stress  histories  and  crack  concentration 
histories,  respectively,  for  selected  computational  cells  for  this  shot.  These  figures  were  con¬ 
structed  from  some  of  the  plots  in  Appendix  I.  They  indicate  in  a  graphic  manner  the  time 
evolution  of  these  selected  properties. 
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Figure  26.  Computed  stress-time  profiles  for  selected  computational  cells  for  the  HF-1  steel, 
heat  treatment  B,  specimen  in  Shot  144.  The  stress  is  parallel  to  the  shock 
propagation  direction.  The  cell  numbers  are  indicated  on  the  profiles.  This  figure 
was  constructed  from  plots  in  Appendix  I. 


Figure  27.  Computed  crack  concentration-time:  prof iles  for  those  computational  cells  in  which 
tensile  fracture  occurred  for  the  HF-1  steel,  heat  treatment  B,  specimen  in  Shot  144. 
The  cell  numbers  are  indicated  on  the  profiles.  This  figure  was  constructed  from 
plots  in  Appendix  I. 
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In  Tabic  7  of  Section  V,  it  was  pointed  out  that  for  some  of  the  shots  for  both  heat  treat¬ 
ments  of  HF-1  steel  some  cracks  occurred  along  the  stringer  impurities  in  the  material.  A  total  of 
44  stringer  cracks  were  counted  for  heat  treatment  A  material.  This  corresponds  to  about  5%  of 
the  total  number  of  digitized  cracks  for  this  heat  treatment.  The  stringer  cracks  were  in  most  cases 
parallel  to  the  shock  wave  propagation  direction.  We  were  interested-in  determining  the  signifi¬ 
cance  of  the  stringer  cracks  on  the  computed  dynamic  fracture  parameters.  The  fracture  para¬ 
meters  for  heat  treatment  A  material  with  stringer  cracks  excluded  were  therefore  determined  for 
comparison  with  the  parameters  in  Table  1 1.  Appendix  J  contains  the  tables  and  figures  for  this 
computation.  With  the  exception  of  an  appreciable  decrease  in  the  nucleation  threshold  rate  N0, 
the  exclusion  of  the  stringer  cracks  did  not  significantly  affect  the  dynamic  fracture  parameters  for 
heat  treatment  A  material. 


VII.  ARMCO  IRON  FRAMING  CAMERA  EXPERIMENT 


In  this  section,  the  results  of  a  framing  camera  experiment  for  an  explosive-filled  Armco 
iron  cylinder  are  presented.  The  cylinder  had  inside  and  outside  diameters  of  76.2  and  1 14.3  mm, 
respectively,  and  was  381  mm  long.  It  was  filled? with  4.10  kg  of  cast-in-place  composition  Bex- 
plosive  that  extended  beyond  the  cylinder  to  minimize  end  effects.  Figure  28  shows  the 
explosive-filled  cylinder  prior  to  detonating  the  explosive  on  the  right  end.  The  explosive  extended 
64  and  90  mm  beyond; the  cylinder  on  the  right:  and  left  ends,  respectively.  Large  and  small  grid 
patterns  were  placed  on  the  cylinder  surface  in  an  attempt  to  measure  the  strain  history  ofsthe 
r  ending  cylinder  using  the  Moire  interference  technique.  This  procedure  (described  in  detail  in 
rences  5  and  6)  was  not  successful  because  there  was  insufficient  detail  in  the  film  to  resolve 
the  interference  patterns.  The  Moire  grid  patterns,  which  extended  191  mm  along  the  cylinder 
axis,  can  be  seen  in  the  middle  of  the  cylinder  in-Figure  28.  The  gross  strain  history  of  the  cylinder 
was  obtained  by  measuring  its  diameter  in  =the;framing  camera  photographs.  A  square  grid  pattern 
wjth  a  25.4-mm  spacingiwas  placed  behind  the  cylinder  for  reference. 
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Figure  28.  Framing  camera  photograph  of  Armco  iron  cylinder  prior  to  detonation  of  the 
composition  B  explosive. 

The  Armco  iron  framing  camera  photographs  are  given  in  Appendix  K.  The  cylinder  outer 
radius  was  measured  at  a  point  near  its  mjd-length  on  21  photographs.  Table  12  summarizes 
the  measurements.  Figure  29  is  a  plot  of  these  values  as  a  function  of  time.  The  slope  of  the  radius¬ 
time  plot  increased  during  the  initial  stages  oficylinder  expansion  (from  about  25  to  65  ps).  During 
the  final  stages  of  cylinder  expansion  (from  about  65  to  95  ps),  the  slope  was  essentially  constant 
at  a  value  of  about  1  km/s.  The  decreased  cylinder  acceleration  that  occurred  during  the-final 
stages  of  expansion  may  be  due  to  explosive  relief  wave  effects  from  the  ends  of  the  cylinder  and 
possibly  around  the  partially  formed  fragments. 


TIME  AFTER  EXPLOSIVE  DETONATION  FRONT  ENTERS  CVUNDCR  M 

Figure  29.  Armco  iron  cylinder  outer  radius  versus  time  plot.  The  data  points  were  obtained 
from  Table  12.  Visible  cracks  first  appeared  on  the  cylinder  surface  at  a  radius  of 
about  65  mm  (14%  increase  in  radius).  This  corresponds  to  about  37  ps  after  the 
detonation  front  entered  the  cylinder  and  about  14  ps  after  the  detonation  front 
passed'the  position  where  the  outer  radius  was  measured: 
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Table  12.  Armco  iron  cylinder  outer  radius  versus  time  values  near  the  mid-length  position  of 
the  cylinder. 


Frame 

No.* 

Time  After 

Explosive  Detonation 
Front  Enters  Cylinder* 

(Ms) 

Axial  Position  of 
Explosive  Detonation 

Front  in  Cylinder* 

(mm) 

Cylinder 

Outer  Radius’* 
(mm) 

1 

6 

47 

57.2 

2 

10 

80 

57.2 

3 

14 

114 

57.2 

4 

18 

147 

57.2 

5 

23 

180 

57.2 

6 

27 

214 

58.9 

7 

31 

247 

60.3 

8 

35 

281 

63.5 

9 

39 

314 

65.7 

10 

43 

348 

69.4 

11 

48 

381 

73.3 

12 

52 

— 

76.3 

13 

56 

— 

— 

14 

60 

— 

84.4 

15 

64 

— 

87.9 

16 

68 

— 

92.2 

17 

73 

— 

96;5 

18 

77 

— 

100.6 

19 

81 

— 

104.3 

20 

85 

— 

108.9 

21 

89 

— 

114.3 

22 

93 

— 

117.2 

23 

98 

— 

— 

24 

102 

— 

— 

25 

106 

— 

— 

a  Twenty-five  framing  camera  photographs  were  obtainedifor  the  explosive-filled  Armco  iron  cylinder  experiment.  Thetime  be¬ 
tween  frames  is  4.17  ps.  The  zero  for  the  time  scale  was  taken  at  the  time  when  theexpiosive  detonation  front  was  located  at 
the  entrance  end  of  the  cylinder.  The  detonation  frontiwas  assumed  to  be  located  at  the  exit  end  of  the  381-mm-long  cylinder 
in  frame  11  (see  Appendix  K).  Using  this  information  and  the  8.02-km/s  detonation  velocity  for  composition  B  explosivc^O  the 
time  corresponding  to  each  frame  and  to  the  axial  position-of  the  detonation  front  for  frames  1  through  10  was  obtained,  The 
composition  B  explosive  detonation  front  entered  the  end  of  the  cylinder  about  8  ps  after  initiation  of  the  detonation.  The  ex¬ 
plosive  extended  about  64  mm  beyond  the  end  of  thecylinder  on  the  detonation  side.  The  photograph  for  frame  Dismissing. 
bThe  cylinder  outer  radius.was  measured  at  a  position  of  178  mm  (or  47%)  from  the  detonation-front  entrance  end  of  the  cylinder 
(the  right  end  in  Figure-28).  This  position  corresponds  to  the  seventh  vertical  gradicule  line  that  intersects  the  top  of  the  cylinder 
in  Figure  28.  Gradicule:lines  O  and  15  correspond  to  the  right  and  left  ends  of  the  cylinder,  respectively,  in  this  figure.  Two 
methods  were  used  to  measure  the  expanding  cylinder-radius:  The  diameter  was  measured  directly  until  the  lower  surface  of 
tire  cylinder  became  obscured:  in  frame  20.  The  position  of  the  upper  surface  was  measured  from  a  fixed  horizontal  gradicule 
line  until  the  upper  surface  became  obscured  by  smoke  in  frame  23.  The  radius  values  listed  for  frames  6  through  19  (excluding 
frame  13)  are  average  values. for-these  two  methods-The  valucs  listed  for  frames  20,  21,  and  22  were  obtained  with.the  second 
second  method.  The  gradicule  lines  are  25.4  mm  apart. 
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Figure  30  gives  the  Arinco  iron  cylinder  outer  radius  for  various  times  after  explosive  de¬ 
tonation.  The  maximum  observed  increase  in  radius  was  about  150%.  The  angle  0  between  the 
expanding  outer  surface  and  the  initial  position  of  the  cylinder  wall  was  observed  to  increase  with 
time.  This  angle  increased  from  0.07  to  0.14  rad  in  a  continuous  manner  as  the  outer  surface 
moved  from  its  position  at  14  ps  to  its  position  at  106  ps.  The  Tayjor  formula,21  V  =  2D  sin(0/2), 
where  D  is  the  explosive  detonation  velocity,  can  be  used  to  estimate  the  cylinder  velocity  V  dur¬ 
ing  the  final  stages  of  cylinder  expansion.  Using  0  =  0.14  rad  and  D  =  8.02  km/s  for  composition  B 
explosive20  gives  V  =  1.1  km/s.  This  value  is  in  good  agreement  with  the  velocity  value  of  1  km/s 
estimated  from  Figure  29. 


381  318  254  191  127  64  0 

DISTANCE  ALONG  CYLINDER  (mm) 


Figure  30.  Armco- iron  cylinder  outer  surface  at  various  times  after  explosive  detonation.  All 
times  are  measured  with  respect  to  when  the  detonation  front  entered  the  cylinder. 
The  position  of  the  deton ation ;f font  in  the  cylinder  is  also  given.  The  cylinder  surface 
is  represented  by  a  dashed  line  when  cracks  were  observed  at  that  position.  The 
surface  positions  at  the  times  14,  31,  48,  64,  81,  98,  and  106  ps  correspond  to 
framing  camera  photographs  3,  7,  11,  15,  19,  23,  and  25,  respectively,  which  are 
presented  in  Appendix  K. 


Figure  31  shows  the  initiation  and  growth  of  a  series  of  Armco  iron  cylinder  surface  cracks 
as  a  function  of  time.  The  cracks  were  first  observed  about  40  to  50  ju s  after  the  detonation 
front  entered  the  cylinder.  They  may  have  been  initiated  sooner  but  were  not  visible.  About  20 
surface  cracks  were  initiated  in  the  region  of  observation.Thc  cracks  were  about  10  to  20  mm  in 
length.  After  about  60  jus  of  growth  and  coalescence,  the  number  of  surface  cracks  decreased  to 
about  half  the  original  number,  but  the  average  crack  length  increased  to  about  30  to  50  mm.  This 
type  of  crack  data  may  be  more  appropriate  for  obtaining  dynamic  fracture  parameters  for  use  in 
exploding  cylinder  simulations  than  the  one-dimensional  spall  fracture  experiments  described  in 
the  previous  sections.  The  dynamic  fracture  parameters  may  be  very  different  for  the  two  types  of 
experiments.  Fewer  cracks  would  probably  be  produced  for  framing  camera  experiments  (tens  of 
cracks  over  a  large  area  in  Figure  31)  as  compared  to  spall  fracture  experiments  (in  some  cases, 
hundreds  of  cracks  in  a  small  area,  see  Appendix  12)  due  to  the  different  strain  rates.  It  may  also 
be  possible  to  use  framing  camera  cylinder  experiments  to  observe  the  growth  and  coalescence  of 
surface  shear  cracks  for  those  materials  in  which  shear  fractures  occur  through  the  entire  v/all 
of  theicylindcr. 


VIII.  EXPERIMENTAL  DETERMINATION  OF  FRAGMENT  MASS  DISTRIBUTIONS 

A  series  of  Armco  iron  and  I1F-1  steel  exploding  cylinder  experiments  and  an  HF-1  steel 
projectile  experiment22  were  performed  to  provide  fragment  mass  distributions  for  comparison 
with  the  computations  and  for  future  model  development  of  the  fracture  processes.  Theseex- 
periments  were  conducted  at  the  NSWC  Sawdust  Pit  Facility.23  Table  13  gives  the  configuration 
details  for  the  experiments.  A  schematic  of  an  explosive-filled  cylinder  experiment  is  shown 
in  Figure  32.  The  cast-in-place  composition  B  explosive  extendedibeyond  the  cylinder  to  minimize 
end  effects.  Figure  33  shows  the  explosive-filled  HF-1  steel  cylinder  for  experiment  5.  A  pro¬ 
jectile  similar  to  the  one  used  in  experiment  7  is  shown  in  Figure  34.  The  projectile  was  sectioned 
to  show  its  inner  surface  and  wall  thickness. 


DISTANCE  ALONG  CYUNOEA  (mm) 

Figure  31.  Evolution  of  Armco  iron  cylinder  surface  cracks  as  a  function  of  time.  The  times 
after  the  detonation  wave  has  enteredthe  cylinder  are  given  in  the  upper  right  corner 
of  each  frame.  The  times  at  127  and  247  mm  on  each  frame  correspond  to  the  times 
after  the  detonation  front  has  reached  these  positions,  respectively.  The  calculated 
position  of  the  detonation  front  is  shown  in  the  first  three  frames.  The  position  of  the 
cylinder  outer  surface  is  shown  in  each  frame.  A  dashed  line  represents  the  surface 
after  cracks  were  observed  at  that  position.  The  inner  surface  of  the  cylinder  is  also 
shown  prior  to  shock  compression.  The  sloped  dashed  line  in  the  cylinder  wall  in  the 
first  three  frames  separates  shocked  and  unshocked  material.  The  evolution  of  a  group 
of  cracks  in  the  center  of  the  cylinder  is  shown  as  they  become  visible  and  grow. 
TheseTrames  (beginning  at  the  bottom)  correspond  to  framing  camera  photographs  7, 
9,  11,  15,  17,  19,  21,  23,  and  25,  respectively.  The  scale  for  the  outer  radius  of  the 
cylinder  is  labeled  only  in  the  first  frame. 
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Table  13.  Configuration  detaiis  for  explosive-filled  cylinder  and  projectile  experiments. 


Experiment 

No.a 

Configuration 
and  Material 

Length 

(mm) 

Inside 

Diameter 

(mm) 

Outside 

Diameter 

(mm) 

Wall 

Thickness  Mass 
(mm)  (kg) 

Hardness 

Explosive 

Mass’* 

(kg) 

Fraction  of 
Explosive  Mass 
Inside  Cylinder® 

Explosive  Mass 
to  Metal 
Mass  Ratio** 
(C/M) 

1 

Armco  Iron  Cylinder 

203.0 

76.2 

1 14.3 

19.1 

9.13 

RB40 

2.75 

0.57 

0.17 

2 

Armco  Iron  Cylinder 

203.0 

76.2 

114.3 

19.1 

9.16 

rb40 

2.77 

0.57 

0.17 

3 

HIM- Steel  Cylinder, 
Heat  Treatment  A 

203.0 

76.2 

120.7 

22.3 

10.94 

Rc39 

2.75 

0.57 

0.14 

4 

HIM  Steel  Cylinder, 
HeatTieatmcm  A 

‘'O.O 

76.2 

120.7 

22.3 

10.96 

Rc39 

2.75 

0.57 

0.14 

5 

HIM  Steel  Cylinder, 
Heat  Treatment  B 

202.3 

75.9 

120.0 

22.1 

10.57 

Rc36e 

2.76 

0.57 

0.15 

6 

HIM  Steel  Cylinder, 
HeauTrcatmcnt  B, 

202.4 

75.9 

120.4 

22.3 

10.77 

Rc36c 

2.68 

0.57 

0.14 

7f 

HF-1  Steel  Projectile 

565.4s 

90.  lh 

126. 1-1 

18.011 

25.21* 

Rc38-39 

3.11 

Not  Applicable 

0.12* 

a  Explosive  experiments  1  through  7  correspond,  respectively,  to  NSWC/DL  Fragmentation  Tests  2872, 2873,  2874, 2875, 2924,  2925,  and  2824. 

"  Cast-in-place  composition  B  explosive  was  used  in  experiments  1  through  6.  PBXN-106  explosive  cast  in  a  0.20-kg  mass,  1.52-mm-thick 
polyethylene  encapsulating  beaker  was  used  in  experiment  7. 

e  To  minimize  end  effects  the  explosive  extended  approximately  90  mm  beyond  the  end  of  the  cylinder  on  the  detonator  end  and  approximately 
64  mm  beyond  the  end  of  the  cylinder  on  the  opposite  cud.  The  diameter  of  the  explosive  was  the  same  both  inside  and  outside  the  cylinder. 
The  inside  surface  of  each  cylinder  was  painted  with  MIL-C-450  cavity  paint  prior  to  inserting  the  composition  B  explosive.  The  explosive  was 
detonated  with  an  Engineer's  Special  blasting  cap  that  was  placed  in  a  7.14-mm-diameter  by  6.35-mrn-decp  hole  in  the  end  of  the  explosive. 

“  For  tlie  cylinder  experiments,  calculated  using  the  estimated  explosive  mass  inside  the  cylinder. 

c  Average  of  nine  measurements  on  each  end  of  the  cylinder.  The  hardness  measurements  varied  from  about  Rq  35  to  37  from  the  inside  to  the 
outside  cylinder  wall. 

The  projectile  wasnosc-initiatcd  with  an  M564  fuse  that  liadEbccn  modified  for  still  detonation.  The.  fuse  mass  was  0.84  kg.  In  addition  to  the 
hardness  values  listed  in  this  table,  other  mechanical  properties  of  the  projectile  material  include  the  following:  1.06-GPa  0.2%  yield  strength, 
1.2-GPa  ultimate  tensile  strength,  8-9%  elongation,  and  14-1 8%  reduction  in  area.  The  projectile  dimensions  listed  in  this  table  were  not  measured 
for  the  projectileiuscd  in  experiment  7  but  for  a  projectile  from  the  same  lot.  The  diameter  and  thickness  tolerances  between  projectiles  arc  a 
few  percent. 

®  Length  of  the  assembled  two-piece  projectile,  but  not  including  the  length  of  the  fuse. 

.  These  values  correspond  to  the  102-mm-Iong  uniform  region  of  the  aft  projectile  section  beginning  at  the  joint. 

!  Sum  of  masses  of  the  forward  (1 1.59  kg)  and  aft  (13.62  kg)  projectile  sections. 

1  The  C/M  ratio  for  the  102-mm-long  uniform  region  of  the  aft  projectile  section  is  0.21;  this  value  was  obtained  by  using  7.77  and  1.64  Mg/m  , 
respectively,  for  the  densities  of  HIM  steel  and  PBXN-106  explosive. 
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METAL  CYLINDER 


After  detonation  the  fragments  for  each  experiment  were  captured  in  a  large  volume  of 
sawdust  and  then  recovered  by  mechanical  screen  and  magnetic  separation  techniques.  Mass 
distributions  were  determined  by  weighing  and  counting  the  fragments.  Fragments  smaller,  than 
1  gr  (0.0648  g)  were  not  counted  due  to  their  large  numbers.  At  least  98.9%  of  the  metal  mass 
was  recovered  in  all  the  experiments.  For  these  thick-wajjedicylinder  and  projectile  experiments 
both  brittle  and  shear  fracture-surfaces  were  observed. 

The  observation  of  brittle  and  shear  fractures  in  cylinder  fragmentation  was  reported  in  the 
eaiiy  work  of  Mott.24  For  the  purpose  of  the  present  work  a  fragment  classification  scheme  was 
devised  to  further  characterize  the  recovered  fragments  according  to  their  cylindrical  andifracture 
surfaces.  The  brittle  fracture  surfaces  usually  occurred  near  the  cylinder  (or  projectile)  outer 
surface;  the  shear  fracture  surfaces  usually  occurred  near  the  cylinder  (or  projectile)  inner  surface. 
Table  14  defines  the  eight  fragment  types.  Schematic  examples  of  these  fragment  types  are-given 
in  Figure  35.  Photographs  of  the  fragment  types  from  the  Armco  iron  and  HF-1  steel  cylinder 
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Figure  33.  Explosive-filled  HF-1  steel  cylinder  with  heat  treatment  B  for  experiment  5.  (a)  After 
fabrication  and  (b)  in  styrofoam  box  with  lid  removed  prior  to  filling  of  sawdust 
recovery  pit. 


Figure  34.  Photographs  of  :{a)  outer  surface  and  (b)  inner  surface  of  sectioned  HF-1  steel  pro¬ 
jectile  with  heat  treatment  A. 

experiments  are  given  in  Appendix  L.  Figure  36  is  a  schematic  of  a  cross  section  through  the 
cylinder  or  projectile  wa|!  showing  the  possible  formation  of  the  fragment  types  during  fragmen¬ 
tation.  It  is  not  known  whether  all  the  fragment  types  were= formed  at  about  the  same  time  after 
explosive  detonation  or  whether  they  were  formed  at  different  times.  It  may  be  possible  that 
many  type  1  and  type  2 A  fragments  formed  initially  and-thenlater  broke  up  into  fragmenUtypes 
2B,  3A,  4A,  and  4B.  A  fragment  characterization  scheme  that  is  similar  to  the  present  one,  but 
limited  to  four  fragment  types,  has  been  independently  devised.25 

Appendix  M  contains  the  detailed  results  for  the  experiments.  Tables  M-l  through  M-7  give 
the  mass  distribution  results  for  the  typed  fragments  for  each  experiment.  Tables  M-8  through 
M-l 4  give  dimension  measurements  and  strain  calculationsTor  selected  type  1  fragments  from  each 


Table  14.  Definitions  of  fragment  types  used  for  characterization  of  recovered  fragments 
from  the  explosive-filled  cylinder  and  projectile  experiments. 


Fragment 

Type 

Definition* 

1 

Has  both  the  inner  and  outer  surfaces  of  the  cylinder 
(or  projectile).  Usually  has  brittle  fracture  surfaces 
near  the  outer  surface  and  shear  fracture  surfaces  near 
the  inner  surface. 

2A 

Has  the  outer  surface  of  the  cylinder  (or  projectile), 
and  both  brittle  fracture  (usually  near  the  outer  surface) 
and  shear  fracture  surfaces. 

2B 

Has  the  outer  surface  of  the  cylinder  (or  projectile), 
and  only  brittle  fracture  surfaces. 

3A 

Has  the  inner  surface  of  the  cylinder  (or  projectile), 
and  both  brittle  fracture  and  shear  fracture  (usually 
near  the  inner  surface)  surfaces. 

3B 

Has  the  inner  surface  of  the' cylinder  (or  projectile) 
and  only  shear  fracture  surfaces. 

4A 

Has  only  brittle  and  shear  fracture  surfaces. 

4B 

Has  only'brittle  fracture  surfaces. 

4C 

Has  only  shear  fracture  surfaces. 

a  These  definitions  do  not  include  the  fragment  end  surfaces  (the  end  surfaces  were  usually  perpendicular  to  the  cylinder  axis).  The 
fragment  end  surfaces  usually  were  brittle  fracture  surfaces.  Sometimes  a  fragment  end  surface  contained  the  end  surface  of  the 
cylinder  or  a  shear  fracture  surface. 
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Figure  35  Schematic  examples  of  fragment  types  from  the  (a)  Armco  iron  and  (b)  HF-1  steel 
cylinder  experiments.  These  examples  may  not  correspond  to  the  most  common 
fragments  of  a  given  type.  The  inner,  outer,  and  end  surfaces  of  the  cylinder  are 
indicated  by  the  letters  I,  0,  and  E,  respectively.  Brittle  and  shear  fracture  surfaces 
are  indicated  by  the  letters  B  and  S,  respectively.  Two  views  are  shown  for  each 
fragment;  a  side  view  (looking  perpendicular  to  the  cylinder  axis)  and  an  end  view 
(looking;  parallel  to  the  cylinder  axis).  The  shear  fracture  surfaces  made  angles  of 
approximately  45°  with  the  inner  cylinder  surface  (fragment  end  view).  No  type 
4C  fragments  were  observed  for  the  Armco  iron  experiments. 


CYLINDER  OR  PROJECTILE 
OUTER  SURFACE 


INNER  SURFACE 

Figure  36.  Schematic  through  cylinder  or  projectile  wall  showing  the  possible  formation  of 
fragment  types  from  the  intersections  of  brittle  and  shear  fracture  surfaces.  The 
fragment  types  are  defined  in  Table  14. 

Figure  37  is  a  semilogarithmic  plot  of  the  number  of  fragments  with  mass  greater  than  m 
versus  the  square  root  of  the  fragment  mass  for  each  experiment.  Mott  and  Einfoot26  found 
experimentally  that  the  fragment  mass  distributions  for  many  projectiles  and  cylindrical  warheads 
could  be  approximated  by  an  empirical  equation  of  the  form 

N(>  in)  =  N0  exp  [-(m/n)v‘]  (4) 

where  N  (>m)  is  the  number  of  fragments  with  mass  greater  than  m,  NQ  is  the  total  number  of 
fragments,  m  is  the  fragmentimass,  and  n  is  a  fragmentation  size  parameter.  Exponents  other  than 
1/2  have  been  used  in  Equation  (4)  for  representing  fragment  mass  distribution  data.6  Figure  37 
suggests  that  Equation  (4)  may  not  adequately  represent  the  fragment  mass  distribution  data  for 
the  present  experiments,  except  possibly  for  the  Armco  iron  experiments,  if  the  small  fragment 
masses  are  excluded.  Fragment  mass  distribution  results  are  plotted  versus  the  fragment  mass  in 
Figure  38.  Figure  38(a)  is  a  plot  of  the  number  of  fragments  with  mass  greater  than  m  and 
Figure  38(b)  is  a  plot  of  the:=average  number  of  fragments  per  1-gr  interval.  To  construct  Figure 
38(b)  some  of  the  mass  groups  for  an  experiment  (given  in  Tables  M-l  through  M-7)  were 
combined  into  larger  mass  groups  before  calculating  the  average  number  of  fragments,  especially  if 
the  smaller  mass  groups  hadtonly  a  few  fragments.  The  values  in  Figure  38(b)  were  plotted  at  a 
fragment  mass  correspondingito  the  midpoint  of  the  appropriate  mass  group.  Figure  38(a)  suggests 
that  for  the  HF-1  steel  cylinder  experiments  an  exponent  of  1  rather  than  1/2  in  Equation  (4)  is  a 
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more  adequate  representation  of  the  fragment  mass  distribution  data.  Table  15  summarizes  the 
average  fragment  mass  results  for  the  cylinder  and  projectile  experiments.  The  average  fragment 
masses  arc  97  gr  for  the  Arnico  iron  cylinder  experiments,  15  gr  for  the  HF-1  steel,  heat  treatment 
A,  cylinder  experiments,  12  gr  for  the  HF-1  steel,  heat  treatment  B,  cylinder  experiments,  and  21 
gr  for  the  HF-1  steel,  heat  treatment  A,  projectile  experiment.  The  HF-1  steel  cylinder  results 
indicate  that' heat  treatment  B  material  is  more  brittle  than  heat  treatment  A  material. 


Figure  37.  Fragment  mass  distributions  (versus  the  square  root  of  the  fragment  mass)  for  the 
cylinder  and  projectileexperiments. 
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Figure  38.  Fragment  mass  distributions  (versus  the  fragment  :mass)  for  the  cylinder  and  pro¬ 
jectile  experiments.  (a)  Cumulative  number  of  fragments,  (b)  Average  number  of 
fragments. 
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It  can  be  shown  that  the  fragmentation  size  parameter  m  in  Equation  (4)  is  related  to  the 
average  fragment  mass  m,  where 


J-00 

I  mdN 
_  0 

m  = -  .  (5) 

/OO 

dN 


Here  dN  is  the  number  of  fragments  between  m  and  m  +  dm.  From  Equation  (4)  dN  =  N0e"xdx 
wiiere  x  has  been  substituted  for  (m .  Equation  (5)  then  becomes 


m 


H  J e"x  x2  dx 


dx 


(6) 


Evaluating  the  integrals  gives  n  =  m/2.  It  is  possible  to  derive  a  general  relationship  between 
H  and  Th  to  account  for  the  different  exponents  that  can  be  used  in  Equation  (4).  Rewriting 
Equation  (4)  for  a  general  exponent  gives 


N  ( >  m)  =  N0  exp  [  -  (m//*)1  /n  ],  (7) 

where  n  =  1,2, 3, ... .  Using  x  =  (m/|i)1/n ,  Equation  (5)  becomes 

,/e-x  xn  dx 

m  =  — -  .  (8) 

/V  dx 
0 

Evaluating  Equation  (8)  gives  n  =  m/n!  .  Therefore,  for  exponents  of  1  and  1/3  in  Equation  (7) 
M  =  i h  and  rii/6,  respectively. 
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Table  15.  Average  fragment  mats  results  for  cylinder  and  projectile  experiments. 


Experiment 

No. 

Total  Number 
of  Fragments 
Greater  than  1  gr* 

Fragment  Mass  Greater 
than  1  gr  * 

(kg) 

Average  Fragment 
Mass  for  Fragments 
Greater  than  1  gr 
(gr)  (g) 

Total  Recovered 
Fragment  Mass*1 
(kg) 

Fraction  of  Metal 
Mass  Recovered 

1 

1395 

8.99 

99.5 

6.45 

9.08 

0.995 

2 

1457 

9.01 

95.4 

6.18 

9.12 

0.996 

3 

10181 

10.33 

15.7 

1.01 

10.89 

0.996 

4 

10740 

10.32 

14.8 

0.96 

10.93 

0.997 

5 

12132 

9.70 

12.3 

0.80 

10.46 

0.989 

6 

12927 

9.88 

11.8 

0.76 

10.65 

0.989 

f 

16403 

22.23 

20.9 

1.36 

25.07 

0.994 

*  Fragments  with  a  mass  less  than  I  gr  (0.0648  g)  were  not  counted. 

“Approximately  1,  5,  and  7%  of  the  total  recovered  mass  was  in  theO  to  1-gr  mass  group  for  cylinder  experiments  1  and  2,  3  and 
4,  and  5  and  6,  respectively.  Approximately  3%  of  the  total  recovered  fragment  mass  was  in  the  0  to  1-gr  mass  group  for  the 
projectile  experiment  7. 

c  Includes  only  fragments  from  the  forward  and  aft  projectile  sections  and  not  the  fuse  fragments. 


Figures  39  through  42  give  the  fragment  mass  distribution  results  for  the  typed  fragments 
for  the  cylinder  and  projectile  experiments.  Fragment  types  2,  3,  and  4  are  the  sums  of  types 
2A  and  2B,  types  3A  and  3B,  and  types  4A,  4B,  and  4C,  respectively.  Table  16  summarizes  the 
results  for  the  typed  fragments  for  the  experiments.  The  average  fragment  mass  results  for  the 
typed  fragments  are  given  in  Table  17. 


For  the  Armcoiron  cylinder  experiments  only  type  1  fragments  were  observed  for  the  larger 
fragment  masses.  Fopexperiment  2,  fragment  types  2,  3,  and  4  were  observed  for  fragment  masses 
less  than  about  400,  200,  and  25  gr,  respectively.  No  type  4C  fragments  were  observed  for  this 
experiment.  The  number  of  type  1  fraginents  increased  with  decreasing  fragment  mass  except  for 
the  smallest  masses.  About  25%  of  therecovered  fragments  were^type  1.  For  the  HF-1  steel,  heat 
treatment  A,  cylinder  experiment  4,  fragment  types  1,  2,  and:  3  were  observed  for  the  largest 
fragment  masses.  Type  4  fragments  were  observed  for  masses  less  than  about  125  gr.  ^Fragment 
types  1  and  2  were  observed  for  the  largcstTragment  masses  for  the  heat  treatment  B,  experiment 
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5.  Fragment  types  3  and  4  were  observed  for  fragment  masses  less  than  about  175  and  105  gr, 
respectively.  Only  about  1%  of  the  fragments  from  the  HF-1  steel  cylinder  experiments  were  type  1. 
About  35 %  of  the  fragments  were  type  4  as  compared  with  about  4%  type  4  fragments  from  the 
Armco  iron  experiment  2. 

For  the  projectile  experiment  7  only  type  1  fragments  were  observed  for  the  larger  fragment 
masses.  Fragment  types  2,  3,  and  4  were  observed  for  fragment  masses  less  than  about  550, 450, 
and  150  gr,  respectively.  Almost  all  of  the  large  type  1  fragments  originated  from  the  nose  region 
of  the  projectile.  The  reconstructed  nose  is  shown  in  Figure  43.  The  reconstructed  nose  mass  was 
about  2.5  kg,  which  corresponds  to  about  10%  of  the  initial  projectile  mass.  Since  the  nose  region 
of  the  projectile  was  not  in  direct  contact  with  the  explosive  charge,  this  region  was  probably 
subjected  to  a  lower  strain  rale  than  the  remainder  of  the  projectile.  Figure  44  shows  that  the 
projectile  fragment  mass  distribution  is  changed  appreciably  if  the  large  nose  fragments  are 
excluded.  The  mass  of  the  largestTragments  decreases  from  about  1600  to  800  gr.  The  mass  of  the 
largest  fragments  for  the  cylinder  experiments  was  about  400  gr. 
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Figure  39.  Fragment  mass  distributions  for  the  Armed  iron  cylinder  experiments. 

(a)  Cumulative  number  of  fragments,  (b)  Average  number  of  fragments. 
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OF  FRA6MINTS  WITH  MASS  CREATER  THAN  M 
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Figure  40.  Fragment  mass  distributions  for  the  HF-1  steel,  heat  treatment  A,  cylinder  experi¬ 
ments.  (accumulative  number  of  fragments,  (b)  Average  number  of  fragments. 
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Table  16.  Summary  of  results  for  typed  fragments  from  cylinder  and  projectile  experiments. 


Average  fragment  mass  results  for  typed  fragments  from  cylinder  and  projectiles  experiments 


The  dimension  measurements  given  in  Tables  M-8  through  M-14  for  selected  type  1 
fragments  from  the  cylinder  and  projectile  experiments  include  the  fragment  mass,  wall  thick¬ 
ness,  maximum  fragment  width  and  length,  and  average  fraction  of  wall  thickness  with  brittle 
fracture.  Radial  and  plastic  strain  values  were  calculated  using  the  initial  and  final  wall  thickness 
values.  The  measurements  and  calculations  are  summarized  in  Table  18.  The  values  for  the 
average  fraction  of  the  wall  thickness  with  brittle  fracture  are  0.56, 0.76,  and  0.80  for  the  Armco 
iron,  HF-1  steel  with  heat  treatment  A,  and  HF-1  steel  with  heat  treatment  B  cylinder  experi¬ 
ments,  respectively.  This  value  is  0.69  for  fragments  from  the  central  uniform  region  of  the 
HF-1  steel  projectile.  These  results  agree  with  the  decreased  ductility  of  these  materials  as  given 
by  the  percent  elongation  values  in  Table  3.  Average  engineering  radial  strain  values  for  the  type  1 
fragments  from  the  Armco  iron,  HF-1  steel  with  heat  treatment  A,  and  HF-1  steel  with  heat 
treatment  B  cylinder  experiments  are  0.35,  0.12,  and  0.13,  respectively.  This  value  is  0.11 
for  fragments  from  the  central  uniform  region  of  the  HF-1  steel  projectile.  It  is  possible  to 
estimate  the  final  outside  diameters  corresponding  to  these  residual  strain  values  using  the 
formula  (dQ  -  aQ)  aQ  =  (df  -  af)  af.  Here  dQ  and  df  are  the  initial  and  final  outside  diameters, 
respectively,  and  ao  and  af  are  the  initial  and  final  wall  thicknesses,  respectively.  This  formula 
was  obtained  by  assuming  that  the  cross-sectional  area  of  the  cylinder  or  projectile  remained 
constant  during  explosive  expansion.  Using  diameter  and  thickness  values  from  Tables  13  and  18, 
the  average  percent  increase  in  the  outside  diameters  corresponding  to  the  residual  strain  values  is 
39,  8,  and  10%  for  the  Armco  iron,  HF-1  steel  with  heat  treatment  A,  and  HF-1  steel  with  heat 
treatment  B  cylinders,  respectively.  The  percent  increase  for  the  HF-1  steel  projectile  is  9%. 
These  results  indicate  that  during  explosive  expansion  the  Armco  iron  cylinders  were  at  i  much 
larger  diameter  than  the  HF-1  steel  cylinders  and  projectile  when  the  increase  in  the  plastic  radial 
strain  effectively  ceased. 


Table  18.  Summary  of  measurements  on  type  1  fragments  from  cylinder  and  projectile  experiments.* * 


Average  Maximum 


Experiment 

No. 

Total  Number 
of  Type  1 
Fragments 

Number  of 

Type  1  Fragments 
Measured 

Fragment 

Mass 

(|»)  (I) 

Wall 

Thickness 

(mm) 

Fraction  of  Wall 
•  Thickness  with 
Brittle  Fracture1* 

Fragment 

Width' 

(mm) 

Fragment 

Length1* 

(mm) 

Engineering 

Radial 

Stralne 

True 

Radial 

Homogeneous 

Plastic 

Strain* 

1 

206h 

25 

550.11 

35.64* 

12.6 

05  j 

9 

66 

0.34 

0.42 

0.72 

17 

2 

195 

25 

575.71 

37.3 11 

12.3 

0.57 

10 

72 

0.36 

0.45 

0  78 

17 

3 

33h 

33 

191.9 

12.43 

19,7 

0.77 

8.3 

21.8 

0.12 

0.12 

0.22 

9.9 

4 

25 

25 

175.9 

11.40 

19.9 

0,75 

6.9 

21.0 

0.11 

0,11 

0.19 

8.7 

5 

55 

55 

146.1 

9.46 

19.0 

0.80 

7.3 

19.4 

0.14 

0.15 

0.26 

9.2 

6 

49 

49 

161.3 

10.48 

19.5 

0.80 

7.9 

21.4 

0.12 

0.13 

0.22 

10.0 

7 

472 

10* 

203.71 

13.20* 

16.0 

0  69 

11.0 

40.9 

0.11 

C.12 

0.20 

* 

8.4 

5k 

8.7 

0.90 

6.6 

19.3 

0.10 

0.10 

0.17 

6.5 

5« 

7.3 

0.78 

5.5 

19.1 

0.16 

0.18 

0.30 

5.6 


?  The  valuta  lilted  in  this  toble  arc  average  values  for  the  fragment  measurements  given  in  Tables  M*8  through  M*|4, 

®  These  values  are  the  average  of  measurements  from  the  cylinder  (or  projectile)  outer  surface  to  the  beginning  of  the  shear  fracture  surface  on  each  side  of  the  fragment. 
The  estimated  uncertainty  in  these  values  Is  10-20%. 

cThe  two  values  listed  arc  the  average  maximum  width  of  the  brittle  fracture  region  near  the  cylinder  (or  projectile)  outer  surface  and  the  average  maximum  width  of  the 
shear  fracture  region  near  the  cylinder  (or  projectile)  Inner  surface. 

“The  fragment  length  is  parallel  to  the  cylinder  (or  projectile)  axis, 

e  The  engineering  radial  strain  is  (a0-af)/a0  where  a0  and  af  are  the  initial  and  final  wall  thicknesses,  respectively, 

*  The  true  radial  strain^  Is  -In  (*{/i0). 

®  The  homogeneous  plastic  strain*  is  -  v3  ln(af/a0), 

j'Six  fragments  were  removed  from  the  distribution  prior  to  the  characterisation  of  the  fragments;  therefore  the  total  number  of  type  1  fragments  may  be  greater, 
j  Average  value  for  total  number  of  type  I  fragments. 

*  fragments  from  the  102-mm-Iong  uniform  region  of  the  aft  projectile  section. 

*  Fragments  from  the  joint  region  of  the  forward  projectile  section. 

'  Fragments  from  the  joint  region  of  the  aft  projectile  section. 
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IX.  COMPUTATION  OF  FRAGMENT  MASS  DISTRIBUTIONS 


In  this  section,  computed  fragment  mass  distribution  results  are  presented  for  Armco  iron 
and  HF-1  steel  explosive-filled  cylinders  and  an  HF-1  steel  projectile.  The  computations  are 
compared  with  the  experiments  of  the  last  section.  Cylindrical  PUFF  and  TROTT  computations 
were  performed  to  simulate  the  exploding  cylinders  and  the  projectile,  respectively.  The  fracture 
processes  were  simulated  with  the  BFRACT2  brittle  fracture  model  and,  the  SHEAR2  shear  band 
model.  Brittle  and  shear  fracture  crack  distribution  results  were  obtained  from  the  computations. 
Fragment  mass  distributions  were  then  calculated  from  the  crack  distributions. 

Computer  listings  of  the  input  parameters  and  the  output  results  for  the  PUFF  and  TROTT 
computations  are  given  in  Apperidix  N.  Four  computations  (Figures  N-l,  N-2,  N-3,  and  N-6) 
were  performed  using  the  BFRACT2  model1  to  simulate  the  fracture  processes  (the  recovered 
fragments  had  more  brittle  than  shear  fracture  surface  area,  see  Table  18).  These  computations 
were  for  explosive-filled  cylinders  of  Armco  iron  and  HF-1  steel  with  heat  treatments  A  and  B,  and 
for  an  HF-1  steel  projectile  with  heat  treatment  A.  Two  computations  (Figures  N-4  and  N-5)  were 
performed  using  the  BFRACT2  and  SHEAR2  fracture  models  together  (the  recovered  fragments 
had  both  brittle  and  shear  fracture  surfaces).  These  computations  were  for  explosive-filled 
cylinders  of  HF-1  steel  with  heat  treatments  A  and  B. 

Input  parameters  for  the  composition  B  explosive,  for  Armco  iron  (excluding  the  yield 
strength,  nucleation  threshold  stress,  and  strain  hardening  parameter),  and  for  the  HF-1  steel 
SHEAR2  subroutine  were  obtained  from  Reference  1.  The  PBXN-106  and  CH6  explosive  input 
parameters  and  the  cell  and  coordinate  layout  for  the  projectile  computation  were  obtained  from 
Reference  3.  The  HF-1  steel  input  parameters  for  the  Hugoniot  equation  of  state  and  for  the 
BFRACT2  subroutine  (excluding  the  nucleation  threshold  stress)  were  obtained  from  Sections  IV 
and  VI,  respectively.  The  longitudinal-direction  yield  strength  values  of  0.18  and  1.04  GPa  for 
Armco  iron  and  HF-1  steel  with  heat  treatment  A,  respectively,  obtainedTrom  Table  3,  were  used 
in  the  computations.  The  transverse-direction  yield  strength  values  in  Table  3  were  not  used,  since 
only  isotropic  material  properties  can  be  treated  in  the  computations.  The  longitudinal  tensile- 
pull  direction  is  parallel  to  the  cylinder  axis  for  the  exploding  cylinder  experiments.  Since  a 
longitudinal-direction  yield  strength  value  was  not  determined  for  HF-1  steel  with  heat  treatment 
B,  the  heat  treatment  A  value  was  used  for  all  the  HF-1  steel  computations.  The  larger  1.36-GPa 


yield  strength  value  for  HF-1  steel  determined  from  the  Hugoniot  measurements  (Section  IV)  was 
used  only  in  the  PUFF  computations  for  the  high-strain-rate  spall  fracture  experiments 
(Appendix  G).  Armco  iron  and  HF-1  steel  nucleation  threshold  stress  values  of -0.18  and  -1.04 
GPa,  respectively,  were  used  in  the  computations.  For  these  simulations,  it  was  assumed  that 
cracks  nucleated  when  the  tensile  stress  exceeded  the  yield  strength. 

Strain  hardening  parameters  for  Armco  iron  and  HF-1  steel  with  heat  treatment  A  were 
obtained  from  the  true  stress-true  strain  curves  in  Figures  4  and  5,  respectively.  The  strain  harden¬ 
ing  parameter  is  defined  as  the  increase  in  the  yield  strength  from  /z  =  /liel  to  n  =  0.2,  where 
the  true  strain  n  =  V0/V  -  1  and  the  elastic  limit  strain  =  Y/(2G).*7  Here  VQ  and  V  are  the 
initial  and  final  specific  volumes,  respectively,  Y  is  the  initial  yield  strength,  and  G  is  the  shear 
modulus.  An  Armco  iron  strain  hardening  parameter  of  0.2}  GPa  was  obtained  from  Figure  4  from 
the  difference  in  the  true  stress  value  of  0.39  GPa  at  g  =  0.2  and  the  elastic  stress  value  of  0.18 
GPa  at  pEL  =  0.001.  An  HF-1  steel  strain  hardening  parameter  of  0.94  GPa  was  obtained  from 
Figure  5  from  the  difference  in  the  extrapolated  true  stress  value  of  1.98  GPa  at  n  =  0.2  and  the 
elastic  stress  value  of  1.04  GPa  at  /iEL  =  0.007.  (The  extrapolated  stress  value  was  obtained  by 
fitting  a  straight  line  through  the  elastic  stress  and  strain  values  and  the  1.41-GPa  stress  and  0.0804 
strain  values  at  fracture.) 

The  crack  distribution  output  results  for  the  computations  are  given  in  Figures  N-l(b), 
N-2(b),  N-3(b),  N-4(b),  N-5(b),  N-6(b),  and  N-6(c).  Results  are  given  for  each  computational  cell 
containing  Armco  iron  or  HF-1  steel.  These  results  are  summarized  and  average  crack  distribution 
values  for  the  metal  cylinder  or  projectile  are  given  in  Tables  O-l,  0-3,  0-5,  0-7,  0-9,  and  0-11. 
For  each  simulation  the  average  crack  density  (NQ)  and  the  average  crack  radius  (R, )  are  obtained 
for  the  activated  cracks  with  each  orientation  by  averaging  the  crack  distribution  values  for  the 
computational  cells.  These  average  values  are  obtained  fromthe  equations 

2  Nn.  V, 


-  SNcV.R, 
R,  =  -i - - — - 1 


2  N0.  V, 


where  N0,  and  R] ,  are  the  crack  density  and  average  crack  radius,  respectively,  for  each  computa¬ 
tional  cell  and  V.  is  the  cell  volume.  The  summation  is  over  the  appropriate  computational  cells. 
The  NQ  and  Rj  values  obtained  for  each  computation  are  given  as  footnotes  in  the  tables  in 
Appendix  0.  Figure  45  gives  plots  of  the  brittle  axial  and  brittle  circumferential  crack 
distributions  for  the  Armco  iron  and  HF-1  steel  cylinder  computations,  and  the  brittle 
circumferential  crack  distribution  for  the  HF-1  steel  projectile  computation  using  the  BFRACT2 
subroutine.  Figure  46  gives  plots  of  the  brittle  axial,  brittle  circumferential,  and  shear  crack 
distributions  for  the  HF-1  steel  cylinder  computations  using  the  BFRACT2  and  SHEAR2 
subroutines.  For  each  computation,  more  axial  cracks  were  activated  than  circumferential  or  shear 
cracks.  Figure  46  indicates  that  for  the  HF-1  steel  cylinder  computations  the  number  of  shear 
cracks  is  a  few  orders  of  magnitude  less  than  the  number  of  axial  or  circumferential  cracks,  but  the 
average  crack  radius  is  larger. 


Figure  45.  Brittle  fracture  crack  distributions  for  the  Armco  iron  and  HF-Tsteellcylinder  PUFF 
computations  and  for  the  HF-1  steel  projectile  TROTT  computation  using  the 
BFRACT2  subroutine.  These  average  crack  distributions  were  determined  from  the 
results  in  Tables  0-1,  0-3,  0-5,  and  0-11. 

71 


Figure  46.  Brittle  and  shear  fracture  crack  distributions  for  the  HF-1  steel  cylinder  PUFF  com¬ 
putations  using  the  BFRACT2  and  SHEAR2  subroutines.  These  average  crack  dis¬ 
tributions  were  determined  from  the  results  in  Tables  0-7,  and  0-9, 


In  the  SHEAR2  subroutine,  shear  bands  nucleate  when  the  plastic  shear  strain  exceeds  a 
threshold  value.  For  these  HF-1  steel  computations,  a  threshold  value  of  0.17  was  used  since 
“contained”  exploding  HF-1  steel  cylinder  experiments  performed  at  SRI  indicate  that  shear 
bands  begin  to  nucleate  at  an  equivalent  plastic  strain  of  0.29  (and  therefore  a  shear  strain  of 
0.29/\/y=  0.17).3  For  the  HF-1  steel  cylinder  computations  using  the  BFRACT2  and  SHEAR2 
subroutines,  brittle  fracture  began  at  the  outer  cylinder  surface  and  propagated  inward  and  shear 
bands  began  at -the  inner  cylinder  surface  and  propagated  outward.  Brittle-fracture  was  nucleated 
in  40%  of  the  HF-1  steel  computational  cells;  shear  bands  were  nucJeatedUn  the  remaining  cells. 
For  the  recovered  HF-1  steel  type  d  fragments,  brittle  fracture  extendedifrom  the  outer  surface 
inward  for  about  75  to  80%  of  the  wall  thickness;  the  remaining  20  to  25%  of  the  wall  thickness 
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was  shear  fracture  (Table  18,  Figures  L-9  and  L-17).  To  nucleate  brittle  fracture  in  80%  of  the 
computational  cells,  it  was  necessary  to  increase  the  plastic  shear  strain  threshold  to  0.45  for  the 
computations;  this  corresponds  to  an  equivalent  plastic  strain  of  0.78.  Average  homogeneous 
plastic  strain  values  of  0.21  and  0.24  were  obtained  for  the  HF-T  steel  cylinder  fragments  with  heat 
treatments  A  and  B,  respectively  (Table  18).  These  values  are  somewhat  smaller  than  the  equivalent 
plastic  strain  threshold  value  of  0.29  used  in  the  computations.  The  Armco  iron:  fragments  had 
about  75%  homogeneous  plastic  strain. 

The  computed  crack  distributions  in  Figures  45  and  46  were  used  to  calculate  fragment  mass 
distributions  for  the  cylinders  and  projectile.  All  computed  fragments  were  assumed  to  be  type  1. 
For  the  computations  using  the  BFRACT2  subroutine,  each  type  1  fragment  contained  the  inner 
and  outer  cylinder  (or  projectile)  surfaces,  two  fracture  surfaces  with  brittle  axial  cracks  perpen¬ 
dicular  to  the  cylinder  (or  projectile),  axis,  and  two  fracture  surfaces  with  brittle  circumferential 
cracks  perpendicular  to  cylinder  (or  projectile)  circumference.  For  the  HF-1  steel  cylinder  com¬ 
putations  using  the  BFRACT2  and;SHEAR2  subroutines  together,  each  fragment  had  two  addi¬ 
tional  fracture  surfaces  with  shear  cracks  on  planes  45°  between  the  radial  and  circumferential 
axes.  Ali  fragments  were  assumed  to  have  no  internal  cracks  and  to  be  formed  from  the  same 
number  of  surface  cracks.  In  order  to  use  all  the  cracks  to  form  fragments,  each  fragment  was  also 
assumed  to  have  the  same  ratio  of  the  number  of  cracks  of  each  orientation  as  determined  by  the 
crack  distribution  results.  Tables  0-2,  0-4,  0-6,  0-8,  0-10,  and  0-1 1  give  the  calculated  fragment 
mass  distribution  results.  The  detailed  procedure  for  obtaining  the  distributions  is  given  in  foot¬ 
notes  in  the  tables.  For  these  simulations,  an  iterative  procedure  was  used  to  select  the  number 
of  cracks  per  fragment.  This  parameter  was  varied  until  the.calculated  fragment  mass  distribution 
gave  the  best  agreement  with  experiment.  Figures  47,  48,  and  49  are  comparisons  of  .the 
calculations  and  experiments. 


FRAGMENT  MASS  M  (GRAIN) 

Figure  47.  Comparison  of  measured  and  calculated  fragment  mass  distributions  for  the  Armco 
iron  cylinder.  The  brittle  fracture  crack  distributions  in  Figure  45  were  used  to  obtain 
the  calculated  fragment  mass  distribution.  The  calculations  are  given  in  Table  0-2. 
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(a) 


FRAGMENT  MASS  M  (GRAIN) 


(b) 

Figure  48.  Comparison  of  measured  and  calculated  fragment  mass  distributions  for  the  HF-1 
steel  cylinders  with  heat  treatments  A  and  B.  (a)  The  brittle  fracture  crack 
distributions  in  Figure  45  were  used  to  obtain  the  calculated  fragment  mass 
distributions.  The  calculations  are  given  in  Tables  0-4  and  0-6.  (b)  The  brittle  and 
shear  fracture  crack  distributions  in  Figure  46  were  used  to  obtain  the  calculated 
fragment  mass  distributions.  The  calculations  are  given  in  Tables  0-8  and  0-10. 
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FRAGMENT  MASS  M  [GRAM) 


Figure  49.  Comparison  of  measured  and  calculated  fragment  mass  distributions  for  the  HF-1 
steel  cylinders  with  heat  treatment  A.  The  brittle  fracture  crack  distribution  in  Figure  45 
was  used  to  obtain  the  calculated  fragment  mass  distribution.  The  calculations  are 
given  in  Table  0-11. 

The  Armco  iron  computation  in  Figure  47  is  not  in  good  agreement  with  experiment.  Many 
more  cracks  were  activated  than  were  actually  needed  to  form  fragments.  Better  agreement  between 
computation  and  experiment  may  be  achieved  if  dynamic  fracture  parameters  could  be  determined 
from  framing  camera  photographs  that  show  the  growth  and  coalescence  of  exploding 
cylinder  surface  cracks  (see  Figure  3 1 ).  A  computation  using  the  BFRACT2  and  SHEAR2  sub¬ 
routines  together  was  also  performed  for  an  explosive-filled  Armco  iron  cylinder.  It  was  not  possible 
to  calculate  a  fragment  mass  distribution  for  this  case  since  type  1  fragments  could  not  be  formed 
because  the  total  computed  shear  band  crack  area  was  insufficient  to  form  fragments.  This  shear 
band  crack  area  was  only  about  0.8  and  0.2%  of  the  total  axial  and  circumferential  crack  areas, 
respectively. 

The  HF-1  steel  cylinder  computations  in  Figure  48  are  in  reasonable  agreement  with 
experiment.  About  twice  as  many  cracks  were  used  to  form  heat  treatment  A  fragments  than  heat 
treatment  B  fragments.  This  is  due  to  a  larger  average  brittle  crack  radius  for  heat  treatment  B 
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material  compared  to  heat  treatment  A  material  resulting  from  the  different  dynamic  fracture 
parameters  for  the  two  heat  treatments.  A  comparison  of  Figures  48(a)  and  48(b)  indicates  that  the 
use  of  the  SHEAR2  subroutine  in  the  PUFF  computations  had  a  negligible  affect  on  the  computed 
fragment  mass  distributions.  Appendix  P  gives  the  computed  stress,  particle  velocity,  and  crack 
concentration  versus  time  plots  for  each  computational  cell  for  the  HF-1  steel,  heat  treatment  B, 
cylindrical  PUFF  computation  using  the  BFRACT2  subroutine.  Input  parameters  and  output  results 
for  this  computation  arc  given  in  Figure  N-3.  This  appendix  is  included  so  that  selected  shock 
histories  can  be  observed  in  all  computional  cells  for  a  cylindrical  PUFF  computation. 

The  calculated  and  measured  HF-1  steel  projectile  fragment  mass  distributions  in 
Figure  49  are  in  reasonable  agreement  with  each  other.  Since  only  circumferential  cracks  were 
activated  for  this  computation,  an  axial  crack  distribution  equivalent  to  the  circumferential  crack 
distribution  was  assumed  in  order  to  form  fragments.  A  TROTT  computation  using  the  BFRACT2 
and  SHEAR2  subroutines  together  was  attempted  but  was  not  successful  due  to  some  unresolved 
computer  code  and  subroutine  incompatabilities. 

Figure  50  is  a  series  of  computer  plots  for  the  TROTT  simulation  of  the  HF-1  steel  projectile 
after  explosive  detonation.  Brittle  fracture  is  initiated  in  the  cells  nearest  the  outer  wall  and  pro¬ 
pagates  inward  in  agreement  with  proposed  fragmentation  models.7  Fracture  occurs  in  all  cells 
except  those  in  the  fuse  region.  Figure  43  shows  that  the  fuse  region  fragmented  into  a  few  tens 
of  mainly  large  type  1  fragments.  The  joint  region  of  the  projectile  was  not  simulated  in  the  com- 
k  tation.  Figure  51  is  a  framing  camera  photograph  of  an  AISI  3140  steel  projectile  approximately 
72  ps  after  explosive  detonation  that  is  shown  for  comparison  with  the  TROTT  simulation.  (A 
framing  camera  experiment  for  an  HF-1  steel  projectile  has  not  been  performed.)  The  AISI  3140 
material  has  an  Rc  40  hardness  and  a  1100-MPa  yield  strength.28  These  mechanical  properties  are 
similar  to  those  for  HF-1  steel  (Table  3).  The  maximum  expansion  of  the  AISI  3140  steel  projectile 
diameter  is  about  51$.  This  value  is  in  reasonable  agreement  with  a  maximum  computational 
expansion  of  58$  at  73  ns  after  explosive  detonation  for  the  HF-1  steel  projectile  simulation  in 
Figure  50. 
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Figure  50.  TROTT  simulation  of  the  exploding  HF-1  steel  projectile  at  selected  times.  The  CH6 
primary  explosive  in  the  fuse  is  detonated  at  0  jus  and  subsequently  detonates  the 
PBXN-106  secondary  explosive.  Prior  to  detonation  the  light  and  dark  cells  are  HF-1 
steel  and  explosive,  respectively.  The  dotted  cells  represent  reacted  explosive.  As  the 
detonation  front  propagates  from  right  to  left  in  the  projectile,  brittle  fracture  is 
initiated  in  HF-1  steel  cells  in  the  outer  wall  region  and  moves  inward.  The  dark  HF-1 
steel  cells  indicate  that  brittle  fracture  has  occurred  in  these  cells.  The  HF-1  steel 
cells  in  the  fuse  region  did  not  fracture  and  were  not  used  in  the  calculation  of  the 
fragment  mass  distribution. 
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Figure  51.  Framing  camera  photograph  of  AISI  3140  steel  projectile  approximately  78  ps  after 
explosive  detonation.  The  mechanical  properties  for  this  material  are  similar  to  those 
for  HF-1  steel  with  heat  treatment  A.  Cracks  cover  the  entire  surface  of  the  projectile 
in  agreement  with  the  computer  plot  at  73  ps  after  explosive  detonation  for  the  HF-1 
steel  projectile  simulation  in  Figure  50. 
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X.  SUMMARY  AND  RECOMMENDATIONS 
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The  fragmentation  of  Armco  iron  and  HF-1  steel  explosive-filled  cylinders  and  an  HF-1  steel 
projectile  has  been  simulated  with  the  SRI  PUFF  8  and  TROTT  computer  programs.  The 
BFRACT2  brittle  fracture  model  and  the  SHEAR2  shear  band  model  simulated  the  fracture 
processes.  Cylinder  computations  were  performed  for  two  heat  treatments  of  HF-1  steel.  Gas  gun 
impact  Hugoniot  and  soft  recovery  experiments  were  performed  for  both  HF-1  steel  heat 
treatments  to  determine  input  parameters  for  the  computations.  A  series  of  exploding  cylinder 
experiments  and  a  projectile  experiment  were  performed  to  provide  fragment  mass  distributions  for 
comparison  with  the  computations.  A  scheme  was  devised  for  characterizing  the  recovered 
fragments  from  the  experiments.  A  framing  camera  experiment  for  an  explosive-filled  Armco  iron 
cylinder  was  performed  in  which  the  initiation  and  growth  of  cylinder  surface  cracks  was  observed 
as  a  function  of  time.  This  type  of  brittle  crack  data  may  be  useful  for  the  direct  determination  of 
dynamic  fracture  parameters  for  use  in  exploding  cylinder  simulations. 


■w  -  -  umrii  •• 
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While  considerable  progress  has  been  made  in  the  computational  simulation  of  naturally 
fragmenting  munitions,  further  refinements  are  needed  before  a  detailed  computational  design  tool 
is  available.  For  example,  future  computational  efforts  should  consider  the  fragment  characterization 
results  of  Section  VI II  and  the  assumptions  in  Section  IX  that  lead  to  the  formation  of  fragments. 
Specifically,  attempts  should  be  made  to  calculate  fragment  mass  distributions  that  include  more 
than  one  fragment  type  (Figures  39  through  42)  and  to  have  fragments  that  contain  internal  cracks 
and  shear  bands  (Appendix  L).  The  mechanisms  of  axial  crack  formation  should  be  further  investi¬ 
gated  since  these  cracks  are  probably  formed  after  circumferential  crack  formation  (Figure  31). 
Analytical  techniques  should  also  be  extended  to  permit  the  use  of  the  BFRACT2  and  SHEAR2 
subroutines  together  in  two-dimensional  TROTT  computations. 
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Table  A-1.  Tensile-pull  measurements  for  transverse-direction  specimens  of  HF-1  steel.*,b 


Heat  Treatment  A  Heat  Treatment  B 

Experiment  1  Experiment  2  Experiment  3  Experiment  4 

Stress  Strain  Stress  Strain  Stress  Strain  Stress  Strain 

(MPa)  (pm/m)  (MPa)  (pm/m)  (MPa)  (pm/m)  (MPa)  (pm/m) 


0 

25.7 

51.3 
77.0 

102.6 

128.3 

153.9 
180.0 

205.3 

230.9 

256.7 

282.2 

307.9 

333.5 

359.2 

384.8 

410.5 

436.2 

461.8 

487.5 

513.1 

538.8 

564.4 

590.2 

615.8 

641.4 

667.1 

692.7 

718.4 
744.0 

769.7 

795.4 


0 

152 

313 

492 

660 

813 

962 

1094 

1236 

1368 

1494 

1633 

1776 

1902 

2036 

2180 

2304 

2440 

2562 

2684 

2814 

2944 

3076 

3200 

3334 

3460 

3590 

3733 

3874 

4028 

4160 

4310 


0 

25.6 

51.1 

76.7 

102.3 

127.9 

153.4 
179.0 

204.6 

230.2 

255.7 

281.3 

306.9 

332.5 
358.0 

383.6 

409.2 

434.8 

460.3 

485.9 

511.5 
537.0 

562.6 

588.2 

613.8 

639.3 

664.9 

690.5 

716.1 

741.6 

767.2 

792.8 


0 

80 

210 

318 

452 

580 

708 

840 

984 

1125 

1265 

1395 

1515 

1652 

1796 

1930 

2068 

2192 

2325 

2454 

2588 

2715 

2834 

2962 

3110 

3250 

3388 

3525 

3662 

3792 

3928 

4070 


0 

25.7 

51.3 
77.0 
102.6 

128.3 

153.9 
180.0 

205.3 

230.9 

256.6 

282.2 

307.9 

333.5 

359.2 

384.8 

410.5 

436.2 

461.8 

487.5 

513.1 

538.8 

564.4 

590.1 

615.8 

641.4 

667.1 

692.7 
744.0 

769.7 

795.4 
821  0 


0 

205 

312 

408 

500 

632 

766 

905 

1044 

1168 

1286 

1414 

1530 

1670 

1800 

1925 

2054 

2195 

2340 

2500 

2638 

2790 

2950 

3090 

3240 

3404 

3590 

3770 

4280 

4655 

5150 

5675 


0 

25.6 

51.2 

76.8 

102.4 
128.0 

153.6 

179.2 

204.8 

230.4 
256.0 

281.6 

307.2 

332.8 

358.4 
384.0 

409.6 

435.2 

460.8 

486.4 
512.0 

537.6 

563.2 

588.8 

614.4 
640.0 

665.6 

691.2 

716.8 

742.4 
768.0 

793.6 


0 

100 

235 

356 

496 

642 

795 

945 

1118 

1292 

1458 

1595 

1730 

1875 

2045 

2202 

2350 

2490 

2630 

2788 

2940 

.3092 

3235 

3390 

3536 

3695 

3874 

4064 

4284 

4555 

4900 

5395 
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Table  A-1.  Tensile-pull  measurements  for  transverse-direction  specimens  of  HF-1  steel.*’b  (Cont'd) 


Heat  Treatment  A 

Heat  Treatment  B 

Experiment  1 

Experiment  2 

Experiment  3 

Experiment  4 

Stress 

Strain 

Stress 

Strain 

Stress 

Strain 

Stress 

Strain 

(MPa) 

(pm/m) 

(MPa) 

(pm/m) 

(MPa) 

(pm/m) 

(MPa) 

(pm/m) 

821.0 

4452 

818.4 

4210 

846.7 

6225 

819.2 

5940 

846.7 

4642 

843.9 

4348 

872.3 

6850 

844.8 

6555 

872.3 

4886 

869.5 

4465 

885.2C 

7170° 

870.4 

7295 

898.0 

5268 

895.1 

4650 

879.4C 

7525c 

923.6 

5786 

920.6 

4844 

949.3 

7110 

946.0 

5084 

975.0 

11700 

971.8 

5416 

987.8 

12960 

984.6 

5710 

1000.6 

14790 

997.4 

6090 

1013.4 

16960 

1010.2 

6858 

1026.3 

18830 

1022.9 

8670 

1039.1 

21020 

1035.7 

11200 

1051.9 

22500 

1048.5 

13370 

1064.7 

23980 

1061.3 

15^20 

1090.4C 

— 

1074.1 

17382 

1086.9  19360 

1099.7  21300 

1112.5  23180 

1 120.1 c  25  325c 


a  The  transverse  tensile  specimens  vcrc  fabricated  from  6.35-mm-thick  disks  of  HF-1  steel  (the  same  group  of  disks  that  was  used 
for  the  gas  gun  experiments).  The  disks  were  cut  from  a  bar  of  HF-1  steel.  The  cylindrical  axis  of  the  specimens  was  transverse  to 
the  axis  of  the  bar.  The  specimens  were  dumbcll  shaped  with  a  central-region  diameter  of  4.70  mm.  For  experiments  1  and  3,  the 
total  length  and  central-region  length  of  the  specimens  was  35.6  and  19.1  mm,  respectively.  For  experiments  2  and  4,  the  total 
length  and  central-region  length  of  the  specimens  was  27.9  and  12.7  mm,  respectively. 

*’T1ic  specimens  were  pulled  in  a  tensile-testing  machine  at  a  constant  crosshead  speed.  The  strain  measurements  were  made  with  a 
portable  strain  indicator  (Vishay  Instruments,  Inc.  Model  P-350).  Strain  gauges  (Micro-Measurements  Gage  Type 
FA-06-1 25BT-1 20, 120.0  i  0.3%  O  resistance,  2,11  ±  0.5%  gage  factor,  3.18-mm  long,  1.57-mm  grid  width,  were  used  for 
experiments  1  and  3;  Micro-Measurements  Gage  Type  EA-06-062AP-120,  120.0  i  0.3%  n  resistance,  2.085  i  0.5%  gage  factor, 
1.57-mm  long,  1.57-mm  grid  width,  were  used  for  experiments  2  and  4)  were  attached  to  the  central  region  of  the  specimens  with 
epoxy. 

c  Stress  and  strain  at  fracture.  For  experiment  1  the  strain  of  fracture  was  not  recorded. 
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APPENDIX  B 

OSCILLOSCOPE  RECORDS  FOR  HF-1  STEEL  HUGONIOT  EXPERIMENTS 


(Figures  B-l,  B-2,  and  B-3  are  the  oscilloscope  records  for  heat  treatment  A  material.  Figures  B-4 
and  B-5  are  oscilloscope  records  for  heat  treatment  B  material.  All  the  experiments  were  shock 
transit  time  experiments  except  Shot  120  (Figure  A-3),  which  was  a  direct  impact  experiment.) 


B-l 


(b) 


Figure  B-1.  Oscilloscope  records  for  Shot  117.  Time  increases  from  left  to  right,  (a)  Tilt  data  re¬ 
cord.  The  vertical  scale  is  4  V/div  and  the  horizontal  scale  is  50  ns/div.  The  middle 
trace  is  an  initial  time  reference  square  pulse.  The  lower  trace  is  a  10-ns-period  time 
calibration  signal,  (b)  Current  record  from  back-surface  quartz  gauge.  The  vertical 
scale  is  100  mA/div  and  the  horizontal  scale  is  50  ns/div.  The  middle  trace  is  a  time 
reference  square  pulse  delayed  by  1120  ns  with  respect  to  the  initial  time  reference 
square  pulse  in  (a)  for  wave  velocity  measurements.  The  lower  trace  is  a  10-ns-period 
time  calibration  signal.  Only  the  beginning  of  the  plastic  wave  was  recorded  for  this 
shot. 
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Figure  B-2.  Oscilloscope  records  for  Shot  118.  Time  increases  from  left  to  right,  (a)  Tilt  data  re¬ 
cord.  The  vertical  scale  is  4  V/div  and  the  horizontal  scale  is  50  ns/div.  The  middle 
trace  is  an  initial  time  reference  square  pulse.  The  lower  trace  is  a  10-ns-period  time 
calibration  signal,  (b)  Current  record  from  back-surface  quartz  gauge.  The  vertical 
scale  is  100  mA/div  and  the  horizontal  scale  is  50  ns/div.  The  middle  trace  is  a  time 
reference  square  pulse  delayed  by  535  ns  with  respect  to  the  initial  time  reference 
square  pulse  in  (a)  for  wave  velocity  measurements.  The  lower  trace  is  a  10-ns-period 
time  calibration  signal. 
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Figure  B-3.  Oscilloscope  record  for  Shot  120.  Current  record  from  direct  impact  onto  a  quartz 
gauge.  Time  increases  from  left  to  right.  The  vertical  scale  is  100  mA/div  and  the 
horizontal  scale  is  50  ns/div.  The  lower  trace  is  a  10-ns-period  time  calibration  signal. 


Figure  B-4.  Oscilloscope  records  for  Shot  152.  Time  increases  from  left  to  right,  (a)  Tilt  data  re¬ 
cord.  The  vertical  scale  is  4  V/div  and  the  horizontal  scale  is  50  ns/div.  The  middle 
trace  is  an  initial  time  reference  square  pulse.  Thedower  trace  is  a  10-ns-period  time 
calibration  signal,  (b)  Current  record  from  back-surface  quartz  gauge.  The  current 
calibration  trace  (upper  horizontal  line)  has  an  amplitude  of  100  mA.  The  horizontal 
scale  is  100  ns/div.  The  middle  trace  is  a  time  reference  square  pulse  delayed  by 
1140  ns  with  respect  to  the  initial  time  reference  square  pulse  in  (a)  for  wave  velocity 
measurements.  The  lower  trace  is  a  20-ns-period  time  calibration  signal. 
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Figure  B-5.  Oscilloscope  records  for  Shot  153.  Time  increases  from  left  to  right,  (a)  Tilt  data  re¬ 
cord.  The  vertical  scale  is  4  V/div  and  the  horizontal  scale  is  50  ns/div.  The  middle 
trace  is  an  initial  time  reference  square  pulse.  The  lower  trace  is  a  10-ns-period  time 
calibration  signal,  (b)  Current  record  from  back-surface  quartz  gauge.  The  current 
calibration  trace  (upper  horizontal  line)  has  an  amplitude  of  100  mA.  The  horizontal 
scale  is  100  ns/div.  The  middle  trace  is  a  time  reference  square  pulse  delayed  by 
1180  ns  with  respect  to  the  initial  time  reference  square  pulse  in  (a)  for  wave  velocity 
measurements.  The  lower  trace  is  a  20-ns-period  time  calibration  signal. 
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APPENDIX  C 

EXPRESSIONS  FOR  TRANSFORMING  THE  HUGONIOT  EQUATION  OF 
STATE  FROM  A  LINEAR  SHOCK  VELOCITY-PARTICLE  VELOCITY 
RELATIONSHIP  INTO  A  STRESS-COMPRESSION  RELATIONSHIP 


e-i 


In  this  appendix,  expressions  are  derived  for  converting  the  Hugoniot  equation  of  state 
from  a  stress-particle  velocity  (a„-up)  relationship  into  a  stress-compression  (0,,-m)  relation¬ 
ship.  Here  the  compression 


P  = 


(C-l) 


where  VQ  and  V  are  initial  and  final  specific  volumes,  respectively. 

In  the  clastic  region  the  momentum  (oH-up)  and  mass  conservation  equations0-1  are, 
respectively, 


and 


aH  =  P 


o  *a  up» 


(C-2) 


(C-3) 


where  pQ  =  1/V0  is  the  initial  density,  CL  is  the  longitudinal  wave  velocity,  and  uc  is  the  particle 
velocity  at  the  elastic  limit.  Eliminating  V/V0  from  Equations  (C-l )  and  (C-3)  gives 


Up  _  P 

cj  =  I+m’ 


P  <  Pc' 


(C-4) 


where  pc  is  the  compression  at  the  elastic  limit.  Substituting  this  equation  into  Equation  (C-2) 
gives 


=  177’  (C'5) 

for  the  stress-compression  relationship. 

This  equation  can  be  transformed  further  by  using  expressions  for  the  adiabatic  bulk 
modulus  K  and  shear  modulus  G  for  an  isotropic  material, 


’m-m.~i.rn  — - . .  ,  ,  u 


u 


C-3 


r 


^  =  ^0  ~  P0  “  T  ^  * 


(C-6) 


and 


G  P0C s  , 


(C-7) 


where  C'1(  and  Cs  are  the  bulk  and  shear  wave  velocities,  respectively.  Substituting  these  equations 
into  Equation  (C-5)  gives 


4  u  u 

a„  =  -G  —  +  K-i— ,  i u  <  u 
11  3  l+/i  1  +u  c 


(C-8) 


In  this  form  of  the  stress-compression  relationship  the  stress  is  the  sum  of  shear  and  hydrostatic 
components.  If  the  deviator  stress  c‘*  oD  is  defined  as 


„  4  r  ^ 
0  ‘  3° 


(C-9) 


then  at  the  elastic  limit  (p  =  pe)  o(|  becomes 


o  =  a.,  +  K  - 

c  D  1  +  P, 


The  deviator  stress  is  a  measure  of  the  ability  of  a  material  to  support  shear. 


(C-10) 


For  the  plastic  region,  an  ideal  elastic-plastic  wave  structure  is  assumed  with  equilibrium 
initial  and  final  states.  The  momentum  and  mass  conservation  relations  C’1  become,  respectively, 


°H  =  °c+Pc  <Us  “  UcXUn  ~  Uc>’  »f>  >  V 


and 


V 


=  1- 


Us  -  up 


C1  Aus“  Uc 


Up>Ue. 


(C-ll) 


(C-12) 


where  =  pQ  ( I  +  pe)  is  the  density  at  the  elastic  limit  and  Us  is  the  shock  velocity.  Substituting 
V^/Vp  =  1  -  ue/CL  from  Equation  (C-3)  into  Equation  (C-12)  and  simplifying  gives 


C-4 


-  rr- 


3 


O 


i 


I 


V  up-uc 

UP>U«  ■ 

c  e 

Using  V/Ve  =  (1  +  juc)/(l  +  p)  from  Equation  (C-l)  in  Equation  (C-13)  gives 


(C-13) 


U,,  -  Uc  H-  He 

\  ~  Ue  = 


,  P>PC 


Inserting  up  -  u  from  this  equation  into  Equation  (C-l  1)  gives 


(C-l  4) 


°H  =  ae+MUi“  ue>2X  .  P>Pe 


(C-l  5) 


where 


P“  Pc  ,  /  P  Pe 

X  =  T—  =  (1  +PC)  (  77 - — ~ 

1  +  p  c  l  1  +  P  1  +  /r 


(C-l  6) 


Equations  (C-l 5)  can  be  reduced  further  by  assuming  an  analytical  form  for  the  Us~up 
relationship.  A  linear  relationship 


Us  =  C0  +  sup 


(C-l  7) 


is  taken  because  it  has  been  found  that  this  form  represents  the  Hugoniot  for  many  materials.  c'2 
The  parameter  C0  for  many  materials  is  usually  very  nearly  equal  to  CB ,  the  bulk  sound  velocity 
at  zero  stress.  c'3  Substituting  Equation  (C-l 7)  into  Equation  (C- 15)  and  using  Equation  (C-l 0) 
and  the  equation  C0  -  uc  +  sup  =  (CQ  -  ug  +  sue )/( 1  -  sx)  gives 


Me  X 

o,i  -  an  +  K - +  p„  (Cn  -  u  +  su  )2  - : 

li  D  1+M  ''e  v  0  e  e'  (,  _  s%y 


2  -  P>PC- 


(C-18) 


for  the  stress-compression  relationship  in  the  plastic  region.  At  the  elastic  limit,  Equation  (C-l 8) 
becomes an  =  a,  since  x  =  0  at  p  =  pe. 


For  ix  >  pe  the  deviatoric  stress  aD  can  be  defined  as  the  difference  between  the  stress- 
compression  relationship  and  the  hydrostat.  If  the  material  strain  hardens,  the  deviatoric  stress 
will  be  a  function  of  p.c"' 
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t. 


Equation  (C-18)  can  be  expressed  as  a  series  expansion  in  terms  of  the  compression  p. 
The  term  x/(i  -  sx)2  in  Equation  (C-l  8)  can  be  written  as 

X  _  ( I  +  aO(M ~ 


(I  -  sx)‘  (1  +s/ie  -  (s-  l)p)' 


Expanding  the  denominator  gives 


(C-l  9) 


X  _ -pc  +  (l  +  ^  ^  0 

— —  __  A  1  T  2 

_  »  «  \  2  e  i  i  \L 


2  /  s  -  1  \  /  s  -  1 


(1  -  sx)' 


(1  +S 


s-  1  V 


[ i 3  +  5 


s-  1  V 


M  + ...  . 


Regrouping  terms  gives 


(C-20) 


0  -lx)2  =  [-  A*c  +  l-2p^e  +  (I  -  +  l-  3P2Me  +  2p(l  -  /ic)  +  1  )m2 


+  [- 4p2juc +  3p(l  - /ie)  +  2|pM3  +  [- 5p2Me +  4p(l  - /ae)  +  3ipV +...  J  ,  (C-2 1 ) 

where  p  -  (s  -  I )/( 1  +  s/ij.  Substituting  Equation  (C-2 1 )  into  Equation  (C-18)  and  taking 
K  =  p0  Cq  gives  after  reducing 

/  [i-M.(-L-i)+sMc-L-T\ 


(1  +spc)2 


(1  -  SMe) 

(1  +spe) 


[(2s-  1)(1  -  s/ic)  +  sjuc(s+ 1)]  ,  (s-  l)[(3s-  1)(1  -  sm0)  +  2spc(s+ 1)] 

(l+s/O2  11  (I  +  sp  )3  11 


(s  -  I  )2 [(4s  -  DM  ••  spc)  +  3sjuc (s  +  1 )] 


(i  f-sMcr 


■/i  + . 


^  r.. 


(C-22) 


Equation  (C-22)  can  be  simplified  because  for  most  materials  pc  <<  1,  s  =2  1.5,  and  CL/C0 
s  l.c'2  Therefore,  neglecting  thepc  terms  gives 


a„  =  oD  +  Km  +  K(2s -  1  )p2  +  K(s -  l)(3s-  1  )m3  +  K(s -  l)2(4s-  1)m4+...,  P  >  Pe-  (C-23) 

This  equation  gives  a  simpler  form  for  the  stress-compression  relationship  than  Equation  (C-22). 
Only  the  initial  density  pQ  and  the  two  parameters  in  the  linear  shock  velocity-particle  relationship 
U$  =  C„  +  su()  are  necessary  to  determine  the  coefficients.  For  HF-1  steel  the  calculated 
coefficients  for  this  equation  only  differed  by  a  few  percent  from  the  more  complicated 
coefficients  in  Equation  (C-22).  For  most  materials  the  coefficients  in  Equation  (C-23)  are 
evaluated  through  the /a3  term.  Christman  et  al.0'1  calculated  a  zero-pressure  Gruneisen  parameter 
70  =  2s- 1  using  the  Dugdale  and  MacDonald  Gruneisen  formula.0'4  The  cofficient  of  the  p2  term 
in  Equation  (C-23)  can  therefore  be  represented  by  f<7o  where  7o  is  a  zero-pressure  Gruneisen 
parameter  estimate  and  K  is  the  low-pressure  bulk  modulus.  This  observation  was  initially  made  by 
Rice  et  al.c'5 

The  results  of  this  section  could  have  been  expressed  in  terms  of  the  compression  17  where 

V  ~  1  •  (C-24) 

vo 

Using  Equation  (C-l)  for  the  compression  p,  it  can  be  shown  that  (1  +  p)(l  -  r/)  =  1,  which  gives 
the  relations  ij  =  p/{  1  +  p)  and  p  =  i 7/(1  -  77). 
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APPENDIX  D 


PHOTOGRAPHS  OF  THE  SPALL  FRACTURE  SPECIMENS  OF  HF-1  STEEL 

(The  impactor  velocity,  initial  impactor  thickness,  and  initial  specimen  thickness  are  given  for 
each  shot.  Cracks  were  digitized  for  all  the  shots  except  the  no-damage  Shots  89  and  147.  Cracks 
were  measured  in  the  approximate  central  one-third  region  of  the  sectioned  specimen  face.  The 
bottom  edge  of  each  specimen  is  the  impact  surface.  The  specimen  magnification  is  approximately 
5X). 


(a)  Shot  95,  0.130-km/s  impactor  velocity,  2.37-mm-thick  impactor,  6.31-mm-thick  specimen. 


(b)  Shot  110,  0.149-km/s  impactor  velocity,  2.37-mm-thick  impactor,  6.35-mm-thick  specimen. 


4  mm 


(c)  Shot  98,  0.160-km/s  impactor  velocity,  2.37-mm-thick  impactor,  6.35-mm-thick  specimen. 


Figure  D-1.  Sectioned  HF-1  steel  specimens  with  heat  treatment  A  for  (a)  Shot  95,  (b)  Shot  110, 
and  (c)  Shot  98. 


D-3 


(a)  Shot  92,  0.180-km/s  impactor  velocity,  2.37-mm-thick  impactor,  6.34-mm-thick  specimen. 


(b)  Shot  91,  0.183-km/s  impactor  velocity,  2.37-mm-thick  impactor,  6.34-mm-thick  specimen. 


Figure  D-2.  Sectioned  HF-1  steel  specimens  with  heat  treatment  A  for  (a)  Shot  92,  (b)  Shot  91, 
and  (c)  Shot  90. 
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(a)  Shot  94,  0.200-  <m/s  impactor  velocity,  2.37-mm-thick  impactor,  6.36-mm-thick  specimen. 


(b)  Shot  89,  0.120-km/s  impactor  velocity,  1.61-mm-thick  impactor,  3.18-mm-thick  specimen. 


Figure  D-3.  Sectioned  HF-1  steel  specimens  with  heat  treatment  A  for  (a)  Shot  94,  (b)  Shot  89, 
and  (c)  Shot  113. 


I 


(a)  Shot  143,  0.169-km/s  impactor  velocity,  2.37-mm-thick  impactor,  6.35-mm-thick  specimen. 


(b)  Shot  142,  0.191-km/s  impactor  velocity,  2.37-mm-thick  impactor,  6.36-mm-thick  specimen. 


4  mm 


(c)  Shot  150,  0.152-km/s  impactor  velocity,  1.58-mm-thick  impactor,  6.36-mm-thick  specimen. 


Figure  D-6.  Sectioned  HF-1  steel  specimens  with  heat  treatment  B  for  (a)  Shot  143,  (b)  Shot  142, 
and  (c)  Shot  150. 
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(c)  Shot  147,  0.147-km/s  impactor  velocity,  1.15-mm-thick  impactor,  3.18-mm-thick  specimen. 


Figure  D-7.  Sectioned  HF-1  steel  ioecimens  with  heat  treatment  D  for  (a)  Shot  148,  (b)  Shot  149, 
and  (c)  Shot  147. 
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(a)  Shot  146,  0.167-km/s  impactor  velocity,  1.16-mm-thick  impactor,  3.18-mm-thick  specimen 


4  mm 


(b)  Shot.  151,  0.190-km/s  impactor  velocity,  1.15-mm-thick-impactor,  6.36-mm-thick  specimen. 


Figure  D-8.  Sectioned  HF-1  steel  specimens  with  heat  treatment  B  for  (a)  Shot  146  and 
(b)  Shot  151. 
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APPENDIX  E 


CRACK  DISTRIBUTIONS  FOR  SPALL  FRACTURE  SPECIMENS  OF  HF-1  STEEL 


(Figures  are  presented  for  eleven  heat  treatment  A  specimens  (Shots  95  through  1 12)  and  ten  heat 
treatment  B  specimens  (Shots  154  through  151).  Figures  are  not  presented  for  the  no-damage 
Shots  S9  and  147.  Eacli  figure  contains  five  plots.  Plots  (a)  and  (b)  are  distributions  for  the 
digitized  cracks.  Cracks  were  measured  in  the  approximate  central  one-third  region  of  the  sectioned 
specimen  face.  A  crack  was  considered  to  be  in  the  central  region  if  its  midpoint  was  in  the  region. 
Plots  (c),  (d),  and  (e)  are  crack-size  distributions  that  were  calculated  from  the  digitized  cracks 
using  the  SRI  statistical  transformation  computer  program  BABS3.  The  dashed  line  in  plot  (d)  of 
each  figure  is  a  fit  of  the  equation  Ng  =  Nfl  exp(-  R/Rj )  to  the  volume  crack-size  distribution  for 
the  zone  of  maximum  damage.  Here  Ng  is  the  number  of  cracks  per  unit  volume  with  radius 
greater  than  Rj .  The  calculated  values  for  NQ  and  Rj  for  the  zone  of  maximum  damage  are  given 
in  the  plot.  The  spall  plane  in  plots  (b)  and  (e)  is  taken  as  the  difference  in  the  specimen  and 
impactor  thicknesses.) 
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Figure  E-4.  (Continued) 
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Figure  E-7.  (Continued) 
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Figure  E-7.  (Continued) 
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Figure  E-8.  (Continued) 
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N0=  7.92X104  CRftCKS/Ctt3 
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Figure  E-9.  Crack  distributions  for  Shot  88.  The  cracks  from 
only  one  half  of  this  full-spalled  specimen  were 
digitized.  This  shot  was  not  used  for  calculating 
dynamic  fracture  parameters  for  heat  treatment 
A  material,  (a)  Digitized  cracks  on  sectioned 
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Figure  E-10.  (Continued) 
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Figure  E-17.  (Continued) 
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Figure  E-20.  (Continued) 
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APPENDIX  F 

DIGITIZED  CRACK  COORDINATES  FOR  THE  SPALL  FRACTURE 
SPECIMENS  OF  HF-1  STEEL 


(Digitized  crack  coordinates  are  presented  for  all  the  counted  cracks  for  heat  treatments  A  and  B 
of  HF-1  steel.  Plot  (a)  of  each  figure  in  Appendix  E  shows  the  digitized  crack  distribution  for 
each  shot.  The  upper  left  corner  of  each  plot  is  the  coordinate  system  orgin  for  the  digitized  crack 
values.  Four  values  are  listed  for  each  crack  in  row  form.  The  initial  two  values  are  the  coordinates 
of  one  endpoint  position  of  a  crack  measured  in  the  specimen  length  and  thickness  directions, 
respectively.  The  final  two  values  are  the  coordinates  of  the  other  endpoint  position  of  the  crack 
measured  in  the  specimen  length  and  thickness  directions,  respectively.  These  values  are  given  in 
inches.  The  cracks  given  in  the  eight  columns  on  the  left  side  of  a  page  are  continued  at  the  top  of 
the  eight  columns  on  the  right  side.  The  specimen  length  and  thickness  magnification  factors  are 
given  in  Table  7.  The  cracks  were  measured  on  enlarged  photographs  of  the  sectioned  specimen 
surfaces.  The  heat  treatment  A  stringer  cracks  (three  for  Shot  98,  four  for  Shot  92,  two  for  Shot 
91,  three  for  Shot  90,  fourteen  for  Shot  94,  nine  for  Shot  1 13,  and  nine  for  Shot  112)  are  listed 
last  for  the  shots  that  contain  them.) 
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Figure  F-1.  Listing  of  digitized  crack  coordinates  for  the  spall 

fracture  specimens  of  HF-1  steel.  Figure  F-1.  (Continued) 
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APPENDIX  G 

COMPUTER  LISTINGS  OF  THE  INPUT  PARAMETERS 
USED  IN  THE  PUFF  COMPUTATIONS  FOR  THE 
HF-1  STEEL  SPALL  FRACTURE  EXPERIMENTS 

(The  nucleation  and  growth  values  in  Tables  9  and  10  were  obtained  from  PUFF  computations 
using  the  input  parameters  in  this  appendix.  Figures  G-l  through  G-5(a)  give  the  input  parameters 
for  the  computations  for  the  heat  treatment  A  experiments.  Figures  G-5(b)  through  G-IO  give  the 
input  parameters  for  the  computations  for  the  heat  treatment  B  experiments.  The  input 
parameters  for  the  BFRACT2  fracture  model  are  given  in  the  TSR1  and  TSR2  lines  of  each  listing. 
The  HF-1  steel  equation  of  state  data  is  given  in  the  EQST  and  Y0  lines.  An  explanation  of  the 
these  parameters  is  giver/  in  Section  VI  and  Reference  1 .) 
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Figure  G-1.  PUFF  input  parameters  for  (a)  Shot  110  and  (b)  Shot  98. 


Figure  G-2.  PUFF  input  parameters  for  (a)  Shot  92  and  (b)  Shot  91 
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Figure  G-3.  PUFF  input  parameters  for  (a)  Shot  90  and  (b)  Shot  94. 
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Figure  G-5.  PUFF  input  parameters  for  (a)  Shot  1 12  and  (b)  Shot  154. 
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Figure  G-6.  PUFF  input  parameters  for  (a)  Shot  145  and  (b)  Shot  144, 
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Figure  G-8.  PUFF  input  parameters  for  (a)  Shot  150  and  (b)  Shot  148. 


input  parameters  for  (a)  Shot  149  and  (b)  Shot  146. 


Figure  G-10.  PUFF  input  parameters  for  Shot  151, 


APPENDIX  H 

COMPUTED  STRESS,  PARTICLE  VELOCITY,  AND  CRACK  CONCENTRATION 
VERSUS  TIME  PLOTS  FOR  THE  CELL  OF  MAXIMUM  DAMAGE  FOR  THE 
HF-1  STEEL  SPALL  FRACTURE  SPECIMENS 

(The  input  parameters  in  Appendix  G  (with  the  appropriate  JEDIT  parameters)  were  used  in 
PUFF  8  computations  to  produce  these  plots.  No  plots  are  presented  for  the  threshold-damage 
Shot  95,  the  no-damage  Shots  89  and  147,  and  the  full-spall  Shot  88.  Tables  9  and  10  summarize 
some  of  the  stress  and  crack  density  values.  The  stress  and  crack  concentration  values  are 
computed  at  the  cell  midpoints.  The  particle  velocity  values  are  computed  at  the  cell  edge  nearest 
the  impact  surface.  Each  specimen  was  divided  into  20  computational  cells  of  equal  width, 
numbered  10  through  29  beginning  at  the  specimen  impact  surface.  Plots  (a)  and  (b)  refer  to 
principal  stresses  in  the  direction  of  shock  wave  propagation  and  orthogonal  to  that  direction, 
respectively.  Compressive  and  tensile  stress  values  are  positive  and  negative,  respectively.  Plots  (c) 
and  (d)  give  the  particle  velocity  in  the  shock  propagation  direction  and  the  generated  crack 
concentration  for  cracks  oriented  in  planes  perpendicular  to  that  direction,  respectively.) 
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Computed  shock  histories  in  cell  23  for  Shot  110.  (a)  Parallel  stress,  {b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 


Figure  H-2.  Computed  shock  hisi 
(d)  Crack  concentrati 
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Figure  H-3.  Computed  shock  histories  in  cell  23  for  Shot  92.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Partic  levelocity 
(d)  Crack  concentration. 


Figure  H-4.  Computed  shock  histories  in  cell  23  for  Shot  91.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 
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Computed  shock  histories  in  cell  23  for  Shot  90.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 
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Figure  H-6.  Computed  shock  histories  in  cell  23  for  Shot  94.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 
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Figure  H-7.  Computed  shock  histories  in  cell  20  for  Shot  113.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 


Computed  shock  histories  in  cell  23  for  Shot  112.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 
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Computed  shock  histories  in  cell  23  for  Shot  144.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 
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Figure  H-13.  Computed  shock  histories  in  cell  23  for  Shot  143.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 


Computed  shock  histories  in  cell  23  for  Shot  142.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 


!.  Computed  shock  histories  in  ceil  25  for  Shot  150.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 
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Computed  shock  histories  in  ceil  23  for  Shot  146.  (a)  Parallel  stress,  (b)  Perpendicular  stress,  (c)  Particle  velocity, 
(d)  Crack  concentration. 


Computed  shock  histories  in  cell  26  for  Shot  151.  (a)  Parallel  stress,  (b)  Perpendiclar  stress,  (c)  Particle  velocity 
(d)  Crack  concentration. 


APPENDIX  I 

COMPUTED  STRESS,  PARTICLE  VELOCITY,  AND  CRACK 
CONCENTRATION  VERSUS  TIME  PLOTS  FOR  EACH 
COMPUTATIONAL  CELL  FOR  THE  HF-1  STEEL,  HEAT 
TREATMENT  B,  SPALL  FRACTURE  SHOT  144 

(For  Shot  144,  a  2.37-mm-thick  disk  impacted  a  6.35-mm-thick  specimen  at  0.148  km/s.  The 
input  parameters  in  Figure  G-6  (with  the  appropriate  JEDIT  parameters)  were  used  in  PUFF  8 
computations  to  produce  these  plots.  The  impactor  was  divided  into  eight  computational  cells  of 
equal  width,  numbered  1  through  8  beginning  at  the  impactor  free  surface.  The  specimen  was 
divided  into  20  computational  cells  of  equal  width,  numbered  10  through  29  beginning  at  the 
specimen  impact  surface.  Plots  (a)  and  (b)  of  Figures  1-1  through  1-28  give  the  principal  stresses  in 
the  shock  wave  propagation  direction  and  orthogonal  to  that  direction,  respectively.  Compressive 
and  tensile  stress  values  are  positive  and  negative,  respectively.  The  stresses  are  computed  at  the 
cell  midpoints.  Figures  1-29  through  142  give  the  particle  velocity  in  each  cell  in  the  shock  wave 
propagation  direction.  The  particle  velocity  values  are  computed  at  the  cell  edge  nearest  the 
impact  surface.  Figures  1-43  through  1-49  give  crack  concentrations  in  those  computational  cells  in 
which  tensile  fracture  occurred.  Crack  concentrations  are  computed  at  the  cell  midpoints  and  refer 
to  cracks  oriented  in  planes  perpendicular  to  the  shock  wave  propagation  direction.) 
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histories  in  cell  3  for  Shot  144.  Figure  1-4.  Computed  stress  histories  in  cell  4  for  Snot  144. 
1  jb)  perpendicular  stresses.  (a)  Parallel  and  (b)  perpendicular  stresses. 
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histories  in  cell  17  for  Shot  144. 
id  (b)  perpendicular  stresses. 
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histories  in  cell  18  for  Shot  144.  Figure  1-18.  Computed  stress  histories  in  cell  19  for  Shot  144, 
id  (b)  perpendicular  stresses.  (a)  Parailel  and  (b)  perpendicular  stresses. 


Computed  stress  histories  in  cel!  22  for  Shot  144.  Figure  1-22.  Computed  stress  histories  in  cell  23  for  Shot  144. 
(a)  Parallel  and  (b)  perpendicular  stresses.  (a)  Parallel  and  (b)  perpendicular  stresses. 


•e  1-25.  Computed  stress  histories  in  cell  26  for  Shot  144.  Figure  1-26.  Computed  stress  histories  in  cell  27  for  Shot  144. 
(a)  Parallel  and  (b)  perpendicular  stresses.  (a)  Parallel  and  (b)  perpendicular  stresses. 
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1-32.  Computed  particle  velocity  histories  in  (a)  cell  7 
and  (b)  cell  8  for  Shot  144. 
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-33.  Computed  particle  velocity  histories  in  (a)  cell  10  Figure  1-34.  Computed  particle  velocity  histories 
and  (b)  cell  1 1  for  Shot  144.  and  (b)  cell  13  for  Shot  144. 
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Shot  144. 
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Computed  particle  velocity  histories 
and  (b)  cell  25  for  Shot  144. 


Computed  particle  velocity  histories  in  (a)  cell  26  Figure  1-42.  Computed 
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Computed  crack  concentration  histories  in  Figure  1-46.  Computed  crack  concentration  histories  in 
(a)  cell  18  and  (b)  cell  19  for  Shot  144.  (a)  cell  20  and  (b)  cell  21  for  Shot  144. 


Figure  l-49.Computed  crack  concentration  histories  in 
(a)  cell  26  and  (b)  cell  27  for  Shot  144. 


APPENDIX  J 


COMPUTATION  OF  DYNAMIC  FRACTURE  PARAMETERS  FOR  HEAT  TREATMENT  A 
OF  HF-1  STEEL  WITH  STRINGER  CRACKS  EXCLUDED 

(To  determine  the  importance  of  the  stringer  cracks,  dynamic  fracture  parameters  were  calculated 
for  heat  treatment  A  material  with  these  cracks  excluded  from  the  distributions.  Stringer  cracks 
are  those  cracks  that  occurred  along  the  stringer  impurities  in  the  material.  Stringer  impurities 
are  elongated  impurity  regions  produced  during  the  material  fabrication  process.  For  the  HF-1 
steel  specimens  the  impurities  were  parallel  to  the  shock  wave  propagation  direction.  The  stringer 
cracks  were  therefore  perpendicular  to  the  majority  of  the  cracks  in  a  specimen.  A  total  of  44 
stringer  cracks  were  measured:  three  for  Shot  98,  four  for  Shot  92,  two  for  Shot  91,  three  for 
Shot  90,  fourteen  for  Shot  94,  nine  for  Shot  1 13,  and  nine  for  Shot  112.  The  digitized  stringer 
crack  coordinates  are  listed  in  Appendix  F.  Figures  J-l  through  J-7  give  the  crack  distributions  for 
these  seven  spall  fracture  specimens  with  stringer  crack  excluded.  These  distributions  can  be 
compared  with  those  that  include  all  the  cracks  given  in  Appendix  E.  Table  J-l  summarizes  the 
experimental  volume  crack-size  distribution  parameters  for  the  heat  treatment  A  specimens 
with  stringer  cracks  excluded  .  Six  iterations  with  the  PUFF  8  computer  program  using  the 
BFRACT2  subroutine  were  performed  to  obtain  the  dynamic  fracture  parameters.  Table  J-2 
gives  the  experimental  and  computed  nucleation  and  growth  values  obtained  using  the  final  set  of 
parameters.  These  values  are  plotted  in  Figure  J-8.  The  experimental  and  computed  intercept  and 
slope  values  for  the  nucleation  rate  curves  in  Figure  J-8(a)  are  7.48  and  8.30  cracks/cm3-ns  and 
-14.1  and  -13.4  GPa"1.  respectively.  The  experimental  and  computed  intercept  and  slope  values 
for  the  growth  curves  in  Figure  J-8(b)  are  59.5  and  59.4  jam  and -0.161  and -0.164  GPa- Vs-1, 
respectively.  For  the  experimental  growth  curve,  the  points  for  Shots  110,  91,  and  113  were  out 
of  range  of  the  majority  of  the  points  and  were  therefore  not  used  in  the  least-squares  fit.  The 
dynamic  fracture  parameters  are  given  in  Table  J-3.  The  exclusion  of  the  stringer  cracks  appreci¬ 
ably  affected  only  the  nucleation  threshold  rate  N().  This  parameter  decreased  from  90  to  42 
cracks/cm3 -ns  when  exluding  the  stringer  cracks.  Figures  J-9  through  J-l  3  give  the  input  listings 
for  the  PUFF  8  computations.  Figures  J-l 4  through  J-22  give  computed  shock  histories  for  the  cell 
of  maximum  damage  for  each  of  the  heat  treatment  A  specimens  with  stringer  cracks  excluded. 
Plots  (a)  and  (b)  give  the  principal  stresses  in  the  direction  of  shock  wave  propagation  and 
orthogonal  to  that  direction,  respectively.  Plots  (c)  and  (d)  give  particle  velocity  and  crack 
concentration  histories,  respectively.) 
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Figure  J-1.  (Continued) 
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Figure  J-2.  (Continued) 
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Figure  J-4.  (Continued) 
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Figure  J-4.  (Continued) 
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Figure  J-7.  (Continued) 
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Figure  J-8 
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Experimental  and  computed  nucleation  rate  and  growth  curves  for  HF-1  steel  with 
heat  treatment  A  with  stringer  cracks  excluded,  (a)  Nucleation  rate  curves, 
(b)  Growth  curves.  The  points  in  these  plots  are  from  Table  J-2  and  were  obtained 
by  using  the  final  set  of  dynamic  fracture  parameters  for  this  heat  treatment  of 
HF-1  steel  in  a  PUFF  8  computation. 
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Figure  J-10.  PUFF  input  parameters  for  {a)  Shot  92  and  (b)  Shot  91  with  stringer  cracks  excluded. 


Figure  J-l  1.  PUFF  input  parameters  for  (a)  Shot  90 
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Figura  J-13.  PUFF  input  parameters  for  Shot  1 12  with  stringer  cracks  excluded. 


Computed  shock  histories  in  cell  23  for  Shot  110  with  stringer  cracks  excluded,  (a)  Parallel  stress, 
(b)  Perpendicular  stress,  (c)  Particle  velocity,  (c)  Crack  concentration. 


SHOT  98 
CEU.  23 


i  cell  23  for  Shot  98  with  stringer  cracks  excluded,  (a)  Parallel  stress. 
Particle  velocity,  (c)  Crack  concentration. 
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Computed  shock  histories  in  cell  23  for  Shot  91  with  stringer  cracks  excluded,  (a)  Parallel  stress, 
(b)  Perpendicular  stress,  (c)  Particle  velocity,  (c)  Crack  concentration. 
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Figure  J-18.  Computed  shock  histories  in  cell  23  for  Shot  90  with  stringer  cracks  excluded,  (a)  Parallel  stress, 
(b)  Perpendicular  stress,  (c)  Particle  velocity,  (c)  Crack  concentration. 


Figure  J-19.  Computed  shock  histories  in  cell  23  for  Shot  94  with  stringer  cracks  excluded,  (a)  Parallel  stress 
(b)  Perpendicular  stress,  (c)  Particle  velocity.  <c)  Crack  concentration. 


Computed  shock  histories  in  cell  20  for  Shot  113  with  stringer  cracks 
(b)  Perpendicular  stress,  (c)  Particle  velocity,  (c)  Crack  concentration. 


Figure  J-21 .  Computed  shock  histories  in  cell  24  for  Shot  1 1 1  with  stringer  cracks  excluded,  (a)  Parallel  stress, 

(b)  Perpendicular  stress,  (c)  Particle  velocity,  (c)  Crack  concentration. 


Computed  shock  histories  in  ceil  24  for  Shot  112  with  stringer  cracks  excluded,  (a)  Parallel  stress, 
(b)  Perpendicular  stress,  (c)  Particle  velocity,  (c)  Crack  concentration. 


Table  J-1.  Summary  of  experimental  volume  crack-size  distribution  parameters  for  spall  fracture 
specimens  of  HF-1  steel  with  heat  treatment  A  with  stringer  cracks  excluded.* 


Distance  from 

Number  of 

Average 

Specimen  Impact 

Cracks  per 

Crack 

Surface  to  Midpoint 

Unit  Volume 

Radius 

of  Zone  of 

N0 

Ri 

Maximum  Damage 

5b 

0.03 

0 

0.74 

8 

3.28 

2 

2.39 

The  parameters  Nq  and  R]  were  obtained  for  each  damage  zone  of  each  specimen  by  fitting  the  equation  Ng  =  Ngexp  (-  R/Rj)  to 
the  volume  crack-size  distribution.  Here  Ng  is  the  number  of  cracks  per  unit  volume  with  radius  greater  than  R.  Plots  of  the  volume 
crack-size  distributions  for  each  damage  zone  of  each  specimen  containing  no  stringer  cracks  are  given  in  Appendix  E  (Shots  95, 
110,  88,  and  111)  and  this  appendix  (Shots  98,  92,  91,  90,  94,  1 13,  and  1 12).  The  parameters  Nq  and  Rj  listed  in  this  table 
are  the  values  for  the  zone  of  maximum  damage. 

This  shot  was  not  used  for  determining  dynamic  fracture  parameters  due  to  the  small  number  of  measured  cracks. 

Cracks  were  not  measured  for  the  no-damage  Shot  89. 

This  full-spall  shot  was  not  used  for  determining  dynamic  fracture  parameters.  Only  cracks  in  one  half  of  the  specimen  were 
digitized. 
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Tabia  J-2.  Summary  of  axparimantal  and  computad  nudaation  and  growth  valuai  for  spall  fractura  spacimara  of  KF-1  staal  with  haat  traatmant  A  with  stringer  cracks  excluded.' 


Table  J-3.  Dynamic  fracture  parameters  for  heat  treatment  A  of  HF-1  steel  with  stringer  cracks 
excluded. 


Parameter3 

Description 

Units 

Value 

T, 

Growth  coefficient 

GPa" 1  ps" 1 

-0.409 

V 

Growth  threshold  stress 

GPa 

-0.1 

Ro 

Nucleation  size  parameter 

pm 

58.7 

N0 

Nucleation  threshold  rate 

no.cm-3  ns-1 

42.2 

ano 

Nucleation  threshold  stress 

GPa 

-2.5 

Nucleation  sensitivity  parameter 


GPa 


-0.0577 


APPENDIX  K 

FRAMING  CAMERA  PHOTOGRAPHS  FOR  EXPLOSIVE-FILLED 
ARMCO  IRON  CYLINDER  EXPERIMENT 

(This  appendix  contains  a  scries  of  24  framing  camera  photographs  for  the  explosive-filled  Armco 
iron  cylinder  experiment.  The  frames  are  numbered  1  through  25;  frame  13  is  micsing.  The 
composition  B  explosive  detonation  front  was  located  inside  the  cylinder  in  frames  1  through 
10.  The  detonation  front  was  assumed  to  be  located  at  the  left  end  of  the  cy'inder  in  frame  11. 
The  calculated  axial  position  of  the  detonation  front  for  these  frames  is  given  in  Table  12.  In 
these  figures  the  time  corresponding  to  each  frame  is  taken  as  the  time  after  the  explosive  detonat¬ 
ion  front  has  entered  the  cylinder.  These  times  are  probably  uncertain  by  a  few  microseconds  since 
the  detonation  front  was  probably  not  located  exactly  at  the  end  of  the  cylinder  in  frame  11. 
The  initial  explosive  detonation  occurred  about  8  ps  prior  to  the  entrance  of  the  detonation  front 
into  the  cylinder.  The  explosive  extended  about  64  mm  beyond  the  metal  cylinder  on  the  initiation 
end.) 


(b) 


Figure  K-2.  Frames  3  and  4  for  Armco  iron  framing  camera  experiment,  (a) 
(b)  frame  4  at  18  jus. 
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Figure  K-8.  Frames  16  and  17  for  Armco  iron  framing  camera  experiment,  (a)  Frame  16  at  68  ps 
and  (b)  frame  17  at  73  ps. 
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APPENDIX  L 

PHOTOGRAPHS  OF  FRAGMENT  TYPES  FROM  ARMCO  IRON  AND 
HF-1  STEEL  CYLINDER  FRAGMENT  RECOVERY  EXPERIMENTS 

(This  appendix  retains  24  figures.  Figures  L-l  through  L-8  show  Armco  iron  fragments  from 
cylinder  experiment  2.  No  type  4C  fragments  were  observed  for  the  Armco  iron  experiments. 
Figures  L-9  through  L-l 6  show  HF-1  steel  fragments  with  heat  treatment  A  from  cylinder  ex¬ 
periment  4.  Figures  L-l 7  through  L-24  show  HF-I  steel  fragments  with  heat  treatment  B  from 
cylinder  experiment  5.  The  fragment  types  are  defined  in  Table  14.  Schematic  examples  of  the 
fragment  types  are  shown  in  Figure  35). 


L-l 


surface,  (b)  cylinder  inner  surface,  and 


Figure  L-2.  Type  1  Armco  iron  cylinder  fragment  showing  (a)  cylinder  outer  surface,  (b)  cylinder  inner  surface  and 
(c)  fragment  end  surface.  Fragment  mass:  437.4  gr  (28.34  g). 


i  cylinder  fragment  showing  (a)  cylinder  outer  surface,  (b)  brittle  and  shear  fracture  surfaces, 
surface.  Fragment  mass:  276.9  gr  (17.94  g). 


pe  2B  Armco  iron  cylinder  fragment  showing  (a)  cylinder  outer  surface,  (b)  brittle  fracture  surfaces  and 
fragment  end  surface.  Fragment  mass:  104.8  gr  (6.75  g). 


Type  3A  Armco  iron  cylinder  fragment  showing  (a)  brittle  and  shear  fracture  surfaces,  (b)  cylinder  inner  surface, 
and  (c)  fragment  end  surface.  Fragment  mass:  161.5  gr  (10.47  g). 


pe  3B  Armco  iron  cylinder  fragment  showing  (a)  shear  fracture  surfaces,  (b)  cylinder  inner  surface,  and 
fragment  end  surface.  Fragment  mess:  188.9  gr  (12.24  g). 


pe  4A  Armco  iron  cylinder  fragment  showing  (a)  brittle  fracture  surfaces,  (b)  shear  fracture  surfaces,  and 
fragment  end  surface.  Fragment  mass:  25.1  gr  (1.63  g). 


Type  4B  Armco  iron  cylinder  fragment  showing  (a)  brittle  fracture  surfaces  and  (b)  fragment  end  surface 
Fragment  mass:  14.9  gr  (0.97  g). 


Figure  L-9.  Type  1  HF-1  steel  cylinder  fragment  with  heat  treatment  A  showing  (a)  cylinder  outer  surface,  (b)  cylinder  inner 
surface,  and  (c)  fragment  end  surface.  Fragment  mass:  347.3  gr  (22.51  g). 


2A  HF-1  steel  cylinder  fragment  with  heat  treatment  A  showing  (a)  cylinder  outer  surface,  (b)  brittle  and 
fracture  surfaces,  and  (c)  fragment  end  surface.  Fragment  mass:  373.2  gr  (24.18  g). 


Figure  L-11.  Type  ?.B  HF-1  steel  cylinder  fragment  with  heat  treatment  A  showing  (a)  cylinder  outer  surface,  (b)  brittle 
fracture  surfaces,  and  (c)  fragment  end  surface.  Fragment  mass:  165.9  gr  (10.75  g). 


Msxtri> 


tee!  cylinder  fragment  with  heat  treatment  A  showing  (a)  brittle  and  shear  fracture  surfaces, 
surface,  and  (c)  fragment  end  surface.  Fragment  mass:  275.7  gr  (17.87  g). 


Type  3B  HF-1  steel  cylinder  fragment  with  heat  treatment  A  showing  (a)  shear  fracture  surfaces,  (b)  cylinder 
inner  surface,  and  (c)  fragment  end  surface.  Fragment  mass:  77.1  gr  (4.99  g). 


Type  4A  HF-1  steel  cylinder  fragment  with  heat  treatment  A  showing  (a)  brittle  fracture  surfaces,  (b)  shear 
fracture  surfaces,  and  (c)  fragment  end  surface.  Fragment  mass:  112.7  gr  (7.30  g). 


ae  4B  HF-1  steel  cylinder  fragment  with  heat  treatment  A  showing  (a)  brittle  fracture  surfaces  and 
fragment  end  surface.  Fragment  mass:  42.8  gr  (2.78  g). 


pe  4C  HF-1  steel  cylinder  fragment  with  heat  treatment  A  showing  (a)  shear  fracture  surfaces  and  (b)  fragment 
d  surface.  Fragment  mass:  17.5  gr  {1.13  g). 


r 


treatment  B  showing  (a)  cylinder  outer  surface,  (b)  cylinder  inner 


2A  HF-1  steel  cylinder  fragment  with  heat  treatment  B  showing  (a)  cylinder  outer  surface,  (b)  brittle  and 
fracture  surfaces,  and  (c)  fragment  end  surface.  Fragment  mass:  231.9  gr  (15.03  g). 


Type  2B  HF-1  steel  cylinder  fragment  with  heat  treatment  B  showing  (a)  cylinder  outer  surface,  lb)  brittle 
fracture  surfaces,  and  (c)  fragment  end  surface.  Fragment  mass:  116.5  gr  (7.55  g). 


Type  3B  HF-1  steel  cylinder  fragment  with  heat  treatment  B  showing  {a)  shear  fracture  surfaces,  (b)  cylinder 
inner  surface,  and  (c)  fragment  end  surface.  Fragment  mass:  72.5  gr  (4.70  g). 


APPENDIX  M 

RESULTS  OF  THE  ARMCO  IRON  CYLINDER,  THE  HF-1  STEEL  CYLINDER, 
AND  THE  HF-1  STEEL  PROJECTILE  FRAGMENT  RECOVERY  EXPERIMENTS 


(This  appendix  contains  tables  of  the  results  of  the  Armco  iron  cylinder,  the  HF-1  steel  cylinder, 
and  the  HF-1  steel  projectile  fragment  recovery  experiments.  The  configuration  details  for  these 
experiments  are  given  in  Table  13.  The  results  are  summarized  in  Section  VIII.  Tables  M-l  through 
M-7  give  the  mass  distribution  results  for  the  typed  fragments  for  the  experiments.  The  fragments 
in  the  0-1  gr  mass  group  of  each  table  were  not  counted  due  to  their  large  numbers  and  small  sizes. 
The  recovered  fragment  mass  was  at  least  98.9%  of  the  initial  metal  mass  in  all  the  experiments. 
Tables  IW-8  through  M-l  4  give  dimension  measurements  and  strain  calculations  for  selected  type  1 
fragments  from  each  experiment.  All  type  1  fragments  from  the  HF-1  steel  cylinder  experiments 
were  measured.  Both  English  and  SI  units  are  used  in  these  tables.  (1  grain  (gr)  =  1/7000  pound 
=  0.06480  gram  (g)). 


M-l 


Table  M-1.  Mass  distribution  and  characterization  of  Armco  iron  cylinder  fragments  from 
experiment  1.* 


Mass 

Group 

(gr) 

Type  1 

Total 

No. 

Mass 

No. 

Mass 

(gr) 

(g) 

(gr) 

(g) 

0-lb 

_ 

1326.7 

85.97 

1-5 

637 

1423.8 

92.26 

5-10 

146 

1055.4 

68.39 

10-15 

84 

1037.8 

67.25 

15-20 

45 

778.4 

50.44 

20-25 

37 

829.5 

53.75 

25-30 

1 

27.3 

1.77 

38 

1047.8 

67.90 

30-40 

2 

70.8 

4.58 

4i 

1403.5 

90.95 

40-50 

6 

278.6 

18.05 

37 

1673.0 

108.41 

50-75 

18c 

1159.4 

75.13 

73 

4514.0 

292.51 

75-100 

13d 

1196.8 

77.55 

40 

3594.8 

232.94 

100-125 

12c 

1351.1 

87.55 

27 

3012.2 

195.19 

125-150 

16f 

2162.3 

140.12 

33 

4484.3 

290.58 

150-200 

10» 

1788.6 

115.90 

17 

2966.0 

192.20 

200-300 

24h 

5911.6 

383.07 

32 

7927.1 

513.68 

300-600 

47* 

21369.3 

1384.73 

48 

21725.0 

1407.78 

600-1000 

21s 

19884.0 

1288.50 

28 

20539.8 

1330.98 

1000-1500 

13* 

16982.1 

1100.44 

15 

19658.0 

1273.86 

1500+ 

17* 

41128.9 

2665.15 

17 

41128.9 

2665.15 

a  For  this  experiment  only  type  1  fragments  were  characterized.  Six  fragments  (one  fragment  from  the  125  •  150  gr  group,  two  frag- 
from  the  200  -  300  gr  group,  one  fragment  from  the  600  -  1000  gr  group,  and  two  fragments  from  the  1000  -  1500  gr  group) 
with  a  total  mass  of  0.216  kg  were  previously  removed  from  the  distribution  for  rr.etallographic  analysis  and  therefore  were  not 
characterized.  Seventy-five  of  the  206  total  type  1  fragments  came  from  the  end  of  the  cylinder. 
bThe  fragments  were  not  counted  in  this  mass  group. 
c  One  fragment  came  from  the  end  of  the  cylinder. 

“  Five  fragments  came  from  the  end  of  the  cylinder. 
e  Four  fragments  came  from  the  end  of  the  cylinder. 

-  Two  fragments  came  from  the  end  of  the  cylinder. 
g  Eight  fragments  came  froih  the  end  of  the  cylinder. 
b  Nine  fragments  came  from  the  end  of  the  cylinder. 

J  Twenty-six  fragments  came  from  end  of  the  cylinder. 

J  Six  fragments  came  from  the  end  of  the  cylinder. 


Table  M-2.  Mass  distribution  and  characterization  of  Armco  iron  cylinder  fragments  from  experiment  2. 


<95*500  |e  <99.9  *  *99  9  32-39 

520*525  I  520  2  33.71  I  *202  JW* 

525*530  2*  1052.8  68.22  2  1052  8  68  22 


ass  distribution  and  characterization  of  Armco  iron  cylinder  fragments  from  experiment  2.a  (Continued) 


Table  M-2.  Mass  distribution  and  characterization  of  Armco  iron  cylinder  fragments  from  experiment  2.a  (Continued) 


Table  M-3.  Mass  distribution  and  characterization  of  HF-1  steel  cylinder 
fragments  with  heat  treatment  A  from  experiment  3,a 


Mass 

Group 

(gr) 

0-lb 

1-5 

5-10 

10-15 

15-20 

20-25 

25-30 

30-40 

40-50 

50-75 

75-100 

100-125 

125-150 

150-200 

200-300 

300-600 


Type  1 


(gr) " 

(g) 

No. 

(gr) 

8568.6 

5521 

12758.9 

1482 

10564.0 

755 

9264.6 

464 

8045.2 

307 

6888.8 

232 

6322.8 

351 

12089.4 

41.5 

2.69 

234 

10502.6 

193.5 

12.54 

370 

22785.1 

257.6 

16.69 

183 

15836.0 

476.2 

30.86 

121 

13562.1 

47 

6451.4 

1261.8 

81.77 

58 

9973.0 

3158.4 

204.66 

47 

11255.3 

1096.6 

71.06 

9 

3141.6 

555.24 
826.78 
684.55 
600.35 

521.33 

446.39 
409.72 

783.39 

680.57 
1476.47 
1026.17 
878.82 
418.05 

646.25 

729.34 

203.58 


I'or  this  experiment  only  type  1  fragments  were  characterized.  Six  fragments  (one 
fragment  from  the  160-200  gr  group,  two  frapnents  from  the  200-300  gr  group, 
and  three  fragments  from  the  300-600  gr  group)  with  a  total  mass  of  0.1 17  kg  were 
previously  removed  from  the  distribution  for  metaJIographic  analysis  and  therefore 
were  not  characterized.  Twelve  of  the  33  total  type  1  fragments  came  from  the  end 
of  the  cylinder. 

b  The  fragments  were  not  counted  in  this  mass  group.  Only  type  I  fragments  were 
characterized. 

c  One  fragment  came  from  the  end  of  the  cylinder. 
dTwo  fragments  came  from  the  end  of  the  cylinder. 
c  Eight  fragments  came  from  the  end  of  the  cylinder. 


Table  M-4.  Mass  distribution  and  characterization  of  HF-1  steel  cylinder  fragments  with  heat  treatment  A  from  experiment  4.a  (Continued) 
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Table  M-5.  Mass  distribution  and  characterization  of  HF-1  steel  cylinder  fragments  with  heat  treatment  B  from  experiments  5.a  (Continued) 


Table  M-6.  Mass  distribution  and  characterization  of  HF-1  steel  cylinder  fragments  with  heat 
treatment  B  from  experiment  6.a 


Mass 

Group 

(gr) 

Type  1 

Total 

No. 

Mass 

No. 

Mass 

(gr)  (g) 

(gr) 

(g) 

0-lb 

11926.0 

772.80 

1-2.5 

4680 

7458.3 

483.30 

2.5-5 

2616 

9316.4 

603.70 

5-7.5 

1268 

7858.6 

509.24 

7.5-10 

735 

6375.4 

413.13 

10-12.5 

546 

6104.9 

395.60 

12.5-15 

421 

5757.0 

373.05 

15-17.5 

327 

5307.9 

343.95 

17.5-20 

259 

4833.7 

313.22 

20-22.5 

241 

5102.1 

330.62 

22.5-25 

193 

4593.6 

297.67 

25-27.5 

190 

5008.4 

324.54 

27.5-30 

128 

3667.4 

237.65 

30-32.5 

131 

4096.1 

265.43 

32.5-35 

101 

3422.5 

221.78 

35-37.5 

94 

3411.1 

221.04 

37.5-40 

110 

4252.5 

275.56 

40-42.5 

75 

2082.9 

199.77 

42.5-45 

63 

2756.5 

178.62 

45-47.5 

56 

2585.7 

167.55 

47.5-50 

62 

3025.6 

196.06 

50-52.5 

46 

2362.1 

153.06 

52.5-55 

46 

2484.1 

160.97 

55-57.5 

lc 

55.1  3.57 

46 

2587.2 

167.65 

57.5-60 

41 

2415.8 

156.54 

60-62.5 

1 

61.1  3.96 

36 

2206.2 

142.96 

62.5-65 

28 

1776.6 

115.12 

65-67.5 

36 

2388.9 

154.80 

67.5-70 

26 

1787.2 

115.81 

70-72.5 

21 

1492.2 

97.02 

72.5-75 

14 

1032.9 

66.93 

75-77.5 

18 

1376.8 

89.22 

77.5-80 

1° 

77.9  5.05 

22 

1731.6 

112.21 

80-82.5 

22 

1792.7 

116.17 

82.5-85 

12 

1005.7 

65.17 

85-87.5 

6 

519.0 

33.63 
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Table  M-6.  Mass  distribution  and  characterization  of  HF-1  steel  cylinder  fragments  with  heat 
treatment  B  from  experiment  6.1  (Continued) 


87.5-90 

13 

1152.0 

74.65 

90-92.5 

12 

1095.6 

70.99 

92.5-95 

9 

837.7 

54.28 

95-97.5 

9 

865.4 

56.08 

97.5-100 

8 

789.4 

51.15 

100-105 

1 

101.4 

6.57 

16 

1634.7 

105.03 

105-110 

1 

107.3 

6.95 

17 

1824.1 

118.20 

110-115 

3 

334.0 

21.64 

16 

1800.3 

116.66 

115-120 

1 

115.2 

7.47 

13 

1521.8 

98.61 

120-125 

1 

122.9 

7.96 

11 

1351.9 

87.60 

125-130 

4C 

508.5 

32.95 

12 

1530.4 

99.17 

130-135 

2 

264.4 

17.13 

6 

793.0 

51.39 

135-140 

4d 

549.1 

35.58 

8 

1097.1 

71.09 

140-145 

2 

283.4 

18.36 

8 

1137.3 

73.70 

145-150 

3e 

442.7 

28.69 

6 

882.7 

57.20 

150-155 

2C 

305.6 

19.80 

5 

763.1 

49.45 

155-160 

lc 

157.4 

10.20 

4 

631.8 

40.94 

160-165 

1 

164.1 

10.63 

3 

486.2 

31.51 

165-170 

3C 

503.1 

32.60 

5 

841.2 

54.51 

170-175 

1 

174.2 

11.29 

2 

347.2 

22.50 

175-180 

lc 

178.0 

11.53 

2 

354.2 

22.95 

180-185 

3C 

546.5 

34.51 

4 

728.4 

47.20 

185-190 

1 

185.4 

12.01 

3 

563.3 

36.50 

190-195 

2e 

384.3 

24.90 

3 

574.3 

37.21 

195-200 

1 

199.3 

12.91 

205-210 

2e 

415.7 

26.94 

3 

625.7 

40.55 

210-215 

1 

210.7 

13.65 

220-225 

lc 

220.1 

14.26 

2 

440.9 

28.57 

230-235 

1° 

234.7 

15.21 

1 

234.7 

15.21 

235-240 

2 

471.0 

30.52 

240-245 

1 

242.2 

15.69 

245-250 

2° 

496.0 

32.14 

2 

496.0 

32.14 

290-295 

lc 

291.8 

18.91 

1 

293.8 

19.04 

315-320 

1 

319.2 

20.68 

1 

319.2 

20.68 

325-330 

1 

,  328.1 

21.26 

1 

328.1 

21.26 

a  If  a  mass  group  is  not  listed  then  no  fragments  were  observed  in  that  group.  For  this  experiment  only  type  1  fragments  were 
characterized.  Twenty-three  of  the  49  total  type  1  fragments  came  from  the  end  of  the  cylinder. 
b  The  fragments  were  not  counted  in  this  mass  group. 

®  One  fragment  came  from  the  end  of  the  cylinder. 

“  Four  fragments  came  front  the  end  of  the  cylinder. 
c  Two  fragments  came  from  the  end  of  the  cylinder. 


Table  M-7.  Mass  distribution  and  characterization  of  HF-1  steel  projectile  fragments  with  heat  treatment  A  from  experiment  7.a 
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Table  M-7.  Mass  distribution  and  characterization  of  HF-1  steel  projectile  fragments  with  heat  treatment  A  from  experiment  7.a  (Continued) 


Table  M-7.  Mass  distribution  and  characterization  of  HF-1  steel  projectile  fragments  with  heat  treatment  A  from  experiment  7.a  (Continued) 
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us  nuss  group. 


Table  M-8.  Measurements  on  twenty-five  selected  type  1  Armco  iron  cylinder  fragments  from  experiment  1 


Fragment 

No. 

(p) 

Fragment 

Maas 

<l) 

Wall 

Thkkne«*> 

(mm) 

Average 

Fraction  of  Wail 
Thkkneaa  with 
Brittle  Fracture® 

.Maximum 

Fragment 

Width** 

(mm) 

Maximum 

Fragment 

Lengthc 

(mm) 

Engineering 

Radial 

Strain* 

True 

Radial 

Strata* 

Homogeneous 
Pi. is  tic 
Strain11 

' 

69.1 

4.48 

13.1 

0.60 

7 

17 

0.31 

0.38 

0.65 

2 

93.0 

6.02 

12.6 

0.59 

6 

11 

22 

0.34 

0.41 

0.72 

3 

120.2 

7.79 

11.7 

0.66 

5 

17 

25 

0.39 

0.49 

0.85 

4 

142.3 

9.24 

11.4 

0.65 

5 

15 

21 

0.40 

0.52 

0.89 

5* 

181.5 

11.76 

12.1 

C.55 

5 

14 

23 

0.36 

0  45 

0.78 

6 

254.9 

16.52 

12.4 

0.56 

8 

11 

37 

0.35 

0.43 

0.74 

7* 

267.4 

17.33 

12.8 

0.63 

6 

11 

23 

0.33 

0.40 

0.69 

8 

282.5 

18.31 

12.3 

0.66 

9 

15 

34 

0.36 

0.44 

0.76 

9 

375.6 

24.34 

12.3 

0.46 

5 

19 

38 

0.36 

0.44 

0.76 

10 

399.7 

25.90 

12.7 

0.67 

6 

12 

47 

0.33 

0.41 

0.70 

11 

476.8 

30.89 

13.2 

0.61 

6 

20 

51 

0.31 

0.36 

0.63 

12* 

543.1 

35.19 

13.0 

0.45 

8 

22 

37 

0.32 

0.39 

0.67 

13 

613.1 

39.73 

12.7 

0.53 

9 

20 

41 

0.33 

0.41 

0.71 

14* 

642.6 

41.64 

IM 

0.44 

12 

17 

51 

0.40 

0.51 

0.88 

13 

686.1 

44,46 

11.7 

0.50 

11 

18 

51 

0.39 

0.49 

0.85 

16* 

863.8 

55.97 

12.3 

0.49 

7 

18 

69 

0.35 

0.44 

0.76 

17 

1127.1 

73.04 

13.6 

0.51 

8 

19 

89 

0.29 

0.34 

0.59 

18* 

1175.5 

76.18 

12.8 

0.61 

n 

15 

102 

0.33 

0.40 

0.69 

19 

1271.3 

82.38 

13.0 

0.59 

6 

17 

111 

0.32 

0.39 

0.67 

20* 

1419.9 

92.01 

12.4 

0.47 

9 

22 

114 

0.35 

0.43 

0.74 

2!1 

1700.4 

110.18 

13.0 

0.47 

7 

14 

130 

0.32 

0.39 

0.67 

22* 

2003.6 

129.83 

12.8 

0.55 

8 

22 

114 

0.33 

0.40 

0.70 

23 

2241.7 

145.26 

13.9 

0.47 

14 

20 

127 

0.27 

0.31 

0.54 

24 

3041,5 

197.09 

13.4 

0.58 

18 

14 

143 

0.30 

0.35 

0.54 

25 

3118.5 

202.08 

12.7 

0.55 

20 

23 

121 

0.33 

0.41 

0.71 

*  The  25  fragments  were  chosen  to  represent  the  entire  mass  range  for  type  1  fragments. 

^The  initial  wall  thickness  was  19.1  mm. 

c  These  values  are  the  average  of  measurements  from  the  outer  cylindrical  surface  to  the  beginning  of  the  shear  fracture  surface  on  eath  side  of  the  fragment  The  estimated  un¬ 
certainly  in  these  values  is  J0-20%, 

^  The  two  values  listed  for  a  frament  are,  respectively,  the  maximum  width  of  the  brittle  fracture  region  near  the  outer  cylindrical  surface  and  the  rnaaimum  width  of  the  shear 
fracture  region  near  the  Inner  cylindrical  surface. 

0  Tho  fragment  tength  is  parallel  to  the  cylindrical  axis. 

*  The  engineering  radial  strain  is  (a0  •  af )/a0  where  a0  and  af  are  the  initial  and  final  wall  thicknesses,  respectively. 

8  The  true  radial  strain^  1$  •  In(af/a0). 

n  The  homogeneous  plastic  strain^  is  -  v  3  ln(a//a0). 

*  This  fragment  came  from  the  end  of  the  cylinder. 
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Table  M-9.  Measurements  on  twenty-five  selected  type  1  Armco  iron  cylinder  fragments  from  experiment  2.a 


Fragment 

No. 

Fragment 

Maaa 

( P ) 

(8) 

Wall 

Thkkne»b 

(mm) 

Average 

Fraction  ol  Wall 

Thkkneaa  with 
Brittle  Fracture® 

Maximum 

Fragment 

Width0 

(mra) 

Maximum 

Fragment 

Length® 

(mm) 

Engineering 

Radial 

Strain® 

True 

Radial 

Strain® 

Homogeneous 

Plastic 

Strain0 

1 

32.1 

2.08 

9.7 

0.82 

6 

10 

0.49 

0.68 

1.18 

2 

43.9 

2.85 

8.2 

0.69 

7 

12 

11 

0.57 

0.84 

1.45 

3 

63.8 

4.13 

9.9 

0.62 

7 

13 

14 

0.48 

0.66 

1.14 

4 

82.1 

5.32 

10.1 

0.63 

13 

11 

13 

0.47 

0.64 

1.11 

5 

107.1 

6.94 

12.5 

0.46 

7 

12 

18 

0.35 

0.43 

0.74 

6 

130.4 

9.75 

13.3 

0.54 

6 

14 

22 

0.30 

0.36 

0.63 

7 

206.3 

13.37 

12.1 

0.56 

9 

13 

27 

0.36 

0.45 

0.78 

8d 

282.3 

18.29 

12.2 

0.51 

10 

18 

31 

0.36 

0.44 

0.77 

9 

334.2 

21.66 

12.2 

0.56 

9 

16 

43 

0.36 

0.45 

0.77 

10 

412.3 

26.73 

13.4 

0.57 

3 

15 

38 

0.30 

0.35 

0.61 

11 

481.3 

31.20 

12.6 

0.56 

6 

15 

57 

0.34 

0.41 

0.71 

12 

677.2 

43.88 

12.5 

0.50 

14 

23 

57 

0.34 

0.42 

0.73 

13 

ms 

50.43 

13.3 

0.57 

11 

21 

51 

0.30 

0.36 

0.62 

14 

904 S 

58.61 

13.3 

0.63 

5 

16 

83 

0.30 

0J6 

0.62 

13 

2074.0 

134.39 

12.8 

0.55 

11 

21 

83 

0.33 

0.40 

0.69 

I6d 

1322.8 

85.72 

11.5 

0.53 

11 

21 

79 

0.40 

0.51 

0.88 

17 

1402.3 

90.87 

13.1 

0.49 

h 

17 

102 

0.31 

0.38 

0.65 

18 

1569.5 

101.71 

13.3 

0.59 

13 

21 

108 

0.30 

0.36 

0.63 

19 

1776.1 

115.09 

13.1 

0.65 

13 

22 

89 

0.31 

0.38 

0.65 

20d 

1912.9 

123.95 

13.1 

0.60 

16 

19 

95 

0.31 

0.37 

0.65 

21 

2167.3 

140.44 

12.7 

0.60 

14 

25 

102 

0.33 

0.41 

0.70 

22d 

2402.8 

155.70 

13.3 

0.61 

8 

14 

159 

0.30 

0.36 

0.63 

33d 

2622,1 

169.91 

12.7 

0.46 

12 

12 

203 

0.33 

0.41 

0.70 

24 

2774.6 

179.79 

14.9 

0.53 

14 

25 

127 

0.22 

0.24 

0.42 

23d 

3333.0 

215.98 

11.0 

0,53 

24 

24 

184 

0.42 

0.54 

0.94 

*  The  25  fraftnents  were  chosen  i  __  represent  the  entire  mass  range  for  type  1  fragments. 
b,nie  Initial  wall  thickness  was  19.1  mm. 

®  Defined  in  footnotes  In  Table  M*8. 

®'fliis  fragment  came  from  the  end  of  the  cylinder. 
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Table  M-10.  Measurements  on  the  type  1  HF-1  steel  cylinder  fragments  with  heat  treatment  A 
from  experiment  3. 


Tiapnent 

No. 

<**> 

Fragment 

Mas 

(1) 

Will 

Ifckknea* 

(mm) 

Average 

Fraction  of  Mail 

lkkkness  with 
Brittle  Fracturtb 

Maximum 

Fragment 

Width*' 

(mm) 

Maximum 

Fragment 

Length*3 

(mm) 

Engineering 

Radial 

Strain** 

True 

Radial 

St»ainb 

Homoft  neons 
Haile 
Su.ln*> 

\c 

415 

2  69 

18  6 

084 

3  0 

7.1 

15  2 

016 

018 

031 

2 

62  8 

4  07 

195 

077 

5  8 

5  7 

116 

0.12 

013 

0  23 

3 

64  9 

4  20 

19  3 

0.75 

29 

45 

106 

0.13 

014 

0  25 

4 

65  8 

4  26 

19  9 

088 

50 

6  1 

12  6 

on 

0.11 

019 

5 

78  0 

503 

19.4 

0  73 

40 

74 

189 

013 

0.14 

0  24 

6 

83  3 

540 

197 

0  83 

59 

67 

16  5 

012 

012 

021 

7 

96,4 

6  25 

201 

061 

39 

7.6 

16.1 

0.10 

010 

018 

8 

101.2 

6.56 

19  5 

0  87 

58 

47 

201 

012 

0.13 

0  23 

9 

104  8 

6  79 

20  2 

0  80 

7.6 

11.0 

20  2 

0.09 

0.10 

0.17 

10 

115  9 

741 

201 

0  73 

33 

102 

289 

0.10 

0.10 

018 

II 

151  1 

9  79 

201 

0  83 

55 

107 

28.1 

0.10 

0.10 

018 

I2C 

171.4 

tl.l 

19.3 

0  87 

11.9 

74 

184 

0.13 

0.14 

0.25 

13 

179  8 

116 

204 

0.76 

80 

98 

24.8 

008 

009 

015 

14 

183  0 

119 

19.7 

0  78 

84 

130 

19.6 

0.11 

012 

021 

15' 

185.4 

12.0 

19  0 

063 

96 

68 

23  0 

015 

0.16 

027 

16 

191,4 

12.4 

19.7 

091 

59 

7.9 

25  0 

0.11 

012 

021 

17 

199  9 

130 

19  9 

081 

69 

92 

210 

011 

on 

019 

18" 

2004 

13.0 

19  2 

068 

83 

14  7 

17  6 

014 

015 

026 

I9C 

2080 

134 

19  3 

0.78 

8.9 

84 

194 

013 

0.14 

025 

20 

219  2 

14  2 

201 

0.77 

90 

82 

27.6 

0.10 

010 

018 

21 

220  7 

14  7 

19.7 

0  71 

99 

136 

17  8 

Oil 

0.12 

021 

22' 

2239 

144 

19  2 

0.78 

106 

99 

17.7 

014 

0.15 

026 

23 

226.0 

14  6 

19  6 

065 

98 

11.7 

214 

0.12 

013 

0.22 

24' 

227  8 

14.8 

19  8 

071 

136 

14  3 

227 

Oil 

0.12 

021 

2J' 

248.4 

16.1 

19  6 

084 

94 

10  6 

201 

012 

0.13 

022 

24' 

259  2 

16  8 

19  2 

0  81 

166 

124 

234 

014 

0.15 

0  25 

27 

263  7 

17.1 

19  9 

0  80 

74 

99 

261 

on 

0.1 1 

019 

28' 

277,5 

180 

19  6 

0  85 

137 

105 

26  2 

012 

0.13 

0,22 

29c 

2860 

184 

19  7 

084 

13  5 

16  8 

26.7 

0.12 

0.12 

0  21 

30 

*91  6 

189 

19  7 

078 

104 

11.7 

269 

0.12 

012 

021 

31' 

3*1  8 

20  2 

19  7 

075 

II.I 

13.4 

27,6 

012 

012 

021 

32 

?24  6 

210 

19,8 

060 

8,7 

122 

313 

Oil 

012 

020 

33 

460,2 

29.8 

20,3 

042 

89 

374 

0.09 

009 

0.16 

m 
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*  The  initial  wall  thickness  wu  22  3  mm 
b  Defined  in  footnote*  in  Table  M  8. 
c  Tris  frapnent  time  from  the  end  of  the  cylinder 


Table  M-1 1.  Measurements  on  the  type  1  HF-1  steel  cylinder  foments  with  heat  treatment  A 
from  experiment  4. 


Fragment 

Fragment 

No.  (p)  (|) 


Average  Maximum  Maximum 

_ 'V*"i  l^””***  P'**™b"<  Fn*ntll  Eofrxfnt  Trot  H«nom«». 

UkkMu  Tbkknea  with  Wki*b  Unphb  R*1U1  R»dW  n»tk 

(n»m>  Rritde  Frtcfure0  (mm)  (mm)  Stnll,t>  g^b  Str^b 


■w--  - 
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Table  M-12.  Measurements  on  the  type  1  HF-1  steel  cylinder  fragments  with  heat  treatment  B 
from  experiment  5. 


Fragment 

No. 

Fragment 

Mui 

( V ) 

(1) 

Wan 

•ntkkjw** 

(mm) 

Average 

Fraction  of  Wall 

Thickness  with 
Brittle  Fracture* 

Maximum 

Fragment 

Width** 

(mm) 

Kfaximum 

Fragment 

Length** 

(mm) 

Engineering 

Radial 

Strain'* 

True 

Radial 

St«Inb 

Hofflofeneoui 

Rank 

Slt.l<ib 

1 

554 

3.60 

19.2 

0  86 

31 

64 

16.3 

013 

014 

0  24 

2 

38  3 

3.78 

19.8 

081 

7.4 

7.8 

104 

0.10 

0.11 

019 

3 

60.2 

3.90 

18.7 

090 

58 

58 

92 

0.15 

0.17 

0.29 

4 

62.7 

4.06 

19  2 

0  89 

8  1 

6.7 

88 

0.13 

0.14 

024 

5 

62.7 

4  07 

19.1 

084 

36 

57 

17.0 

0.13 

0.14 

0.25 

6 

80  3 

3  20 

19.9 

0  97 

6.7 

5.2 

11.1 

0.10 

010 

0.18 

7 

85  9 

5.57 

19.1 

0.74 

56 

7.9 

154 

014 

015 

0  25 

949 

615 

184 

0  85 

4.3 

10.7 

16.7 

0.16 

0.18 

0.30 

9 

95.2 

6.17 

204 

0.71 

5.1 

129 

11.8 

007 

007 

0.13 

10 

99  3 

6.43 

194 

068 

5.2 

106 

15.4 

0.12 

0.12 

021 

11 

100.4 

6.51 

19.3 

049 

50 

86 

21.3 

0.12 

013 

023 

12' 

103  8 

6.73 

18.7 

0  87 

6.9 

68 

188 

0.15 

0.17 

029 

13 

106.1 

6  88 

194 

0.83 

94 

94 

124 

0.12 

0.12 

0  21 

Id' 

108.1 

7X11 

18.3 

0.83 

5.2 

8.7 

20.1 

0.17 

0.18 

0  32 

IS' 

115.7 

7.50 

18.2 

0  95 

5.4 

5.7 

180 

0.18 

0  20 

0.34 

16 

121.7 

7.89 

19.3 

0.76 

46 

9.0 

220 

0.12 

0.13 

0.23 

17 

124.7 

8j08 

194 

0.80 

7.4 

13.1 

11.2 

0.12 

0.13 

021 

18 

125.6 

8.14 

19.4 

0  70 

60 

7.6 

19.3 

012 

0.13 

0.22 

19c 

125.8 

8.15 

18  4 

092 

11.1 

59 

116 

016 

0.18 

0  31 

20c 

128  8 

8.35 

184 

0  82 

71 

7.7 

21.8 

0.17 

0.18 

0  32 

21 

131.6 

843 

194 

068 

56 

9.7 

186 

0.12 

0.12 

0.21 

22' 

135.4 

8.77 

183 

095 

9.4 

5.4 

14  8 

0.17 

0.19 

0.33 

23 

136.1 

842 

182 

0  82 

8.1 

8.1 

20.3 

0.17 

0.19 

033 

24' 

137.8 

8.9  3 

186 

085 

6.2 

9.2 

20.7 

0.16 

0J7 

0  30 

25 

1406 

9  Ii 

19  2 

084 

7.2 

86 

21.2 

013 

0.14 

0.24 

26' 

142.5 

9.23 

18.7 

084 

10.7 

88 

202 

0.15 

0.17 

029 

27 

143.4 

9.29 

19.1 

074 

4.3 

100 

26.5 

014 

0.15 

025 

28 

143.7 

9.31 

189 

0.73 

54 

85 

22.9 

0.14 

0.15 

027 

2,c 

143.9 

9,32 

18.7 

0.85 

10.0 

8.2 

19,1 

0.15 

016 

0.28 

30 

147X1 

943 

19.2 

073 

3.3 

13.3 

244 

0.13 

0.14 

024 

si' 

147.8 

948 

186 

0  85 

88 

6.1 

200 

0.16 

0.17 

0.30 

32' 

1514 

9  81 

18  3 

0.79 

11.7 

10.7 

1 6  6 

0.17 

0.19 

032 

33' 

155.1 

10.05 

18  8 

088 

11.3 

108 

194 

0.15 

0.16 

0.28 

34 

155.6 

10X18 

193 

0.79 

68 

244 

0.13 

013 

023 

!O.J 
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Table  M-12.  Measurements  on  the  type  1  HF-1  steel  cylinder  fragments  with  heat  treatment  B 
from  experiment  5.(Continued) 


Fragment 

No. 

Fragment 

Mass 

<W 

(8) 

Wall 

Thickness* 

(mm) 

Average 

Fraction  of  Wall 
Thickness  with 
Brittle  Fracture^ 

Maximum 

Fragment 

Width** 

(mm) 

Maximum 

fragment 

Length** 

(mm) 

Engineering 

Radial 

Strain** 

True 

Radial 

Strtinb 

Homogeneous 

Plastic 

Str»lnb 

3  Sc 

157.2 

10.19 

184 

0.87 

7.9 

8.9 

23.8 

0.17 

0.18 

0.31 

36 

159.3 

10.32 

19.2 

0.64 

6.1 

11.1 

20.8 

0.13 

0.14 

0.25 

37 

165.1 

10.70 

19.3 

0.78 

50 

10.2 

22.0 

0.13 

0.13 

0.23 

38 

166.6 

10.80 

18.4 

0.77 

6.3 

8.8 

26.7 

0.16 

0.18 

0.31 

39 

169.4 

10.98 

19.1 

0.72 

8.8 

85 

21.9 

0.14 

C  IS 

0.25 

40 

169-5 

10.98 

19.6 

0.78 

10.1 

11.4 

19.3 

0.11 

0.12 

0.21 

41 

171.5 

11.11 

19.3 

0.65 

45 

8.3 

24.2 

0.13 

0.14 

0.24 

42 

172.3 

11.17 

19.6 

0.78 

6.7 

9.0 

24.9 

0.11 

0.12 

0.21 

43 

177.9 

1157 

19.1 

0.65 

7.8 

10.6 

19.1 

0.14 

0.15 

0.25 

44 

1785 

11.57 

19.1 

0.79 

5.2 

95 

22.0 

0.13 

0.14 

0.25 

45C 

188.8 

12.23 

18.4 

0.81 

10.4 

9.6 

21.0 

0.16 

0.18 

0.31 

46C 

192.6 

12.48 

18.4 

0.83 

9.9 

14.4 

19.0 

0.17 

0.18 

0.32 

47C 

199.8 

12.95 

18.5 

0.74 

11.2 

9.1 

17.6 

0.16 

0.17 

0.30 

48C 

206.1 

13.36 

18.7 

0.86 

8.3 

12.0 

21.6 

0.15 

0.17 

0.29 

49 

211.4 

13.70 

19.3 

0.65 

5.6 

11.8 

29.6 

0.13 

0.13 

0.22 

50 

212.9 

13.80 

19.7 

0.76 

6.4 

12.6 

25.3 

0.11 

0.12 

0.20 

51 

219.1 

14.20 

195 

0.85 

5.4 

11.0 

22.8 

0.12 

0.13 

0.22 

52c 

232.8 

15.09 

19.3 

0.83 

12.9 

9.3 

21.2 

0.12 

0.13 

0.23 

53c 

237.0 

15.36 

18.7 

0.86 

11.6 

9.1 

19.9 

0.15 

0.17 

0.29 

54' 

240.4 

15.58 

18.4 

0.85 

10.1 

11.4 

19.8 

0.16 

0.18 

0.31 

55 

277.4 

17.98 

19.2 

0.72 

7.7 

135 

29.4 

0.13 

0.14 

0.24 

*  The  initial  wall  thickness  was  22.1  mm. 
b  Defined  in  footnotes  in  Table  M-8 
c  This  fragment  came  from  the  end  of  the  cylinder. 
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Table  M-13.  Measurements  on  the  type  1  HF-1  steel  cylinder  fragments  with  heat  treatment  B 
from  experiment  6. 


Fragment 

No. 

Fuimtnt 

Mata 

(p) 

(*> 

Well 

Thickness* 

(mm) 

Average 

Fraction  of  Well 
Thickness  with 
Brittle  Fracture*5 

Maximum 

Fragment 

Width*1 

(mm) 

Maximum 

Fragment 

Length*5 

(mm) 

Engineering 

Radial 

StHln11 

True 

Radial 

Strain*5 

Homogeneous 

Plastic 

Strain*5 

ic 

S5.I 

3.57 

19.5 

0,76 

5.4 

9.1 

17.0 

012 

0.13 

023 

2 

61.1 

3.96 

194 

0  80 

3  0 

5.3 

18.2 

013 

0.14 

0  23 

»c 

77,9 

5  05 

IM 

0.77 

50 

7,6 

18 .5 

012 

0.13 

023 

4 

101.4 

6.57 

19.5 

0.76 

5.4 

9.1 

17.0 

0.12 

0.13 

0  23 

5 

107,3 

6.95 

19.7 

0.67 

5.1 

11.4 

22.5 

0.11 

0.12 

021 

6 

1106 

7,17 

19.3 

0.87 

5.5 

8.0 

21.4 

0.13 

0.14 

0.25 

7 

1113 

7.21 

19  6 

0  82 

5.3 

7.7 

23.5 

0.12 

0.13 

0.22 

8 

112  1 

7.26 

19.4 

0.76 

5.0 

10.3 

18.5 

0.13 

0.14 

024 

9 

I1S.2 

7.47 

19.8 

0  81 

6.6 

11.0 

12.3 

0.11 

0.12 

0.20 

10 

122.9 

7.96 

19  6 

0  82 

3.2 

8.8 

22.7 

0.12 

0.13 

0.22 

II 

125.6 

8.14 

19  7 

0.84 

8.3 

11.7 

17.6 

0.11 

0.12 

0.21 

12 

126  6 

8.20 

19.8 

061 

5.5 

10  2 

22.4 

0.11 

0.11 

020 

13 

J28.I 

8.30 

196 

069 

5.1 

106 

17.4 

0.12 

0.12 

0  21 

14C 

128.2 

8.31 

19.2 

0.90 

100 

90 

19.7 

0.13 

0.14 

0.25 

IS 

130.7 

8  47 

19.2 

0  82 

6.6 

11.2 

19.1 

0.13 

0.14 

02$ 

16 

133  6 

8  66 

19.3 

0  83 

6.5 

10.4 

15.3 

0.13 

0.14 

0.24 

17C 

136.4 

884 

189 

0.87 

7.2 

9.8 

20.9 

0.15 

0.16 

0.28 

18c 

136.7 

8.86 

19.2 

080 

6.9 

109 

20.5 

0.14 

0.15 

025 

I9C 

137.6 

8  92 

20.1 

0.85 

10.5 

10.4 

23.1 

009 

0.10 

0.17 

20c 

138.2 

8  96 

188 

0.86 

10  0 

7.0 

21.9 

0.15 

0.17 

029 

21 

141.2 

9.15 

19.7 

0.59 

4.6 

12.8 

184 

0.11 

0.12 

0.21 

22 

142.2 

9.21 

19  6 

0  85 

6.2 

10.1 

206 

0.12 

0.13 

022 

23c 

145.1 

9.40 

19.4 

084 

10.1 

8.1 

21.0 

0.13 

0.14 

0  23 

24 

148.5 

9.62 

200 

0  80 

11.9 

90 

15.1 

0.10 

0.11 

0.18 

25c 

149.3 

9.67 

19.0 

0.85 

8.3 

102 

20.8 

0.15 

0.16 

0  27 

26 

151.3 

9.80 

19.5 

0  80 

6.2 

11,3 

21.5 

0.12 

0.13 

022 

27* 

154.3 

1000 

19.6 

0  82 

10.9 

9.1 

16.6 

0.12 

0.13 

0.22 

nc 

157^4 

10.20 

19.5 

0  82 

8.3 

80 

24.5 

009 

0.10 

0.17 

29 

164.1 

10.63 

19.7 

0.77 

S3 

1 1.7 

20.6 

0.11 

0.12 

021 

30 

166.2 

10.77 

19.2 

0.77 

7.1 

8.2 

23  9 

0.14 

0.15 

026 

31c 

167,4 

1085 

19.7 

0.80 

7.6 

11.5 

19.9 

0.1 1 

0.12 

0.21 
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Table  M-13.  Measurements  on  the  type  1  HF-1  steel  cylinder  fragments  with  heat  treatment  B 
from  experiment  6.(Continued) 


Fragment 

No. 

Fragment 

Mass 

(gr) 

(8) 

Wall 

Thickness* 

(mm) 

Average 

Fraction  of  Wall 
Thickness  with 
Brittle  Fracture** 

Maximum 

Fragment 

Width6 

(mm) 

Maximum 

Fragment 

Length** 

(mm) 

Engineering 

Radial 

Strain** 

True 

Radial 

Strain6 

Homogeneous 

Plastic 

Strain** 

32c 

169.4 

10.98 

19.4 

0  89 

8.3 

9.0 

23.1 

0.13 

0.14 

0.23 

33 

174.2 

11.29 

19.7 

0.78 

6.9 

10.9 

19.6 

0.11 

0.12 

0.21 

34c 

178.0 

11.53 

19.6 

0.82 

8.1 

8.8 

21.2 

0.12 

0.13 

0.22 

35c 

180.8 

11.72 

19.4 

0.89 

8.8 

1 1.1 

18.7 

0.13 

0.14 

0.23 

36 

181.4 

11.75 

19.3 

0  81 

6.5 

9.9 

22.7 

0.13 

0.14 

0.25 

37 

184.3 

11.94 

19.8 

0.78 

11.2 

6.2 

22.5 

0.11 

0.12 

0.20 

38 

18S.4 

12.01 

19.4 

0.75 

9.7 

7.4 

24,2 

0.13 

0.14 

0.24 

39c 

190.1 

12.32 

19.3 

0.80 

7.2 

8.6 

25.9 

0.13 

0.14 

024 

40c 

194.2 

12.58 

19.9 

0.81 

11.2 

10.4 

24.3 

0.11 

0.1 1 

0.19 

41C 

207.6 

13.45 

20.1 

0-78 

10.1 

14.6 

23.4 

0.09 

0.10 

0.17 

42c 

208.2 

13.49 

19.9 

0.86 

9.4 

12.3 

22.6 

0.11 

0.1 1 

0.19 

43c 

220.1 

14.26 

19.4 

0.87 

13.2 

12.7 

22.6 

0.13 

0.14 

0.24 

44c 

234.7 

15.21 

19.5 

0.89 

14.1 

9.5 

28.1 

0.12 

0.13 

0.23 

45 

247.5 

16.04 

19.7 

0.72 

9.3 

13.2 

21.8 

0.11 

0.12 

0.21 

46c 

248.5 

16.10 

19.4 

0.90 

1S.9 

9.3 

23.0 

0.13 

0.14 

0.23 

47c 

291.8 

18.91 

19.8 

0,82 

11.6 

10“ 

26.7 

0.11 

0.12 

0.20 

48 

319.2 

20  68 

19.7 

0.85 

6.7 

105 

32.1 

0.11 

0.12 

0.21 

49 

328.1 

21.26 

19.8 

0.67 

7.7 

13. 3 

32.4 

0.11 

0.12 

0.20 

*  The  Initial  wall  thickness  was  22.3  mm. 

**  Defined  In  footnotes  In  Table  M*8. 
c  ThU  figment  came  from  the  end  of  the  cylinder. 
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twenty  selected  type 
i  experiment  7. 


Fragment 

No. 

Fragment 

Maas 

(P) 

(1) 

Wall 

Thickness* * 

(mm) 

Average 

Fraction  of  Wall 
Thickness  with 
Brittle  Fracture* 

Maximum 

Fragment 

Width* 

(mm) 

Maximum 

Fragment 

Length* 

(mm) 

Engineering 

Radial 

Striin* 

True 

Radial 

Strain* 

Homogeneous 

Mastic 

Strain* 

1  ,b 

794.3 

51.47 

16.0 

0.56 

16.5 

58.8 

0.11 

0.12 

0.20 

, 

16.3 

■  2 

639.6 

41.44 

15.9 

0.43 

8.2 

67.2 

0.12 

0.12 

0.22 

5.8 

3 

516.8 

33.49 

15.7 

0.47 

7.7 

58.2 

0.13 

0.14 

0.24 

4  466.9 

5  428.4 

6  348.3 

7  298.6 

8  204.9 

9  153.2 

10  104.2 

Uc  152.3 

12  102.9 


16“  202.9 

17  74.71 


*  Defined  in  footnotes  In  Table  M*8. 

b  Fragments  I  through  10  came  from  the  102-mm  long  uniform  region  of  the  aft  projectile  section.  The  initial  wall  thickness  was  18.0  mm, 

*  Fragments  1 1  through  15  came  from  the  joint  region  of  the  forward  projectile  section.  The  Initial  wall  thickness  was  9.7  mm, 
d  Fragments  16  through  20  came  from  the  joint  region  of  the  aft  projectile  section.  Tho  initial  wall  thickness  was  8.7  mm. 


APPENDIX  N 

COMPUTER  LISTINGS  FOR  THE  ARMCO  IRON  AND  HF-1  STEEL 
CYLINDRICAL  PUFF  AND  TROTT  COMPUTATIONS 


(Computer  listings  of  the  PUFF  computational  results  for  the  Arrnco  iron  and  HF-1  steel  cylinders 
and  the  TROTT  computational  results  for  the  HF-1  steel  projectile  are  given  in  this  appendix.  The 
BFRACT2  subroutine  was  used  for  the  computations  given  in  Figures  N-l,  N-2,  N-3,  and  N-6.  The 
BFRACT2  and  SHEAR2  subroutines  were  used  for  the  computations  given  in  Figures  N-4  and 
N-S.  References  1  and  2  provide  detailed  explanations  of  the  parameters  used  in  the  PUFF  and 
TROTT  computer  codes,  respectively.  The  crack  distribution  results  for  these  computations 
are  summarized  in  Appendix  0.  The  fragment  mass  distribution  results  that  are  listed  in  the 
computer  output  for  each  brittle  fracture  computational  cell  were  derived  assuming  spherical 
fragments  and  were  not  used  in  the  present  work.  The  details  are  given  in  Reference  3.) 
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PUFF  computation  using  the  3FRACT2  subroutine. 


Figure  N-1  (Continued) 
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Figure  N-1  (Continued) 
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Figure  N-2  (Continued) 
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Figure  N-2  (Continued) 
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Figure  N-3.  Computer  listing  of  (a)  input  parameters  and  (b)  output  results  at  500  time  cycles  for  the  HF-1  steel,  heat  treat¬ 
ment  B,  cylindrical  PUFF  computation  using  the  BFRACT2  subroutine. 
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Figure  N-3  (Continued) 
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Figure  N-3  (Continued) 


Figure  N-4.  Computer  listing  of  (a)  input  parameters  and  (b)  output  results  at  500  time  cycles  for  the  HF-1  steel,  heat  treat¬ 
ment  A,  cylindrical  PUFF  computation  using  the  BFRACT2  and  SHEAR2  subroutines. 
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Figure  N-4  (Continued) 
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Figure  N-5  (Continued) 
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Figure  N-6  (Continued) 
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Figure  N-6  (Continued) 
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Figure  N-6  (Continued) 
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Figure  N-6  (Continued) 
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Figure  N-6  (Continued) 
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Figure  N-6  (Continued) 
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Figure  N-6  (Continued) 


APPENDIX  O 


SUMMARY  OF  THE  COMPUTATIONAL  CRACK  DISTRIBUTION  RESULTS  AND  THE 
CALCULATED  FRAGMENT  MASS  DISTRIBUTION  RESULTS  FOR  THE  ARMCO  IRON  AND 
HF-1  STEEL  CYLINDER  EXPERIMENTS  AND  THE  HF-1  STEEL  PROJECTILE  EXPERIMENT 


(Tables  0-1,  0-3,  0-5,  0-7,  0-9,  and  0-1 1  summarize  the  crack  distribution  results  for  the  PUFF 
and  TROTT  computations  given  in  Appendix  N.  Tables  0-2,  0-4,  0-6,  0-8,  0-10,  and  0-11  give 
the  calculated  fragment  mass  distribution  results  using  the  crack  distribution  results.) 


■’nP  ■-*  <r  ?~ 


Table  0-1.  Summary  of  brittle  fracture  crack  distribution  results  for  the  Armco  iron  cylindrical  PUFF  computation. 


Summary  of  brittle  fracture  crack  distribution  results  for  the  HF-1  steel,  heat  treatment  A,  cylindrical  PUFF  computation 


Table  O-o.  Summary  of  brittle  fracture  crack  distribution  results  for  the  HF-1  steel,  heat  treatment  B,  cylindrical  PUFF  computation.1 
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Table  0-6.  Calculated  fragment  mass  distribution  results  for  the  HF-1  steel  cylinder  with  heat  treatment  8  using  the  brittle  fracture 
crack  distribution  results  from  Table  0-5.* 


Table  0-7.  Summary  of  brittle  and  shear  fracture  crack  distribution  results  for  the  HF-1  steel,  heat  treatment  A,  cylindrical  PUFF  computation.' 
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APPENDIX  P 

COMPUTED  STRESS,  PARTICLE  VELOCITY,  AND  CRACK 
CONCENTRATION  VERSUS  TIME  PLOTS  FOR  EACH 
COMPUTATIONAL  CELL  FOR  THE  HF-1  STEEL, 

HEAT  TREATMENT  B,  CYLINDER  EXPERIMENTS 

(These  plots  were  produced  using  cylindrical  PUFF  computations  with  the  BFRACT  2  subroutine 
and  the  input  parameters  from  Figure  N-3(a)  (with  the  appropriate  JEDIT  parameters).  The 
JEDIT  parameters  SI,  S2,  S3,  U,  COM9,  and  COM1 1  used  for  these  plots  correspond  to  the  radial 
stress,  circumferential  stress,  axial  stress,  radial  particle  velocity,  circumferential  crack  concentr¬ 
ation,  and  axial  crack  concentration,  respectively.  Cells  1  through  10  contain  composition  B 
explosive.  Cells  12  through  21  contain  HF-1  steel  with  heat  treatment  B.  Plots  (a),  (b),  and  (c) 
of  Figures  P-1  through  P-20  are  the  principal  stresses  in  the  radial,  circumferential,  and  axial 
directions,  respectively.  Compressive  and  tensile  stress  values  are  positive  and  negative, 
respectively.  The  stresses  are  computed  at  the  cell  midpoints.  Figures  P-21  through  P-30  give  the 
radial  particle  velocity  in  the  cells.  A  particle  velocity  refers  to  the  cell  edge.  Plots  (a)  and  (b)  of 
Figures  31  through  40  give  the  circumferential  and  axial  crack  concentrations,  respectively,  for  the 
fractured  HF-1  steel  cells.  Crack  concentrations  are  computed  at  the  cell  midpoints). 


P-1 


Computed  stress  histories  in  cell  5  for  the  HF-1  Figure  P-6.  Computed  stress  histories  in  cell  6  for  the  HF-1 
steel,  heat  treatment  B,  cylinder  experiments.  steel,  heat  treatment  B,  cylinder  experiments, 

(a)  Radial,  (b)  circumferential,  and  (c)  axial  (a)  Radial,  (b)  circumferential,  and  (c)  axial 

stresses.  stresses. 
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i  P-13.  Computed  stress  histories  in  cell  14  for  the  HF-1  Figure  P-14.  Computed  stress  histories  in  cell  15  for  the  HF-1 
steel,  heat  treatment  B,  cylinder  experiments.  steel,  heat  treatment  B,  cylinder  experiments, 

(a)  Radial,  (b)  circumferential,  and  (c)  axial  (a)  Radial,  (b)  circumferential,  and  (c)  axial 

stresses.  stresses. 


Figure  P-19.  Computed  stress  histories  in  ceil  20  for  the  HF-1  Figure  P-20.  Computed  stress  histories  in  cell  21  for  the  HF-1 
steel,  heat  treatment  B,  cylinder  experiments.  steel,  heat  treatment  B,  cylinder  experiments, 

(a)  Radial,  (b)  circumferential,  and  (c)  axial  (a)  Radial,  (b)  circumferential,  and  (c)  axial 

stresses.  stresses. 
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imputed  particle  velocity  histories  in  (a)  cell  9  Figure  P-26.  Computed  particle  velocity  histories  in  (a)  cel! 
d  (b)  cell  10  for  the  HF-1  steel,  heat  treatment  12  and  (b)  cell  13  for  the  HF-1  steel,  heat  treat- 

cylinder  experiments.  ment  B,  cylinder  experiments. 


Computed  particle  velocity  histories  in  (a)  cell  Figure  P-28.  Computed  particle  velocity  histories  in  (a)  cell 
14  and  (b)  cell  15  for  the  HF-1  steel,  heat  16  and  (b)  cell  17  for  the  HF-1  steel,  heat 

treatment  B,  cylinder  experiments.  treatment  B,  cylinder  experiments. 
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gure  P-33.  Computed  crack  concentration  histories  in  cell  Figure  P-34.  Computed  crack  concentration  histories  in  cell 
14  for  the  HF-1  steel,  heat  treatment  B,  cylinder  15  for  the  HF-1  steel,  heat  treatment  B,  cylinder 

experiments,  (a)  Circumferential  and  (b)  axial  experiments,  (a)  Circumferential  and  (b)  axial 

crack  concentrations.  crack  concentrations. 
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Computed  crack  concentration  histories  in  cell  Figure  P-36.  Computed  crack  concentration  histories  in  cell 
16  for  the  HF-1  steel,  heat  treatent  B,  cylinder  17  for  the  HF-1  steel,  heat  treatment  B,  cylinder 

experiments,  (a)  Circumferential  and  (b)  axial  experiments,  (a)  Circumferential  and  (b)  axial 

crack  concentrations.  crack  concentrations. 
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Computed  crack  concentration  histories  in  cell  Figure  P-38.  Computed  crack  concentration  histories  in  cell 
18  for  the  HF-1  steel,  heat  treatment  B,  cylinder  19  for  the  HF-1  steel,  heat  treatment  B,  cylinder 

experiments,  (a)  Circumferential  and  (b)  axial  experiments,  (a)  Circumferential  and  (b)  axial 

crack  concentrations.  crack  concentrations. 
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